Fix a bug in EB2::Level::fillLevelSet

WeiqunZhang · WeiqunZhang · commit 7cf2dd7ec59e · 2024-10-31T15:35:10.000-07:00
We should fill ghost cells too. The linear solver used by WarpX relies on
that information.
diff --git a/Src/EB/AMReX_EB2_Level.cpp b/Src/EB/AMReX_EB2_Level.cpp
@@ -880,7 +880,7 @@ void
 Level::fillLevelSet (MultiFab& levelset, const Geometry& geom) const
 {
     levelset.setVal(-1.0);
-    levelset.ParallelCopy(m_levelset,0,0,1,0,0);
+    levelset.ParallelCopy(m_levelset,0,0,1,IntVect(0),levelset.nGrowVect(),geom.periodicity());
 
     const std::vector<IntVect>& pshifts = geom.periodicity().shiftIntVect();
 

Original file line number	Diff line number	Diff line change
`@@ -880,7 +880,7 @@ void`
`880`	`880`	`Level::fillLevelSet (MultiFab& levelset, const Geometry& geom) const`
`881`	`881`	`{`
`882`	`882`	`levelset.setVal(-1.0);`
`883`		`- levelset.ParallelCopy(m_levelset,0,0,1,0,0);`
	`883`	`+ levelset.ParallelCopy(m_levelset,0,0,1,IntVect(0),levelset.nGrowVect(),geom.periodicity());`
`884`	`884`
`885`	`885`	`const std::vector<IntVect>& pshifts = geom.periodicity().shiftIntVect();`
`886`	`886`