Method fillLevelSet

Src/EB/AMReX_EB2_Level.cpp:906–943

}

void
Level::fillLevelSet (MultiFab& levelset, const Geometry& /*geom*/) const
{
    levelset.setVal(-1.0);
    levelset.ParallelCopy(m_levelset,0,0,1,IntVect(0),levelset.nGrowVect(),
                          -m_shift, m_geom.periodicity());

    std::vector<IntVect> pshifts = m_geom.periodicity().shiftIntVect();
    if (m_shift != 0) {
        for (auto& pit : pshifts) { pit += m_shift; }
    }

    Real cov_val = 1.0; // for covered cells

#ifdef AMREX_USE_OMP
#pragma omp parallel if (Gpu::notInLaunchRegion())
#endif
    if (!m_covered_grids.empty())
    {
        std::vector<std::pair<int,Box> > isects;
        for (MFIter mfi(levelset); mfi.isValid(); ++mfi)
        {
            const auto& lsfab = levelset.array(mfi);
            const Box& ccbx = amrex::enclosedCells(mfi.fabbox());
            for (const auto& iv : pshifts)
            {
                m_covered_grids.intersections(ccbx+iv, isects);
                for (const auto& is : isects) {
                    const Box& fbx = amrex::surroundingNodes(is.second-iv);
                    AMREX_HOST_DEVICE_PARALLEL_FOR_3D(fbx, i, j, k,
                    {
                        lsfab(i,j,k) = cov_val;
                    });
                }
            }
        }
    }
}

void
Level::fillCutCellMask (iMultiFab& cutcellmask, const Geometry&) const