[DEBUG] additional condition if nx is even

QuentinMale · QuentinMale · commit 6f1d0cf80d38 · 2025-03-28T08:58:06.000-04:00
diff --git a/ExampleCodes/FFT/EnergySpectrum/main.cpp b/ExampleCodes/FFT/EnergySpectrum/main.cpp
@@ -59,8 +59,9 @@ int main (int argc, char* argv[])
                     +        amrex::norm(cza[b](i,j,k));
 		// Account for Hermitian symmetry in x-direction
 	        // Hermitian symmetry Y[nx-i,j,k] = Y[i,j,k]*
-                if (i > 0) { // Avoid double-counting at kx = 0
-                    value *= 2.0; // Multiply by 2 because we have kx and -kx
+                if ((i > 0) && (2*i != nx)) {
+		    // Multiply by 2 because we have +ki and -ki
+                    value *= 2.0;
                 }
                 HostDevice::Atomic::Add(pke+di, value);
             }

Original file line number	Diff line number	Diff line change
`@@ -59,8 +59,9 @@ int main (int argc, char* argv[])`
`59`	`59`	`+ amrex::norm(cza[b](i,j,k));`
`60`	`60`	`// Account for Hermitian symmetry in x-direction`
`61`	`61`	`// Hermitian symmetry Y[nx-i,j,k] = Y[i,j,k]*`
`62`		`- if (i > 0) { // Avoid double-counting at kx = 0`
`63`		`- value *= 2.0; // Multiply by 2 because we have kx and -kx`
	`62`	`+ if ((i > 0) && (2*i != nx)) {`
	`63`	`+ // Multiply by 2 because we have +ki and -ki`
	`64`	`+ value *= 2.0;`
`64`	`65`	`}`
`65`	`66`	`HostDevice::Atomic::Add(pke+di, value);`
`66`	`67`	`}`