eb_ Class Reference

#include <ERF_EB.H>

Collaboration diagram for eb_:

Public Member Functions
	~eb_ ()
	eb_ (amrex::Geometry const &a_geom, amrex::FArrayBox const &terrain_fab, amrex::Gpu::DeviceVector< amrex::Real > &a_dz_stretched, bool is_anelastic)
	eb_ ()
void	define (int level, amrex::Geometry const &a_geom, amrex::EB2::Level const *a_eb_level, bool is_anelastic)
void	make_all_factories (int level, amrex::Geometry const &a_geom, amrex::BoxArray const &ba, amrex::DistributionMapping const &dm, amrex::EB2::Level const &a_eb_level)
void	make_cc_factory (int level, amrex::Geometry const &a_geom, amrex::BoxArray const &ba, amrex::DistributionMapping const &dm, amrex::EB2::Level const &a_eb_level)
int	nghost_basic () const
int	nghost_volume () const
int	nghost_full () const
const std::unique_ptr< amrex::EBFArrayBoxFactory > &	get_const_factory () const noexcept
void	set_connection_flags ()
eb_aux_ const *	get_u_const_factory () const noexcept
eb_aux_ const *	get_v_const_factory () const noexcept
eb_aux_ const *	get_w_const_factory () const noexcept

Private Member Functions
void	make_terrain (amrex::Geometry const &a_geom)
template<class F >
void	build_level (amrex::Geometry const &a_geom, amrex::EB2::GeometryShop< F > a_gshop)
	Construct EB levels from Geometry shop. More...
amrex::FabArray< amrex::EBCellFlagFab > &	getNonConstEBCellFlags (const amrex::EBFArrayBoxFactory &ebfact)
amrex::MultiFab &	getNonConstVolFrac (const amrex::EBFArrayBoxFactory &ebfact)
amrex::MultiCutFab &	getNonConstCentroid (const amrex::EBFArrayBoxFactory &ebfact)
amrex::Array< amrex::MultiCutFab *, AMREX_SPACEDIM >	getNonConstAreaFrac (const amrex::EBFArrayBoxFactory &ebfact)
amrex::Array< amrex::MultiCutFab *, AMREX_SPACEDIM >	getNonConstFaceCent (const amrex::EBFArrayBoxFactory &ebfact)

Private Attributes
int	m_has_eb
std::string	m_type
amrex::EBSupport	m_support_level
int	m_write_eb_surface
amrex::FabArray< amrex::EBCellFlagFab > *	m_cellflags = nullptr
amrex::EB2::Level const *	m_eb_level
	EB level constructed from building GeometryShop. More...
std::unique_ptr< amrex::EBFArrayBoxFactory >	m_factory = nullptr
eb_aux_	m_u_factory
eb_aux_	m_v_factory
eb_aux_	m_w_factory

Constructor & Destructor Documentation

~eb_()

eb_::~eb_

(

)

{
    // if (m_factory) { m_factory.reset(nullptr); }
}

eb_() [1/2]

eb_::eb_	(	amrex::Geometry const &	a_geom,
		amrex::FArrayBox const &	terrain_fab,
		amrex::Gpu::DeviceVector< amrex::Real > &	a_dz_stretched,
		bool	is_anelastic
	)

eb_() [2/2]

eb_::eb_

(

)

    : m_has_eb(0),
      m_support_level(EBSupport::full),
      m_write_eb_surface(0)
{ }

Member Function Documentation

build_level()

template<class F >

void eb_::build_level ( amrex::Geometry const &  a_geom, amrex::EB2::GeometryShop< F >  a_gshop  )

inlineprivate

Construct EB levels from Geometry shop.

    {
      int const req_lev(0);
      int const max_lev(2);
 
      amrex::EB2::Build(a_gshop, a_geom, req_lev, max_lev);
      const amrex::EB2::IndexSpace& ebis = amrex::EB2::IndexSpace::top();
      m_eb_level = &(ebis.getLevel(a_geom));
    }

define()

void eb_::define	(	int	level,
		amrex::Geometry const &	a_geom,
		amrex::EB2::Level const *	a_eb_level,
		bool	is_anelastic
	)

get_const_factory()

const std::unique_ptr<amrex::EBFArrayBoxFactory>& eb_::get_const_factory ( ) const

inlinenoexcept

{ return m_factory; }

Referenced by compute_gradp(), erf_slow_rhs_pre(), ERF::estTimeStep(), and make_buoyancy().

