diff options
| author |  Edoardo Zoni <59625522+EZoni@users.noreply.github.com>
| 2020-07-06 11:26:41 -0700 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2020-07-06 11:26:41 -0700 |
| commit | 345feb7faa0647ec52025adb450c2855154e8111 (patch) | |
| tree | 10fc5525a86b9aacfef9ef537d0ee363ec505c55 /Source/WarpX.cpp | |
| parent | 0a146f12b7a18159561bc74595a3853063216d3c (diff) | |
| download | WarpX-345feb7faa0647ec52025adb450c2855154e8111.tar.gz WarpX-345feb7faa0647ec52025adb450c2855154e8111.tar.zst WarpX-345feb7faa0647ec52025adb450c2855154e8111.zip | |
PSATD: add option to update E without using rho (#1128)
* Introduce option to update E with/without rho
* Clean up
* Include equations in docs
* Fix EOL whitespaces error
* Small clean-up
* Clean up
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( |