ERF_SolverUtils.H File Reference

#include "AMReX_LO_BCTYPES.H"

Include dependency graph for ERF_SolverUtils.H:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions
Array< LinOpBCType, AMREX_SPACEDIM >	get_lo_projection_bc (Geometry const &lev_geom, Array< std::string, 2 *AMREX_SPACEDIM > l_domain_bc_type)
Array< LinOpBCType, AMREX_SPACEDIM >	get_hi_projection_bc (Geometry const &lev_geom, Array< std::string, 2 *AMREX_SPACEDIM > l_domain_bc_type)

Function Documentation

get_hi_projection_bc()

Array<LinOpBCType,AMREX_SPACEDIM> get_hi_projection_bc ( Geometry const &  lev_geom, Array< std::string, 2 *AMREX_SPACEDIM >  l_domain_bc_type  )

inline

Define high-side domain boundary conditions for the projection solve.

Parameters

lev_geom	Geometry used to identify periodic high-side boundaries
l_domain_bc_type	Domain boundary-condition names

Returns

High-side linear-operator boundary condition in each direction

{
    amrex::Array<amrex::LinOpBCType,AMREX_SPACEDIM> r;
    for (int dir = 0; dir < AMREX_SPACEDIM; ++dir) {
        if ( lev_geom.isPeriodic(dir) )
        {
            r[dir] = LinOpBCType::Periodic;
        } else {
            auto bc_type = l_domain_bc_type[Orientation(dir,Orientation::high)];
            if (bc_type == "Outflow" || bc_type == "Open") {
                r[dir] = LinOpBCType::Dirichlet;
            } else {
                r[dir] = LinOpBCType::Neumann;
            }
        }
    }
    return r;
}

Referenced by ERF::project_velocity_tb(), solve_with_EB_mlmg(), and solve_with_mlmg().

Here is the caller graph for this function:

get_lo_projection_bc()

Array<LinOpBCType,AMREX_SPACEDIM> get_lo_projection_bc ( Geometry const &  lev_geom, Array< std::string, 2 *AMREX_SPACEDIM >  l_domain_bc_type  )

inline

Define the domain boundary conditions for the (optional) Poisson solve if we want to enforce that the initial conditions satisfy the constraint

Parameters

lev_geom	Geometry used to identify periodic low-side boundaries
l_domain_bc_type	Domain boundary-condition names

Returns

Low-side linear-operator boundary condition in each direction

{
    amrex::Array<amrex::LinOpBCType,AMREX_SPACEDIM> r;
    for (int dir = 0; dir < AMREX_SPACEDIM; ++dir) {
        if ( lev_geom.isPeriodic(dir) )
        {
            r[dir] = LinOpBCType::Periodic;
        } else {
            auto bc_type = l_domain_bc_type[Orientation(dir,Orientation::low)];
            if (bc_type == "Outflow" || bc_type == "Open") {
                r[dir] = LinOpBCType::Dirichlet;
            } else {
                r[dir] = LinOpBCType::Neumann;
            }
        }
    }
    return r;
}

Referenced by ERF::project_velocity_tb(), solve_with_EB_mlmg(), and solve_with_mlmg().

Here is the caller graph for this function: