diff options
Diffstat (limited to 'Source/WarpX.cpp')
| -rw-r--r-- | Source/WarpX.cpp | 6 |
1 files changed, 4 insertions, 2 deletions
diff --git a/Source/WarpX.cpp b/Source/WarpX.cpp index 9999dd40b..22debe82e 100644 --- a/Source/WarpX.cpp +++ b/Source/WarpX.cpp @@ -606,6 +606,7 @@ WarpX::ReadParameters () pp.query("noy", noy_fft); pp.query("noz", noz_fft); pp.query("current_correction", current_correction); + pp.query("update_with_rho", update_with_rho); pp.query("v_galilean", v_galilean); // Scale the velocity by the speed of light for (int i=0; i<3; i++) v_galilean[i] *= PhysConst::c; @@ -906,7 +907,7 @@ WarpX::AllocLevelMFs (int lev, const BoxArray& ba, const DistributionMapping& dm bool const pml=false; spectral_solver_fp[lev].reset( new SpectralSolver( realspace_ba, dm, nox_fft, noy_fft, noz_fft, do_nodal, v_galilean, dx_vect, dt[lev], - pml, fft_periodic_single_box ) ); + pml, fft_periodic_single_box, update_with_rho ) ); # endif #endif m_fdtd_solver_fp[lev].reset( @@ -997,7 +998,8 @@ WarpX::AllocLevelMFs (int lev, const BoxArray& ba, const DistributionMapping& dm # else realspace_ba.grow(ngE); // add guard cells spectral_solver_cp[lev].reset( new SpectralSolver( realspace_ba, dm, - nox_fft, noy_fft, noz_fft, do_nodal, v_galilean, cdx_vect, dt[lev] ) ); + nox_fft, noy_fft, noz_fft, do_nodal, v_galilean, cdx_vect, dt[lev], + pml, fft_periodic_single_box, update_with_rho ) ); # endif #endif m_fdtd_solver_cp[lev].reset( |