diff options
| author |  Remi Lehe <remi.lehe@normalesup.org>
| 2019-04-30 09:41:21 -0700 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2019-04-30 09:41:21 -0700 |
| commit | bb1289f3879155244372f791d6854107150f4266 (patch) | |
| tree | 5c2697eeb93fa75b6f82b5b3603c393b6c57e39a /Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp | |
| parent | 9a7a87fcd64d7faab8c57cbc769840f8ee74454a (diff) | |
| parent | fcab5379ddadb02f81d14c9583000c0aeb5b7db6 (diff) | |
| download | WarpX-bb1289f3879155244372f791d6854107150f4266.tar.gz WarpX-bb1289f3879155244372f791d6854107150f4266.tar.zst WarpX-bb1289f3879155244372f791d6854107150f4266.zip | |
Merge pull request #104 from ECP-WarpX/fix_spectral
Fix bugs in the spectral solver
Diffstat (limited to 'Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp')
| -rw-r--r-- | Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp | 8 |
1 files changed, 6 insertions, 2 deletions
diff --git a/Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp b/Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp index 844ac97a7..6e6cc124f 100644 --- a/Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp +++ b/Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp @@ -159,8 +159,10 @@ SpectralFieldData::ForwardTransform( const MultiFab& mf, if (is_nodal_x==false) spectral_field_value *= xshift_arr[i]; #if (AMREX_SPACEDIM == 3) if (is_nodal_y==false) spectral_field_value *= yshift_arr[j]; -#endif if (is_nodal_z==false) spectral_field_value *= zshift_arr[k]; +#elif (AMREX_SPACEDIM == 2) + if (is_nodal_z==false) spectral_field_value *= zshift_arr[j]; +#endif // Copy field into the right index fields_arr(i,j,k,field_index) = spectral_field_value; }); @@ -212,8 +214,10 @@ SpectralFieldData::BackwardTransform( MultiFab& mf, if (is_nodal_x==false) spectral_field_value *= xshift_arr[i]; #if (AMREX_SPACEDIM == 3) if (is_nodal_y==false) spectral_field_value *= yshift_arr[j]; -#endif if (is_nodal_z==false) spectral_field_value *= zshift_arr[k]; +#elif (AMREX_SPACEDIM == 2) + if (is_nodal_z==false) spectral_field_value *= zshift_arr[j]; +#endif // Copy field into temporary array (after normalization) tmp_arr(i,j,k) = inv_N*spectral_field_value; }); |