diff options
author | 2020-12-02 17:23:14 -0800 | |
---|---|---|
committer | 2020-12-02 17:23:14 -0800 | |
commit | fcc6612bc88928b3c46ae4ac6311faa2163e1fac (patch) | |
tree | 4b20249478532dee20f0240c1bca5acc3d77cc65 /Source/BoundaryConditions/PML.cpp | |
parent | 7fe38c8d13d9a511ab8f2c1243a4c58ca393c813 (diff) | |
download | WarpX-fcc6612bc88928b3c46ae4ac6311faa2163e1fac.tar.gz WarpX-fcc6612bc88928b3c46ae4ac6311faa2163e1fac.tar.zst WarpX-fcc6612bc88928b3c46ae4ac6311faa2163e1fac.zip |
Box to build PML boxarray in domain (#1545)
Diffstat (limited to 'Source/BoundaryConditions/PML.cpp')
-rw-r--r-- | Source/BoundaryConditions/PML.cpp | 5 |
1 files changed, 4 insertions, 1 deletions
diff --git a/Source/BoundaryConditions/PML.cpp b/Source/BoundaryConditions/PML.cpp index 58ab7b9a4..7bf11a227 100644 --- a/Source/BoundaryConditions/PML.cpp +++ b/Source/BoundaryConditions/PML.cpp @@ -442,7 +442,10 @@ PML::PML (const BoxArray& grid_ba, const DistributionMapping& /*grid_dm*/, // In order to implement this, a reduced domain is created here (decreased by ncells in all direction) // and passed to `MakeBoxArray`, which surrounds it by PML boxes // (thus creating the PML boxes at the right position, where they overlap with the original domain) - Box domain0 = geom->Domain(); + // minimalBox provides the bounding box around grid_ba for level, lev. + // Note that this is okay to build pml inside domain for a single patch, or joint patches + // with same [min,max]. But it does not support multiple disjoint refinement patches. + Box domain0 = grid_ba.minimalBox(); for (int idim = 0; idim < AMREX_SPACEDIM; ++idim) { if ( ! geom->isPeriodic(idim)) { if (do_pml_Lo[idim]){ |