FillBoundaryG: Add PML Exchange (#2381)

1 files changed, 10 insertions, 2 deletions

diff --git a/Source/Parallelization/WarpXComm.cpp b/Source/Parallelization/WarpXComm.cpp
index 1c871af5b..d03c9c721 100644
--- a/Source/Parallelization/WarpXComm.cpp
+++ b/Source/Parallelization/WarpXComm.cpp
@@ -821,7 +821,11 @@ void WarpX::FillBoundaryG (int lev, PatchType patch_type, IntVect ng)
 {
     if (patch_type == PatchType::fine && G_fp[lev])
     {
-        // TODO Exchange in PML cells will go here
+        if (do_pml && pml[lev]->ok())
+        {
+            pml[lev]->ExchangeG(patch_type, G_fp[lev].get(), do_pml_in_domain);
+            pml[lev]->FillBoundaryG(patch_type);
+        }
 
         const auto& period = Geom(lev).periodicity();
 
@@ -839,7 +843,11 @@ void WarpX::FillBoundaryG (int lev, PatchType patch_type, IntVect ng)
     }
     else if (patch_type == PatchType::coarse && G_cp[lev])
     {
-        // TODO Exchange in PML cells will go here
+        if (do_pml && pml[lev]->ok())
+        {
+            pml[lev]->ExchangeG(patch_type, G_cp[lev].get(), do_pml_in_domain);
+            pml[lev]->FillBoundaryG(patch_type);
+        }
 
         const auto& cperiod = Geom(lev-1).periodicity();