1 files changed, 7 insertions, 8 deletions
diff --git a/Source/FieldSolver/SpectralSolver/PsatdAlgorithm.cpp b/Source/FieldSolver/SpectralSolver/PsatdAlgorithm.cpp
index 60e9d58c0..900b542be 100644
--- a/Source/FieldSolver/SpectralSolver/PsatdAlgorithm.cpp
+++ b/Source/FieldSolver/SpectralSolver/PsatdAlgorithm.cpp
@@ -9,17 +9,16 @@ PsatdAlgorithm::PsatdAlgorithm(const SpectralKSpace& spectral_kspace,
                          const DistributionMapping& dm,
                          const int norder_x, const int norder_y,
                          const int norder_z, const bool nodal, const Real dt)
-{
-    const BoxArray& ba = spectral_kspace.spectralspace_ba;
-
-    // Allocate the 1D vectors
-    modified_kx_vec = spectral_kspace.getModifiedKComponent(dm, 0, norder_x, nodal);
+// Compute and assign the modified k vectors
+: modified_kx_vec(spectral_kspace.getModifiedKComponent(dm, 0, norder_x, nodal))
 #if (AMREX_SPACEDIM==3)
-    modified_ky_vec = spectral_kspace.getModifiedKComponent(dm, 1, norder_y, nodal);
-    modified_kz_vec = spectral_kspace.getModifiedKComponent(dm, 2, norder_z, nodal);
+  modified_ky_vec(spectral_kspace.getModifiedKComponent(dm, 1, norder_y, nodal))
+  modified_kz_vec(spectral_kspace.getModifiedKComponent(dm, 2, norder_z, nodal))
 #else
-    modified_kz_vec = spectral_kspace.getModifiedKComponent(dm, 1, norder_z, nodal);
+  modified_kz_vec(spectral_kspace.getModifiedKComponent(dm, 1, norder_z, nodal))
 #endif
+{
+    const BoxArray& ba = spectral_kspace.spectralspace_ba;
 
     // Allocate the arrays of coefficients
     C_coef = SpectralCoefficients(ba, dm, 1, 0);