From 964c3dce36c1ed80b7d94f400d2fcbe0ab58dd7f Mon Sep 17 00:00:00 2001
From: Edoardo Zoni <59625522+EZoni@users.noreply.github.com>
Date: Wed, 6 Oct 2021 13:35:00 -0700
Subject: FillBoundaryG: Add PML Exchange (#2381)

---
 Source/Parallelization/WarpXComm.cpp | 12 ++++++++++--
 1 file changed, 10 insertions(+), 2 deletions(-)

(limited to 'Source/Parallelization/WarpXComm.cpp')

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();
 
-- 
cgit v1.2.3