diff options
author | 2017-12-04 15:42:20 -0800 | |
---|---|---|
committer | 2017-12-04 15:42:20 -0800 | |
commit | 2cb9607db613e5f1bff48b2b0de94ed06412e3b0 (patch) | |
tree | 2b1e92decde80055fe9db0dcbd1f35ed68f65a06 /Source/WarpXPML.cpp | |
parent | eab724c5391f29f19cf75e16f1348120cf0ca0a3 (diff) | |
parent | 5a33f6890451b27bb8a2930faeea431beb58910d (diff) | |
download | WarpX-2cb9607db613e5f1bff48b2b0de94ed06412e3b0.tar.gz WarpX-2cb9607db613e5f1bff48b2b0de94ed06412e3b0.tar.zst WarpX-2cb9607db613e5f1bff48b2b0de94ed06412e3b0.zip |
Merge branch 'master' into io
Conflicts:
Source/WarpX.cpp
Diffstat (limited to 'Source/WarpXPML.cpp')
-rw-r--r-- | Source/WarpXPML.cpp | 7 |
1 files changed, 4 insertions, 3 deletions
diff --git a/Source/WarpXPML.cpp b/Source/WarpXPML.cpp index 6db887136..19e36f38a 100644 --- a/Source/WarpXPML.cpp +++ b/Source/WarpXPML.cpp @@ -425,7 +425,7 @@ MultiSigmaBox::ComputePMLFactorsE (const Real* dx, Real dt, const std::string& p PML::PML (const BoxArray& grid_ba, const DistributionMapping& grid_dm, const Geometry* geom, const Geometry* cgeom, - int ncell, int delta, int ref_ratio, int do_dive_cleaning) + int ncell, int delta, int ref_ratio, int do_dive_cleaning, int do_moving_window) : m_geom(geom), m_cgeom(cgeom) { @@ -441,6 +441,7 @@ PML::PML (const BoxArray& grid_ba, const DistributionMapping& grid_dm, int nge = 2; int ngb = 2; + int ngf = (do_moving_window) ? 2 : 0; pml_E_fp[0].reset(new MultiFab(amrex::convert(ba,WarpX::Ex_nodal_flag), dm, 3, nge)); pml_E_fp[1].reset(new MultiFab(amrex::convert(ba,WarpX::Ey_nodal_flag), dm, 3, nge)); @@ -458,7 +459,7 @@ PML::PML (const BoxArray& grid_ba, const DistributionMapping& grid_dm, if (do_dive_cleaning) { - pml_F_fp.reset(new MultiFab(amrex::convert(ba,IntVect::TheUnitVector()), dm, 3, 0)); + pml_F_fp.reset(new MultiFab(amrex::convert(ba,IntVect::TheUnitVector()), dm, 3, ngf)); pml_F_fp->setVal(0.0); } @@ -492,7 +493,7 @@ PML::PML (const BoxArray& grid_ba, const DistributionMapping& grid_dm, if (do_dive_cleaning) { - pml_F_cp.reset(new MultiFab(amrex::convert(cba,IntVect::TheUnitVector()), cdm, 3, 0)); + pml_F_cp.reset(new MultiFab(amrex::convert(cba,IntVect::TheUnitVector()), cdm, 3, ngf)); pml_F_cp->setVal(0.0); } |