diff options
| -rw-r--r-- | Source/EmbeddedBoundary/WarpXInitEB.cpp | 17 |
1 files changed, 5 insertions, 12 deletions
diff --git a/Source/EmbeddedBoundary/WarpXInitEB.cpp b/Source/EmbeddedBoundary/WarpXInitEB.cpp index a51ecd148..aa9d41f53 100644 --- a/Source/EmbeddedBoundary/WarpXInitEB.cpp +++ b/Source/EmbeddedBoundary/WarpXInitEB.cpp @@ -398,17 +398,10 @@ WarpX::ComputeDistanceToEB () { #ifdef AMREX_USE_EB BL_PROFILE("ComputeDistanceToEB"); - amrex::ParmParse pp_warpx("warpx"); - std::string impf; - pp_warpx.query("eb_implicit_function", impf); - if (! impf.empty()) { - auto eb_if_parser = makeParser(impf, {"x", "y", "z"}); - ParserIF pif(eb_if_parser.compile<3>()); - auto gshop = amrex::EB2::makeShop(pif); - amrex::FillImpFunc(*m_distance_to_eb[maxLevel()], gshop, Geom(maxLevel())); - m_distance_to_eb[maxLevel()]->negate(m_distance_to_eb[maxLevel()]->nGrow()); // signed distance f = - imp. f. - } else { - m_distance_to_eb[maxLevel()]->setVal(100.0); // some positive value - } + const amrex::EB2::IndexSpace& eb_is = amrex::EB2::IndexSpace::top(); + const amrex::EB2::Level& eb_level = eb_is.getLevel(Geom(maxLevel())); + auto const eb_fact = fieldEBFactory(maxLevel()); + + amrex::FillSignedDistance(*m_distance_to_eb[maxLevel()], eb_level, eb_fact, 1); #endif } |