added FillBoundary to fluid container Ue calculation

ECP-WarpX · Dec 14, 2024 · 1f06f10 · 1f06f10
1 parent daba314
commit 1f06f10
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 3 deletions.
diff --git a/Source/Fluids/QdsmcParticleContainer.cpp b/Source/Fluids/QdsmcParticleContainer.cpp
@@ -479,8 +479,5 @@ QdsmcParticleContainer::DepositK(int lev, amrex::MultiFab &Kfield)
             Kfield, 0, Kfield.nComp(), Kfield.nGrowVect(), Kfield.nGrowVect(),
             WarpX::do_single_precision_comms, period);
 
-    // Add fill boundary ?
-    //Kfield.FillBoundary(Kfield.nGrowVect(), period);
-
     //amrex::Gpu::streamSynchronize();
 }
diff --git a/Source/Fluids/WarpXFluidContainer.cpp b/Source/Fluids/WarpXFluidContainer.cpp
@@ -1415,6 +1415,10 @@ void WarpXFluidContainer::HybridInitializeUe (
         HybridPICModel const* hybrid_model,
         int lev)
 {
+    WarpX &warpx = WarpX::GetInstance();
+    const amrex::Geometry &geom = warpx.Geom(lev);
+    const amrex::Periodicity &period = geom.periodicity();
+
     using ablastr::fields::Direction;
     using warpx::fields::FieldType;
 
@@ -1493,6 +1497,10 @@ void WarpXFluidContainer::HybridInitializeUe (
 
             });
         }
+
+    m_fields.get(name_mf_NU, Direction{0}, lev)->FillBoundary(m_fields.get(name_mf_NU, Direction{0}, lev)->nGrowVect(), period);
+    m_fields.get(name_mf_NU, Direction{1}, lev)->FillBoundary(m_fields.get(name_mf_NU, Direction{1}, lev)->nGrowVect(), period);
+    m_fields.get(name_mf_NU, Direction{2}, lev)->FillBoundary(m_fields.get(name_mf_NU, Direction{2}, lev)->nGrowVect(), period);
 }