diff options
author | 2022-02-14 16:41:52 -0800 | |
---|---|---|
committer | 2022-02-15 00:41:52 +0000 | |
commit | 2cec61478b9bc33d094b6c8f0ff1f968a149266c (patch) | |
tree | 9cc45afc35a475ee892218f2b44a8e836b7916f4 | |
parent | e39f6e560214152ff279c296505f08b98d438011 (diff) | |
download | WarpX-2cec61478b9bc33d094b6c8f0ff1f968a149266c.tar.gz WarpX-2cec61478b9bc33d094b6c8f0ff1f968a149266c.tar.zst WarpX-2cec61478b9bc33d094b6c8f0ff1f968a149266c.zip |
Only set modified k to 0 for even number of points (#2852)
* Only set modified k to 0 for even number of points
* Update Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp
-rw-r--r-- | Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp b/Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp index 2f11f84ca..6d31116b6 100644 --- a/Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp +++ b/Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp @@ -266,7 +266,7 @@ SpectralKSpace::getModifiedKComponent( const DistributionMapping& dm, } else { // The other axes contains both positive and negative k ; // the Nyquist frequency is in the middle of the array. - if (i == N/2) { + if ( (N%2==0) && (i == N/2) ){ p_modified_k[i] = 0.0_rt; } } |