Here is the caller graph for this function:

get_u_const_factory()

eb_aux_ const* eb_::get_u_const_factory ( ) const

inlinenoexcept

{ return &m_u_factory; }

Referenced by AdvectionSrcForMom_EB(), compute_gradp(), and DiffusionSrcForMom_EB().

Here is the caller graph for this function:

get_v_const_factory()

eb_aux_ const* eb_::get_v_const_factory ( ) const

inlinenoexcept

{ return &m_v_factory; }

Referenced by AdvectionSrcForMom_EB(), compute_gradp(), and DiffusionSrcForMom_EB().

Here is the caller graph for this function:

get_w_const_factory()

eb_aux_ const* eb_::get_w_const_factory ( ) const

inlinenoexcept

{ return &m_w_factory; }

Referenced by AdvectionSrcForMom_EB(), compute_gradp(), and DiffusionSrcForMom_EB().

Here is the caller graph for this function:

getNonConstAreaFrac()

amrex::Array<amrex::MultiCutFab*, AMREX_SPACEDIM> eb_::getNonConstAreaFrac ( const amrex::EBFArrayBoxFactory &  ebfact )

inlineprivate

    {
      auto afrac_const = ebfact.getAreaFrac();
      amrex::Array<amrex::MultiCutFab*, AMREX_SPACEDIM> afrac;
      for (int dir = 0; dir < AMREX_SPACEDIM; ++dir) {
          afrac[dir] = const_cast<amrex::MultiCutFab*>(afrac_const[dir]);
      }
      return afrac;
    }

getNonConstCentroid()

amrex::MultiCutFab& eb_::getNonConstCentroid ( const amrex::EBFArrayBoxFactory &  ebfact )

inlineprivate

    {
      const amrex::MultiCutFab& vcent_const = ebfact.getCentroid();
      return const_cast<amrex::MultiCutFab&>(vcent_const);
    }

getNonConstEBCellFlags()

amrex::FabArray<amrex::EBCellFlagFab>& eb_::getNonConstEBCellFlags ( const amrex::EBFArrayBoxFactory &  ebfact )

inlineprivate

    {
      const amrex::FabArray<amrex::EBCellFlagFab>& flags_const = ebfact.getMultiEBCellFlagFab();
      return const_cast<amrex::FabArray<amrex::EBCellFlagFab>&>(flags_const);
    }

Referenced by set_connection_flags().

Here is the caller graph for this function:

getNonConstFaceCent()

amrex::Array<amrex::MultiCutFab*, AMREX_SPACEDIM> eb_::getNonConstFaceCent ( const amrex::EBFArrayBoxFactory &  ebfact )

inlineprivate

    {
      auto fcent_const = ebfact.getFaceCent();
      amrex::Array<amrex::MultiCutFab*, AMREX_SPACEDIM> fcent;
      for (int dir = 0; dir < AMREX_SPACEDIM; ++dir) {
          fcent[dir] = const_cast<amrex::MultiCutFab*>(fcent_const[dir]);
      }
      return fcent;
    }

getNonConstVolFrac()

amrex::MultiFab& eb_::getNonConstVolFrac ( const amrex::EBFArrayBoxFactory &  ebfact )

inlineprivate

    {
      const amrex::MultiFab& vfrac_const = ebfact.getVolFrac();
      return const_cast<amrex::MultiFab&>(vfrac_const);
    }

make_all_factories()

void eb_::make_all_factories	(	int	level,
		amrex::Geometry const &	a_geom,
		amrex::BoxArray const &	ba,
		amrex::DistributionMapping const &	dm,
		amrex::EB2::Level const &	a_eb_level
	)

{
    Print() << "making EB factory\n";
    m_factory = std::make_unique<EBFArrayBoxFactory>(a_eb_level, a_geom, ba, dm,
        Vector<int>{nghost_basic(), nghost_volume(), nghost_full()}, m_support_level);
 
    // Correct cell connectivity
    eb_::set_connection_flags();
 
    { int const idim(0);
        Print() << "making EB staggered u-factory\n";
        //m_u_factory.set_verbose();
        m_u_factory.define(level, idim, a_geom, ba, dm,
            Vector<int>{nghost_basic(), nghost_volume(), nghost_full()},
            m_factory.get());
    }
 
    { int const idim(1);
        Print() << "making EB staggered v-factory\n";
        //m_v_factory.set_verbose();
        m_v_factory.define(level, idim, a_geom, ba, dm,
            Vector<int>{nghost_basic(), nghost_volume(), nghost_full()},
            m_factory.get());
    }
 
    { int const idim(2);
        Print() << "making EB staggered w-factory\n";
        //m_w_factory.set_verbose();
        m_w_factory.define(level, idim, a_geom, ba, dm,
            Vector<int>{nghost_basic(), nghost_volume(), nghost_full()},
            m_factory.get());
    }
    Print() << "\nDone making EB factory at level = " << level << ".\n\n";
}

