AMReX_EB_utils.H

Go to the documentation of this file.

#ifndef AMREX_EB_UTILS_H_
#define AMREX_EB_UTILS_H_
#include <AMReX_Config.H>
 
#include <AMReX.H>
#include <AMReX_MultiFab.H>
#include <AMReX_EB2.H>
#include <AMReX_EB2_GeometryShop.H>
 
namespace amrex {
    template <typename G>
    void FillImpFunc (MultiFab& mf, G const& gshop, Geometry const& geom)
    {
        AMREX_ALWAYS_ASSERT(mf.is_nodal());
 
        Box bounding_box = geom.Domain();
        bounding_box.surroundingNodes();
        bool extend_domain_face = EB2::ExtendDomainFace();
        for (int idim = 0; idim < AMREX_SPACEDIM; ++idim) {
            if (!extend_domain_face || geom.isPeriodic(idim)) {
                bounding_box.grow(mf.nGrowVect()[idim]);
            }
        }
 
#ifdef AMREX_USE_OMP
#pragma omp parallel if (Gpu::notInLaunchRegion())
#endif
        for(MFIter mfi(mf); mfi.isValid(); ++ mfi) {
            gshop.fillFab(mf[mfi], geom, RunOn::Gpu, bounding_box);
        }
    }
 
    void FillSignedDistance (MultiFab& mf, EB2::Level const& ls_lev,
                             EBFArrayBoxFactory const& eb_fac, int refratio,
                             bool fluid_has_positive_sign=true);
 
    void FillSignedDistance (MultiFab& mf, bool fluid_has_positive_sign=true);
}
 
#endif