diff options
| author |  Edoardo Zoni <59625522+EZoni@users.noreply.github.com>
| 2022-09-20 16:28:39 -0700 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2022-09-20 16:28:39 -0700 |
| commit | 5761b4bf998eab84ea0f7e4b132026593f3ddf9f (patch) | |
| tree | 39a60a20ac3fac7368231b66f49fad063321f624 /Source/BoundaryConditions/PML.cpp | |
| parent | 2fed2828933831ee464f0ca5d02a23dd2df54aad (diff) | |
| download | WarpX-5761b4bf998eab84ea0f7e4b132026593f3ddf9f.tar.gz WarpX-5761b4bf998eab84ea0f7e4b132026593f3ddf9f.tar.zst WarpX-5761b4bf998eab84ea0f7e4b132026593f3ddf9f.zip | |
PSATD: More Options for Time Dependency of J, Rho (#3242)
* Rename PsatdAlgorithm as PsatdAlgorithmJConstantInTime
* Add New Inputs: psatd.J_in_time, psatd.rho_in_time
* Update PSATD Classes
* Fix Bug for FDTD Build
* Fix Warning for RZ PSATD Build
* Fix Multi-J Loop w/ J Constant in Time
* Clean up
* Fix Error Message
* Fix Time of Charge Deposition for Rho Linear
* Add 3D Langmuir CI Tests w/ Multi-J Algo
* Fix Checksums of New CI Tests
* Remove Extra CI Tests (added in #3363)
* Fix CI Tests w/ Multi-J PSATD
* Add Docs for New Input Parameters
* Galilean/Comoving PSATD Not Compatible w/ J Linear
Diffstat (limited to 'Source/BoundaryConditions/PML.cpp')
| -rw-r--r-- | Source/BoundaryConditions/PML.cpp | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/Source/BoundaryConditions/PML.cpp b/Source/BoundaryConditions/PML.cpp index 9c3c17b0b..712e287f0 100644 --- a/Source/BoundaryConditions/PML.cpp +++ b/Source/BoundaryConditions/PML.cpp @@ -548,7 +548,7 @@ PML::PML (const int lev, const BoxArray& grid_ba, const DistributionMapping& gri int ncell, int delta, amrex::IntVect ref_ratio, Real dt, int nox_fft, int noy_fft, int noz_fft, bool do_nodal, int do_moving_window, int /*pml_has_particles*/, int do_pml_in_domain, - const bool do_multi_J, + const int J_in_time, const int rho_in_time, const bool do_pml_dive_cleaning, const bool do_pml_divb_cleaning, const amrex::IntVect& fill_guards_fields, const amrex::IntVect& fill_guards_current, @@ -738,7 +738,7 @@ PML::PML (const int lev, const BoxArray& grid_ba, const DistributionMapping& gri if (WarpX::maxwell_solver_id == MaxwellSolverAlgo::PSATD) { #ifndef WARPX_USE_PSATD - amrex::ignore_unused(lev, dt, do_multi_J); + amrex::ignore_unused(lev, dt, J_in_time, rho_in_time); # if(AMREX_SPACEDIM!=3) amrex::ignore_unused(noy_fft); # endif @@ -759,7 +759,7 @@ PML::PML (const int lev, const BoxArray& grid_ba, const DistributionMapping& gri spectral_solver_fp = std::make_unique<SpectralSolver>(lev, realspace_ba, dm, nox_fft, noy_fft, noz_fft, do_nodal, v_galilean_zero, v_comoving_zero, dx, dt, in_pml, periodic_single_box, update_with_rho, - fft_do_time_averaging, do_multi_J, m_dive_cleaning, m_divb_cleaning); + fft_do_time_averaging, J_in_time, rho_in_time, m_dive_cleaning, m_divb_cleaning); #endif } @@ -878,7 +878,7 @@ PML::PML (const int lev, const BoxArray& grid_ba, const DistributionMapping& gri spectral_solver_cp = std::make_unique<SpectralSolver>(lev, realspace_cba, cdm, nox_fft, noy_fft, noz_fft, do_nodal, v_galilean_zero, v_comoving_zero, cdx, dt, in_pml, periodic_single_box, update_with_rho, - fft_do_time_averaging, do_multi_J, m_dive_cleaning, m_divb_cleaning); + fft_do_time_averaging, J_in_time, rho_in_time, m_dive_cleaning, m_divb_cleaning); #endif } } |