Here is the call graph for this function:

make_cc_factory()

void eb_::make_cc_factory	(	int	level,
		amrex::Geometry const &	a_geom,
		amrex::BoxArray const &	ba,
		amrex::DistributionMapping const &	dm,
		amrex::EB2::Level const &	a_eb_level
	)

{
    Print() << "making EB factory\n";
    m_factory = std::make_unique<EBFArrayBoxFactory>(a_eb_level, a_geom, ba, dm,
        Vector<int>{nghost_basic(), nghost_volume(), nghost_full()}, m_support_level);
 
    Print() << "\nDone making EB factory at level " << level << ".\n\n";
}

Here is the call graph for this function:

make_terrain()

void eb_::make_terrain ( amrex::Geometry const &  a_geom )

private

nghost_basic()

int eb_::nghost_basic ( ) const

inline

{ return 5; } // nghost_eb_basic ()

Referenced by make_all_factories(), and make_cc_factory().

Here is the caller graph for this function:

nghost_full()

int eb_::nghost_full ( ) const

inline

{ return 4; } // nghost_eb_full ()

Referenced by make_all_factories(), and make_cc_factory().

Here is the caller graph for this function:

nghost_volume()

int eb_::nghost_volume ( ) const

inline

{ return 5; } // nghost_eb_volume ()

Referenced by make_all_factories(), and make_cc_factory().

Here is the caller graph for this function:

set_connection_flags()

void eb_::set_connection_flags

(

)

{
    // Get non-const reference to EBCellFlagFab FabArray
    FabArray<EBCellFlagFab>& cellflag = getNonConstEBCellFlags(*m_factory);
 
    const MultiFab& volfrac = m_factory->getVolFrac();
 
    for (MFIter mfi(cellflag, false); mfi.isValid(); ++mfi) {
        const Box& bx = mfi.validbox();
        const Box gbx = amrex::grow(bx, cellflag.nGrow()-1); // Leave one cell layer
 
        Array4<EBCellFlag> const& flag = cellflag.array(mfi);
        Array4<Real const> const& vfrac = volfrac.const_array(mfi);
 
        ParallelFor(gbx, [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
        {
            for(int kk(-1); kk<=1; kk++) {
            for(int jj(-1); jj<=1; jj++) {
            for(int ii(-1); ii<=1; ii++)
            {
                if (vfrac(i+ii,j+jj,k+kk) == 0.0) {
                    flag(i,j,k).setDisconnected(ii,jj,kk);
                }
            }}}
        });
 
        ParallelFor(gbx, [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
        {
            if (vfrac(i,j,k)==0.0) {
                flag(i,j,k).setCovered();
            }
        });
 
    }
}

Referenced by make_all_factories().

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

m_cellflags

amrex::FabArray<amrex::EBCellFlagFab>* eb_::m_cellflags = nullptr

private

m_eb_level

amrex::EB2::Level const* eb_::m_eb_level

private

EB level constructed from building GeometryShop.

Referenced by build_level().

m_factory

std::unique_ptr<amrex::EBFArrayBoxFactory> eb_::m_factory = nullptr

private

Referenced by get_const_factory(), make_all_factories(), make_cc_factory(), and set_connection_flags().

m_has_eb

int eb_::m_has_eb

private

m_support_level

amrex::EBSupport eb_::m_support_level

private

Referenced by make_all_factories(), and make_cc_factory().

m_type

std::string eb_::m_type

private

m_u_factory

eb_aux_ eb_::m_u_factory

private

Referenced by get_u_const_factory(), and make_all_factories().

m_v_factory

eb_aux_ eb_::m_v_factory

private

Referenced by get_v_const_factory(), and make_all_factories().

m_w_factory

eb_aux_ eb_::m_w_factory

private

Referenced by get_w_const_factory(), and make_all_factories().

m_write_eb_surface

int eb_::m_write_eb_surface

private

---

The documentation for this class was generated from the following files:

• Source/EB/ERF_EB.H
• Source/EB/ERF_EB.cpp