aboutsummaryrefslogtreecommitdiff
path: root/Source/BoundaryConditions/PML.cpp
diff options
context:
space:
mode:
authorGravatar Dave Grote <grote1@llnl.gov> 2019-08-16 12:53:47 -0700
committerGravatar Dave Grote <grote1@llnl.gov> 2019-08-16 12:53:47 -0700
commita26f56418b18c98e3784d1d48073cae81b938e79 (patch)
treed609a8692e8dd29530d62c87809a54a6699b090d /Source/BoundaryConditions/PML.cpp
parentc3ce219b9d25e8d28e5a6cc5b878b3c5793cf90a (diff)
parent2d4cf600befb1ef4d5acbefd9d681ff66c11334e (diff)
downloadWarpX-a26f56418b18c98e3784d1d48073cae81b938e79.tar.gz
WarpX-a26f56418b18c98e3784d1d48073cae81b938e79.tar.zst
WarpX-a26f56418b18c98e3784d1d48073cae81b938e79.zip
Merge branch 'dev' into RZgeometry
Diffstat (limited to 'Source/BoundaryConditions/PML.cpp')
-rw-r--r--Source/BoundaryConditions/PML.cpp17
1 files changed, 12 insertions, 5 deletions
diff --git a/Source/BoundaryConditions/PML.cpp b/Source/BoundaryConditions/PML.cpp
index 96bc08af9..21d348482 100644
--- a/Source/BoundaryConditions/PML.cpp
+++ b/Source/BoundaryConditions/PML.cpp
@@ -319,11 +319,12 @@ PML::PML (const BoxArray& grid_ba, const DistributionMapping& grid_dm,
#ifdef WARPX_USE_PSATD
Real dt, int nox_fft, int noy_fft, int noz_fft, bool do_nodal,
#endif
- int do_dive_cleaning, int do_moving_window)
+ int do_dive_cleaning, int do_moving_window,
+ const amrex::IntVect do_pml_Lo, const amrex::IntVect do_pml_Hi)
: m_geom(geom),
m_cgeom(cgeom)
{
- const BoxArray& ba = MakeBoxArray(*geom, grid_ba, ncell);
+ const BoxArray& ba = MakeBoxArray(*geom, grid_ba, ncell, do_pml_Lo, do_pml_Hi);
if (ba.size() == 0) {
m_ok = false;
return;
@@ -399,7 +400,7 @@ PML::PML (const BoxArray& grid_ba, const DistributionMapping& grid_dm,
BoxArray grid_cba = grid_ba;
grid_cba.coarsen(ref_ratio);
- const BoxArray& cba = MakeBoxArray(*cgeom, grid_cba, ncell);
+ const BoxArray& cba = MakeBoxArray(*cgeom, grid_cba, ncell, do_pml_Lo, do_pml_Hi);
DistributionMapping cdm{cba};
@@ -438,12 +439,18 @@ PML::PML (const BoxArray& grid_ba, const DistributionMapping& grid_dm,
}
BoxArray
-PML::MakeBoxArray (const amrex::Geometry& geom, const amrex::BoxArray& grid_ba, int ncell)
+PML::MakeBoxArray (const amrex::Geometry& geom, const amrex::BoxArray& grid_ba, int ncell,
+ const amrex::IntVect do_pml_Lo, const amrex::IntVect do_pml_Hi)
{
Box domain = geom.Domain();
for (int idim = 0; idim < AMREX_SPACEDIM; ++idim) {
if ( ! geom.isPeriodic(idim) ) {
- domain.grow(idim, ncell);
+ if (do_pml_Lo[idim]){
+ domain.growLo(idim, ncell);
+ }
+ if (do_pml_Hi[idim]){
+ domain.growHi(idim, ncell);
+ }
}
}