From 2cec61478b9bc33d094b6c8f0ff1f968a149266c Mon Sep 17 00:00:00 2001 From: Remi Lehe Date: Mon, 14 Feb 2022 16:41:52 -0800 Subject: 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 --- Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp') 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; } } -- cgit v1.2.3