#include <AMReX_MLEBABecLap.H>

Inheritance diagram for amrex::MLEBABecLap:

Public Member Functions
	MLEBABecLap ()=default

	MLEBABecLap (const Vector< Geometry > &a_geom, const Vector< BoxArray > &a_grids, const Vector< DistributionMapping > &a_dmap, const LPInfo &a_info, const Vector< EBFArrayBoxFactory const * > &a_factory, int a_ncomp=1)

	~MLEBABecLap () override

	MLEBABecLap (const MLEBABecLap &)=delete

	MLEBABecLap (MLEBABecLap &&)=delete

MLEBABecLap &	operator= (const MLEBABecLap &)=delete

MLEBABecLap &	operator= (MLEBABecLap &&)=delete

void	define (const Vector< Geometry > &a_geom, const Vector< BoxArray > &a_grids, const Vector< DistributionMapping > &a_dmap, const LPInfo &a_info, const Vector< EBFArrayBoxFactory const * > &a_factory, int a_ncomp=1)

void	setPhiOnCentroid ()

void	setScalars (Real a, Real b)

void	setACoeffs (int amrlev, const MultiFab &alpha)

void	setACoeffs (int amrlev, Real alpha)

void	setBCoeffs (int amrlev, const Array< MultiFab const *, 3 > &beta, Location a_beta_loc)

void	setBCoeffs (int amrlev, const Array< MultiFab const *, 3 > &beta)

void	setBCoeffs (int amrlev, Real beta)

void	setBCoeffs (int amrlev, Vector< Real > const &beta)

void	setEBDirichlet (int amrlev, const MultiFab &phi, const MultiFab &beta)

void	setEBDirichlet (int amrlev, const MultiFab &phi, Real beta)

void	setEBDirichlet (int amrlev, const MultiFab &phi, Vector< Real > const &beta)

void	setEBHomogDirichlet (int amrlev, const MultiFab &beta)

void	setEBHomogDirichlet (int amrlev, Real beta)

void	setEBHomogDirichlet (int amrlev, Vector< Real > const &beta)

int	getNComp () const override
	Return number of components.

bool	needsUpdate () const override
	Does it need update if it's reused?

void	update () override
	Update for reuse.

std::unique_ptr< FabFactory< FArrayBox > >	makeFactory (int amrlev, int mglev) const final

bool	isCrossStencil () const override

void	applyBC (int amrlev, int mglev, MultiFab &in, BCMode bc_mode, StateMode s_mode, const MLMGBndry *bndry=nullptr, bool skip_fillboundary=false) const final

void	apply (int amrlev, int mglev, MultiFab &out, MultiFab &in, BCMode bc_mode, StateMode s_mode, const MLMGBndry *bndry=nullptr) const override

void	compGrad (int amrlev, const Array< MultiFab *, 3 > &grad, MultiFab &sol, Location loc) const final

void	prepareForSolve () override

bool	isSingular (int amrlev) const override
	Is it singular on given AMR level?

bool	isBottomSingular () const override
	Is the bottom of MG singular?

void	Fapply (int amrlev, int mglev, MultiFab &out, const MultiFab &in) const final

void	Fsmooth (int amrlev, int mglev, MultiFab &sol, const MultiFab &rhs, int redblack) const final

void	FFlux (int amrlev, const MFIter &mfi, const Array< FArrayBox *, 3 > &flux, const FArrayBox &sol, Location loc, int face_only=0) const final

void	normalize (int amrlev, int mglev, MultiFab &mf) const final

Real	getAScalar () const final

Real	getBScalar () const final

MultiFab const *	getACoeffs (int amrlev, int mglev) const final

Array< MultiFab const *, 3 >	getBCoeffs (int amrlev, int mglev) const final

std::unique_ptr< MLLinOp >	makeNLinOp (int) const final

void	restriction (int, int, MultiFab &crse, MultiFab &fine) const final

void	interpolation (int amrlev, int fmglev, MultiFab &fine, const MultiFab &crse) const final

void	averageDownSolutionRHS (int camrlev, MultiFab &crse_sol, MultiFab &crse_rhs, const MultiFab &fine_sol, const MultiFab &fine_rhs) final

void	getEBFluxes (const Vector< MultiFab * > &a_flux, const Vector< MultiFab * > &a_sol) const override

void	applyRobinBCTermsCoeffs ()

Public Member Functions inherited from amrex::MLCellABecLapT< MF >
	MLCellABecLapT ()=default

	~MLCellABecLapT () override=default

	MLCellABecLapT (const MLCellABecLapT< MF > &)=delete

	MLCellABecLapT (MLCellABecLapT< MF > &&)=delete

MLCellABecLapT< MF > &	operator= (const MLCellABecLapT< MF > &)=delete

MLCellABecLapT< MF > &	operator= (MLCellABecLapT< MF > &&)=delete

void	define (const Vector< Geometry > &a_geom, const Vector< BoxArray > &a_grids, const Vector< DistributionMapping > &a_dmap, const LPInfo &a_info=LPInfo(), const Vector< FabFactory< FAB > const * > &a_factory={})

void	define (const Vector< Geometry > &a_geom, const Vector< BoxArray > &a_grids, const Vector< DistributionMapping > &a_dmap, const Vector< iMultiFab const * > &a_overset_mask, const LPInfo &a_info=LPInfo(), const Vector< FabFactory< FAB > const * > &a_factory={})

iMultiFab const *	getOversetMask (int amrlev, int mglev) const

void	setDirichletNodesToZero (int amrlev, int mglev, MF &mf) const override

void	getFluxes (const Vector< Array< MF , 3 > > &a_flux, const Vector< MF > &a_sol, Location a_loc) const final

void	getFluxes (const Vector< MF * > &, const Vector< MF * > &) const final

void	applyInhomogNeumannTerm (int amrlev, MF &rhs) const final
	Extra terms introduced when we treat inhomogeneous Nuemann BC as homogeneous.

void	addInhomogNeumannFlux (int amrlev, const Array< MF *, 3 > &grad, MF const &sol, bool mult_bcoef) const final

void	applyOverset (int amrlev, MF &rhs) const override
	for overset solver only

Public Member Functions inherited from amrex::MLCellLinOpT< MF >
	MLCellLinOpT ()

	~MLCellLinOpT () override=default

	MLCellLinOpT (const MLCellLinOpT< MF > &)=delete

	MLCellLinOpT (MLCellLinOpT< MF > &&)=delete

MLCellLinOpT< MF > &	operator= (const MLCellLinOpT< MF > &)=delete

MLCellLinOpT< MF > &	operator= (MLCellLinOpT< MF > &&)=delete

void	define (const Vector< Geometry > &a_geom, const Vector< BoxArray > &a_grids, const Vector< DistributionMapping > &a_dmap, const LPInfo &a_info=LPInfo(), const Vector< FabFactory< FAB > const * > &a_factory={})

void	setLevelBC (int amrlev, const MF levelbcdata, const MF robinbc_a=nullptr, const MF robinbc_b=nullptr, const MF robinbc_f=nullptr) final
	Set boundary conditions for given level. For cell-centered solves only.

void	setGaussSeidel (bool flag) noexcept

virtual bool	isTensorOp () const

void	updateSolBC (int amrlev, const MF &crse_bcdata) const

void	updateCorBC (int amrlev, const MF &crse_bcdata) const

virtual void	applyBC (int amrlev, int mglev, MF &in, BCMode bc_mode, StateMode s_mode, const MLMGBndryT< MF > *bndry=nullptr, bool skip_fillboundary=false) const

BoxArray	makeNGrids (int grid_size) const

void	restriction (int, int, MF &crse, MF &fine) const override
	Restriction onto coarse MG level.

void	interpolation (int amrlev, int fmglev, MF &fine, const MF &crse) const override
	Add interpolated coarse MG level data to fine MG level data.

void	interpAssign (int amrlev, int fmglev, MF &fine, MF &crse) const override
	Overwrite fine MG level data with interpolated coarse data.

void	interpolationAmr (int famrlev, MF &fine, const MF &crse, IntVect const &nghost) const override
	Interpolation between AMR levels.

void	averageDownSolutionRHS (int camrlev, MF &crse_sol, MF &crse_rhs, const MF &fine_sol, const MF &fine_rhs) override
	Average-down data from fine AMR level to coarse AMR level.

void	apply (int amrlev, int mglev, MF &out, MF &in, BCMode bc_mode, StateMode s_mode, const MLMGBndryT< MF > *bndry=nullptr) const override
	Apply the linear operator, out = L(in)

void	smooth (int amrlev, int mglev, MF &sol, const MF &rhs, bool skip_fillboundary, int niter) const final
	Smooth.

void	solutionResidual (int amrlev, MF &resid, MF &x, const MF &b, const MF *crse_bcdata=nullptr) override
	Compute residual for solution.

void	prepareForFluxes (int amrlev, const MF *crse_bcdata=nullptr) override

void	correctionResidual (int amrlev, int mglev, MF &resid, MF &x, const MF &b, BCMode bc_mode, const MF *crse_bcdata=nullptr) final
	Compute residual for the residual-correction form, resid = b - L(x)

void	reflux (int crse_amrlev, MF &res, const MF &crse_sol, const MF &, MF &, MF &fine_sol, const MF &) const final
	Reflux at AMR coarse/fine boundary.

void	compFlux (int amrlev, const Array< MF *, 3 > &fluxes, MF &sol, Location loc) const override
	Compute fluxes.

void	compGrad (int amrlev, const Array< MF *, 3 > &grad, MF &sol, Location loc) const override
	Compute gradients of the solution.

void	applyMetricTerm (int amrlev, int mglev, MF &rhs) const final
	apply metric terms if there are any

void	unapplyMetricTerm (int amrlev, int mglev, MF &rhs) const final
	unapply metric terms if there are any

Vector< RT >	getSolvabilityOffset (int amrlev, int mglev, MF const &rhs) const override
	get offset for fixing solvability

void	fixSolvabilityByOffset (int amrlev, int mglev, MF &rhs, Vector< RT > const &offset) const override

RT	xdoty (int amrlev, int mglev, const MF &x, const MF &y, bool local) const final
	x dot y, used by the bottom solver

RT	dotProductPrecond (Vector< MF const * > const &x, Vector< MF const * > const &y) const final

RT	norm2Precond (Vector< MF const * > const &x) const final

virtual void	Fapply (int amrlev, int mglev, MF &out, const MF &in) const =0

virtual void	Fsmooth (int amrlev, int mglev, MF &sol, const MF &rhs, int redblack) const =0

virtual void	FFlux (int amrlev, const MFIter &mfi, const Array< FAB *, 3 > &flux, const FAB &sol, Location loc, int face_only=0) const =0

RT	normInf (int amrlev, MF const &mf, bool local) const override

void	averageDownAndSync (Vector< MF > &sol) const override

void	avgDownResAmr (int clev, MF &cres, MF const &fres) const override

void	beginPrecondBC () override

void	endPrecondBC () override

void	setInterpBndryHalfWidth (int w)

Public Member Functions inherited from amrex::MLLinOpT< MF >
	MLLinOpT ()=default

virtual	~MLLinOpT ()=default

	MLLinOpT (const MLLinOpT< MF > &)=delete

	MLLinOpT (MLLinOpT< MF > &&)=delete

MLLinOpT< MF > &	operator= (const MLLinOpT< MF > &)=delete

MLLinOpT< MF > &	operator= (MLLinOpT< MF > &&)=delete

void	define (const Vector< Geometry > &a_geom, const Vector< BoxArray > &a_grids, const Vector< DistributionMapping > &a_dmap, const LPInfo &a_info, const Vector< FabFactory< FAB > const * > &a_factory, bool eb_limit_coarsening=true)

virtual std::string	name () const

void	setDomainBC (const Array< BCType, 3 > &lobc, const Array< BCType, 3 > &hibc) noexcept
	Boundary of the whole domain.

void	setDomainBC (const Vector< Array< BCType, 3 > > &lobc, const Vector< Array< BCType, 3 > > &hibc) noexcept
	Boundary of the whole domain.

void	setDomainBCLoc (const Array< Real, 3 > &lo_bcloc, const Array< Real, 3 > &hi_bcloc) noexcept
	Set location of domain boundaries.

bool	needsCoarseDataForBC () const noexcept
	Needs coarse data for bc?

void	setCoarseFineBC (const MF *crse, int crse_ratio, LinOpBCType bc_type=LinOpBCType::Dirichlet) noexcept
	Set coarse/fine boundary conditions. For cell-centered solves only.

void	setCoarseFineBC (const MF *crse, IntVect const &crse_ratio, LinOpBCType bc_type=LinOpBCType::Dirichlet) noexcept

template<typename AMF , std::enable_if_t<!std::is_same_v< MF, AMF >, int > = 0>
void	setCoarseFineBC (const AMF *crse, int crse_ratio, LinOpBCType bc_type=LinOpBCType::Dirichlet) noexcept

template<typename AMF , std::enable_if_t<!std::is_same_v< MF, AMF >, int > = 0>
void	setCoarseFineBC (const AMF *crse, IntVect const &crse_ratio, LinOpBCType bc_type=LinOpBCType::Dirichlet) noexcept

template<typename AMF , std::enable_if_t<!std::is_same_v< MF, AMF >, int > = 0>
void	setLevelBC (int amrlev, const AMF levelbcdata, const AMF robinbc_a=nullptr, const AMF robinbc_b=nullptr, const AMF robinbc_f=nullptr)

void	setVerbose (int v) noexcept
	Set verbosity.

void	setMaxOrder (int o) noexcept
	Set order of interpolation at coarse/fine boundary.

int	getMaxOrder () const noexcept
	Get order of interpolation at coarse/fine boundary.

void	setEnforceSingularSolvable (bool o) noexcept

bool	getEnforceSingularSolvable () const noexcept

virtual BottomSolver	getDefaultBottomSolver () const

virtual int	getNGrow (int=0, int=0) const

virtual void	normalize (int amrlev, int mglev, MF &mf) const
	Divide mf by the diagonal component of the operator. Used by bicgstab.

virtual void	unimposeNeumannBC (int, MF &) const
	This is needed for our nodal projection solver.

virtual bool	scaleRHS (int, MF *) const
	scale RHS to fix solvability

virtual void	fixSolvabilityByOffset (int, int, MF &, Vector< RT > const &) const
	fix solvability by subtracting offset from RHS

virtual void	preparePrecond ()

virtual void	getEBFluxes (const Vector< MF * > &, const Vector< MF * > &) const

virtual bool	supportNSolve () const

virtual void	copyNSolveSolution (MF &, MF const &) const

virtual void	postSolve (Vector< MF * > const &) const

virtual void	avgDownResMG (int clev, MF &cres, MF const &fres) const

bool	isMFIterSafe (int amrlev, int mglev1, int mglev2) const

int	NAMRLevels () const noexcept
	Return the number of AMR levels.

int	NMGLevels (int amrlev) const noexcept
	Return the number of MG levels at given AMR level.

const Geometry &	Geom (int amr_lev, int mglev=0) const noexcept

Public Attributes
Real	m_a_scalar = std::numeric_limits<Real>::quiet_NaN()

Real	m_b_scalar = std::numeric_limits<Real>::quiet_NaN()

Vector< Vector< MultiFab > >	m_a_coeffs

Vector< Vector< Array< MultiFab, 3 > > >	m_b_coeffs

bool	m_scalars_set = false

bool	m_acoef_set = false

Public Attributes inherited from amrex::MLCellLinOpT< MF >
Vector< std::unique_ptr< MF > >	m_robin_bcval

Public Attributes inherited from amrex::MLLinOpT< MF >
Vector< Array< BCType, 3 > >	m_lobc

Vector< Array< BCType, 3 > >	m_hibc

Vector< Array< BCType, 3 > >	m_lobc_orig

Vector< Array< BCType, 3 > >	m_hibc_orig

Protected Member Functions
bool	isEBDirichlet () const noexcept

void	averageDownCoeffsSameAmrLevel (int amrlev, Vector< MultiFab > &a, Vector< Array< MultiFab, 3 > > &b, const Vector< MultiFab * > &b_eb)

void	averageDownCoeffs ()

void	averageDownCoeffsToCoarseAmrLevel (int flev)

bool	supportRobinBC () const noexcept override

Protected Member Functions inherited from amrex::MLCellABecLapT< MF >
bool	supportInhomogNeumannBC () const noexcept override

Protected Member Functions inherited from amrex::MLLinOpT< MF >
const Vector< int > &	AMRRefRatio () const noexcept
	Return AMR refinement ratios.

int	AMRRefRatio (int amr_lev) const noexcept
	Return AMR refinement ratio at given AMR level.

FabFactory< FAB > const *	Factory (int amr_lev, int mglev=0) const noexcept

GpuArray< BCType, 3 >	LoBC (int icomp=0) const noexcept

GpuArray< BCType, 3 >	HiBC (int icomp=0) const noexcept

bool	hasBC (BCType bct) const noexcept

bool	hasInhomogNeumannBC () const noexcept

bool	hasRobinBC () const noexcept

bool	isBottomActive () const noexcept

MPI_Comm	BottomCommunicator () const noexcept

MPI_Comm	Communicator () const noexcept

void	setCoarseFineBCLocation (const RealVect &cloc) noexcept

bool	doAgglomeration () const noexcept

bool	doConsolidation () const noexcept

bool	doSemicoarsening () const noexcept

bool	isCellCentered () const noexcept

virtual IntVect	getNGrowVectRestriction () const

virtual void	make (Vector< Vector< MF > > &mf, IntVect const &ng) const

virtual MF	make (int amrlev, int mglev, IntVect const &ng) const

virtual MF	makeAlias (MF const &mf) const

virtual MF	makeCoarseMG (int amrlev, int mglev, IntVect const &ng) const

virtual MF	makeCoarseAmr (int famrlev, IntVect const &ng) const

virtual void	resizeMultiGrid (int new_size)

bool	hasHiddenDimension () const noexcept

int	hiddenDirection () const noexcept

Box	compactify (Box const &b) const noexcept

template<typename T >
Array4< T >	compactify (Array4< T > const &a) const noexcept

template<typename T >
T	get_d0 (T const &dx, T const &dy, T const &) const noexcept

template<typename T >
T	get_d1 (T const &, T const &dy, T const &dz) const noexcept

Protected Attributes
int	m_ncomp = 1

bool	m_needs_update = true

Location	m_beta_loc

Location	m_phi_loc

Vector< Vector< iMultiFab > >	m_cc_mask

Vector< std::unique_ptr< MultiFab > >	m_eb_phi

Vector< Vector< std::unique_ptr< MultiFab > > >	m_eb_b_coeffs

Vector< int >	m_is_singular

int	m_is_eb_inhomog

Protected Attributes inherited from amrex::MLCellABecLapT< MF >
Vector< Vector< std::unique_ptr< iMultiFab > > >	m_overset_mask

LPInfo	m_lpinfo_arg

Protected Attributes inherited from amrex::MLCellLinOpT< MF >
bool	m_has_metric_term = false

Vector< std::unique_ptr< MLMGBndryT< MF > > >	m_bndry_sol

Vector< std::unique_ptr< BndryRegisterT< MF > > >	m_crse_sol_br

Vector< std::unique_ptr< MLMGBndryT< MF > > >	m_bndry_cor

Vector< std::unique_ptr< BndryRegisterT< MF > > >	m_crse_cor_br

Vector< std::unique_ptr< MLMGBndryT< MF > > >	m_bndry_sol_zero

Vector< Vector< std::unique_ptr< BndryCondLoc > > >	m_bcondloc

Vector< Vector< BndryRegisterT< MF > > >	m_undrrelxr

Vector< Vector< Array< MultiMask, 2 *3 > > >	m_maskvals

Vector< std::unique_ptr< iMultiFab > >	m_norm_fine_mask

Vector< YAFluxRegisterT< MF > >	m_fluxreg

bool	m_use_gauss_seidel = true

Protected Attributes inherited from amrex::MLLinOpT< MF >
int	mg_domain_min_width = 2

LPInfo	info

int	verbose = 0

int	maxorder = 3

bool	enforceSingularSolvable = true

int	m_num_amr_levels = 0

Vector< int >	m_amr_ref_ratio

Vector< int >	m_num_mg_levels

const MLLinOpT< MF > *	m_parent = nullptr

IntVect	m_ixtype

bool	m_do_agglomeration = false

bool	m_do_consolidation = false

bool	m_do_semicoarsening = false

Vector< IntVect >	mg_coarsen_ratio_vec

Vector< Vector< Geometry > >	m_geom
	first Vector is for amr level and second is mg level

Vector< Vector< BoxArray > >	m_grids

Vector< Vector< DistributionMapping > >	m_dmap

Vector< Vector< std::unique_ptr< FabFactory< FAB > > > >	m_factory

Vector< int >	m_domain_covered

MPI_Comm	m_default_comm = MPI_COMM_NULL

MPI_Comm	m_bottom_comm = MPI_COMM_NULL

std::unique_ptr< CommContainer >	m_raii_comm

Array< Real, 3 >	m_domain_bloc_lo {{ 0._rt , 0._rt , 0._rt }}

Array< Real, 3 >	m_domain_bloc_hi {{ 0._rt , 0._rt , 0._rt }}

bool	m_needs_coarse_data_for_bc = false

LinOpBCType	m_coarse_fine_bc_type = LinOpBCType::Dirichlet

IntVect	m_coarse_data_crse_ratio = IntVect(-1)

RealVect	m_coarse_bc_loc

const MF *	m_coarse_data_for_bc = nullptr

MF	m_coarse_data_for_bc_raii

bool	m_precond_mode = false

Additional Inherited Members
Public Types inherited from amrex::MLCellABecLapT< MF >
using	FAB = typename MF::fab_type

using	RT = typename MF::value_type

using	Location = typename MLLinOpT< MF >::Location

Public Types inherited from amrex::MLCellLinOpT< MF >
using	FAB = typename MF::fab_type

using	RT = typename MF::value_type

using	BCType = LinOpBCType

using	BCMode = typename MLLinOpT< MF >::BCMode

using	StateMode = typename MLLinOpT< MF >::StateMode

using	Location = typename MLLinOpT< MF >::Location

Public Types inherited from amrex::MLLinOpT< MF >
using	MFType = MF

using	FAB = typename FabDataType< MF >::fab_type

using	RT = typename FabDataType< MF >::value_type

using	BCType = LinOpBCType

using	BCMode = LinOpEnumType::BCMode

using	StateMode = LinOpEnumType::StateMode

using	Location = LinOpEnumType::Location

Protected Types inherited from amrex::MLCellLinOpT< MF >
using	RealTuple = Array< RT, 2 *3 >

using	BCTuple = Array< BoundCond, 2 *3 >

Static Protected Attributes inherited from amrex::MLLinOpT< MF >
static constexpr int	mg_coarsen_ratio = 2

static constexpr int	mg_box_min_width = 2

Constructor & Destructor Documentation

◆ MLEBABecLap() [1/4]

amrex::MLEBABecLap::MLEBABecLap ( )

default

◆ MLEBABecLap() [2/4]

amrex::MLEBABecLap::MLEBABecLap	(	const Vector< Geometry > &	a_geom,
		const Vector< BoxArray > &	a_grids,
		const Vector< DistributionMapping > &	a_dmap,
		const LPInfo &	a_info,
		const Vector< EBFArrayBoxFactory const * > &	a_factory,
		int	a_ncomp = `1`
	)

◆ ~MLEBABecLap()

amrex::MLEBABecLap::~MLEBABecLap ( )

overridedefault

◆ MLEBABecLap() [3/4]

amrex::MLEBABecLap::MLEBABecLap ( const MLEBABecLap & )

delete

◆ MLEBABecLap() [4/4]

amrex::MLEBABecLap::MLEBABecLap ( MLEBABecLap && )

delete

Member Function Documentation

◆ apply()

void amrex::MLEBABecLap::apply	(	int	amrlev,
		int	mglev,
		MultiFab &	out,
		MultiFab &	in,
		BCMode	bc_mode,
		StateMode	s_mode,
		const MLMGBndry *	bndry = `nullptr`
	)		const

override

◆ applyBC()

void amrex::MLEBABecLap::applyBC	(	int	amrlev,
		int	mglev,
		MultiFab &	in,
		BCMode	bc_mode,
		StateMode	s_mode,
		const MLMGBndry *	bndry = `nullptr`,
		bool	skip_fillboundary = `false`
	)		const

final

◆ applyRobinBCTermsCoeffs()

void amrex::MLEBABecLap::applyRobinBCTermsCoeffs ( )

◆ averageDownCoeffs()

void amrex::MLEBABecLap::averageDownCoeffs ( )

protected

◆ averageDownCoeffsSameAmrLevel()

void amrex::MLEBABecLap::averageDownCoeffsSameAmrLevel	(	int	amrlev,
		Vector< MultiFab > &	a,
		Vector< Array< MultiFab, 3 > > &	b,
		const Vector< MultiFab * > &	b_eb
	)

protected

◆ averageDownCoeffsToCoarseAmrLevel()

void amrex::MLEBABecLap::averageDownCoeffsToCoarseAmrLevel ( int flev )

protected

◆ averageDownSolutionRHS()

void amrex::MLEBABecLap::averageDownSolutionRHS	(	int	camrlev,
		MultiFab &	crse_sol,
		MultiFab &	crse_rhs,
		const MultiFab &	fine_sol,
		const MultiFab &	fine_rhs
	)

final

◆ compGrad()

void amrex::MLEBABecLap::compGrad	(	int	amrlev,
		const Array< MultiFab *, 3 > &	grad,
		MultiFab &	sol,
		Location	loc
	)		const

final

◆ define()

void amrex::MLEBABecLap::define	(	const Vector< Geometry > &	a_geom,
		const Vector< BoxArray > &	a_grids,
		const Vector< DistributionMapping > &	a_dmap,
		const LPInfo &	a_info,
		const Vector< EBFArrayBoxFactory const * > &	a_factory,
		int	a_ncomp = `1`
	)

◆ Fapply()

void amrex::MLEBABecLap::Fapply	(	int	amrlev,
		int	mglev,
		MultiFab &	out,
		const MultiFab &	in
	)		const

final

◆ FFlux()

void amrex::MLEBABecLap::FFlux	(	int	amrlev,
		const MFIter &	mfi,
		const Array< FArrayBox *, 3 > &	flux,
		const FArrayBox &	sol,
		Location	loc,
		int	face_only = `0`
	)		const

final

◆ Fsmooth()

void amrex::MLEBABecLap::Fsmooth	(	int	amrlev,
		int	mglev,
		MultiFab &	sol,
		const MultiFab &	rhs,
		int	redblack
	)		const

final

◆ getACoeffs()

MultiFab const * amrex::MLEBABecLap::getACoeffs	(	int	amrlev,
		int	mglev
	)		const

inlinefinalvirtual

Implements amrex::MLCellABecLapT< MF >.

◆ getAScalar()

Real amrex::MLEBABecLap::getAScalar ( ) const

inlinefinalvirtual

Implements amrex::MLCellABecLapT< MF >.

◆ getBCoeffs()

Array< MultiFab const *, 3 > amrex::MLEBABecLap::getBCoeffs	(	int	amrlev,
		int	mglev
	)		const

inlinefinalvirtual

Implements amrex::MLCellABecLapT< MF >.

◆ getBScalar()

Real amrex::MLEBABecLap::getBScalar ( ) const

inlinefinalvirtual

Implements amrex::MLCellABecLapT< MF >.

◆ getEBFluxes()

void amrex::MLEBABecLap::getEBFluxes	(	const Vector< MultiFab * > &	a_flux,
		const Vector< MultiFab * > &	a_sol
	)		const

override

◆ getNComp()

int amrex::MLEBABecLap::getNComp ( ) const

inlineoverridevirtual

Return number of components.

Reimplemented from amrex::MLLinOpT< MF >.

Reimplemented in amrex::MLEBTensorOp.

◆ interpolation()

void amrex::MLEBABecLap::interpolation	(	int	amrlev,
		int	fmglev,
		MultiFab &	fine,
		const MultiFab &	crse
	)		const

final

◆ isBottomSingular()

bool amrex::MLEBABecLap::isBottomSingular ( ) const

inlineoverridevirtual

Is the bottom of MG singular?

Implements amrex::MLLinOpT< MF >.

Reimplemented in amrex::MLEBTensorOp.

◆ isCrossStencil()

bool amrex::MLEBABecLap::isCrossStencil ( ) const

inlineoverridevirtual

Reimplemented from amrex::MLCellLinOpT< MF >.

Reimplemented in amrex::MLEBTensorOp.

◆ isEBDirichlet()

bool amrex::MLEBABecLap::isEBDirichlet ( ) const

inlineprotectednoexcept

◆ isSingular()

bool amrex::MLEBABecLap::isSingular ( int amrlev ) const

inlineoverridevirtual

Is it singular on given AMR level?

Implements amrex::MLLinOpT< MF >.

Reimplemented in amrex::MLEBTensorOp.

◆ makeFactory()

std::unique_ptr< FabFactory< FArrayBox > > amrex::MLEBABecLap::makeFactory	(	int	amrlev,
		int	mglev
	)		const

finalvirtual

Reimplemented from amrex::MLLinOpT< MF >.

◆ makeNLinOp()

std::unique_ptr< MLLinOp > amrex::MLEBABecLap::makeNLinOp ( int ) const

inlinefinalvirtual

Reimplemented from amrex::MLLinOpT< MF >.

◆ needsUpdate()

bool amrex::MLEBABecLap::needsUpdate ( ) const

inlineoverridevirtual

Does it need update if it's reused?

Reimplemented from amrex::MLCellABecLapT< MF >.

Reimplemented in amrex::MLEBTensorOp.

◆ normalize()

void amrex::MLEBABecLap::normalize	(	int	amrlev,
		int	mglev,
		MultiFab &	mf
	)		const

final

◆ operator=() [1/2]

MLEBABecLap & amrex::MLEBABecLap::operator= ( const MLEBABecLap & )

delete

◆ operator=() [2/2]

MLEBABecLap & amrex::MLEBABecLap::operator= ( MLEBABecLap && )

delete

◆ prepareForSolve()

void amrex::MLEBABecLap::prepareForSolve ( )

overridevirtual

Reimplemented from amrex::MLCellABecLapT< MF >.

Reimplemented in amrex::MLEBTensorOp.

◆ restriction()

void amrex::MLEBABecLap::restriction	(	int	amrlev,
		int	cmglev,
		MultiFab &	crse,
		MultiFab &	fine
	)		const

final

◆ setACoeffs() [1/2]

void amrex::MLEBABecLap::setACoeffs	(	int	amrlev,
		const MultiFab &	alpha
	)

◆ setACoeffs() [2/2]

void amrex::MLEBABecLap::setACoeffs	(	int	amrlev,
		Real	alpha
	)

◆ setBCoeffs() [1/4]

void amrex::MLEBABecLap::setBCoeffs	(	int	amrlev,
		const Array< MultiFab const *, 3 > &	beta
	)

inline

◆ setBCoeffs() [2/4]

void amrex::MLEBABecLap::setBCoeffs	(	int	amrlev,
		const Array< MultiFab const *, 3 > &	beta,
		Location	a_beta_loc
	)

◆ setBCoeffs() [3/4]

void amrex::MLEBABecLap::setBCoeffs	(	int	amrlev,
		Real	beta
	)

◆ setBCoeffs() [4/4]

void amrex::MLEBABecLap::setBCoeffs	(	int	amrlev,
		Vector< Real > const &	beta
	)

◆ setEBDirichlet() [1/3]

void amrex::MLEBABecLap::setEBDirichlet	(	int	amrlev,
		const MultiFab &	phi,
		const MultiFab &	beta
	)

◆ setEBDirichlet() [2/3]

void amrex::MLEBABecLap::setEBDirichlet	(	int	amrlev,
		const MultiFab &	phi,
		Real	beta
	)

◆ setEBDirichlet() [3/3]

void amrex::MLEBABecLap::setEBDirichlet	(	int	amrlev,
		const MultiFab &	phi,
		Vector< Real > const &	beta
	)

◆ setEBHomogDirichlet() [1/3]

void amrex::MLEBABecLap::setEBHomogDirichlet	(	int	amrlev,
		const MultiFab &	beta
	)

◆ setEBHomogDirichlet() [2/3]

void amrex::MLEBABecLap::setEBHomogDirichlet	(	int	amrlev,
		Real	beta
	)

◆ setEBHomogDirichlet() [3/3]

void amrex::MLEBABecLap::setEBHomogDirichlet	(	int	amrlev,
		Vector< Real > const &	beta
	)

◆ setPhiOnCentroid()

void amrex::MLEBABecLap::setPhiOnCentroid ( )

◆ setScalars()

void amrex::MLEBABecLap::setScalars	(	Real	a,
		Real	b
	)

◆ supportRobinBC()

bool amrex::MLEBABecLap::supportRobinBC ( ) const

inlineoverrideprotectedvirtualnoexcept

Reimplemented from amrex::MLLinOpT< MF >.

◆ update()

void amrex::MLEBABecLap::update ( )

overridevirtual

Update for reuse.

Reimplemented from amrex::MLCellABecLapT< MF >.

Reimplemented in amrex::MLEBTensorOp.

Member Data Documentation

◆ m_a_coeffs

Vector<Vector<MultiFab> > amrex::MLEBABecLap::m_a_coeffs

◆ m_a_scalar

Real amrex::MLEBABecLap::m_a_scalar = std::numeric_limits<Real>::quiet_NaN()

◆ m_acoef_set

bool amrex::MLEBABecLap::m_acoef_set = false

◆ m_b_coeffs

Vector<Vector<Array<MultiFab,3> > > amrex::MLEBABecLap::m_b_coeffs

◆ m_b_scalar

Real amrex::MLEBABecLap::m_b_scalar = std::numeric_limits<Real>::quiet_NaN()

◆ m_beta_loc

Location amrex::MLEBABecLap::m_beta_loc

protected

◆ m_cc_mask

Vector<Vector<iMultiFab> > amrex::MLEBABecLap::m_cc_mask

protected

◆ m_eb_b_coeffs

Vector<Vector<std::unique_ptr<MultiFab> > > amrex::MLEBABecLap::m_eb_b_coeffs

protected

◆ m_eb_phi

Vector<std::unique_ptr<MultiFab> > amrex::MLEBABecLap::m_eb_phi

protected

◆ m_is_eb_inhomog

int amrex::MLEBABecLap::m_is_eb_inhomog

mutableprotected

◆ m_is_singular

Vector<int> amrex::MLEBABecLap::m_is_singular

protected

◆ m_ncomp

int amrex::MLEBABecLap::m_ncomp = 1

protected

◆ m_needs_update

bool amrex::MLEBABecLap::m_needs_update = true

protected

◆ m_phi_loc

Location amrex::MLEBABecLap::m_phi_loc

protected

◆ m_scalars_set

bool amrex::MLEBABecLap::m_scalars_set = false

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

/home/runner/work/amrex/amrex/Src/LinearSolvers/MLMG/AMReX_MLEBABecLap.H
/home/runner/work/amrex/amrex/Src/LinearSolvers/MLMG/AMReX_MLEBABecLap.cpp
/home/runner/work/amrex/amrex/Src/LinearSolvers/MLMG/AMReX_MLEBABecLap_F.cpp

Public Member Functions

Public Attributes

Protected Member Functions

Protected Attributes

Additional Inherited Members

Constructor & Destructor Documentation

◆ MLEBABecLap() [1/4]

◆ MLEBABecLap() [2/4]

◆ ~MLEBABecLap()

◆ MLEBABecLap() [3/4]

◆ MLEBABecLap() [4/4]

Member Function Documentation

◆ apply()

◆ applyBC()

◆ applyRobinBCTermsCoeffs()

◆ averageDownCoeffs()

◆ averageDownCoeffsSameAmrLevel()

◆ averageDownCoeffsToCoarseAmrLevel()

◆ averageDownSolutionRHS()

◆ compGrad()

◆ define()

◆ Fapply()

◆ FFlux()

◆ Fsmooth()

◆ getACoeffs()

◆ getAScalar()

◆ getBCoeffs()

◆ getBScalar()

◆ getEBFluxes()

◆ getNComp()

◆ interpolation()

◆ isBottomSingular()

◆ isCrossStencil()

◆ isEBDirichlet()

◆ isSingular()

◆ makeFactory()

◆ makeNLinOp()

◆ needsUpdate()

◆ normalize()

◆ operator=() [1/2]

◆ operator=() [2/2]

◆ prepareForSolve()

◆ restriction()

◆ setACoeffs() [1/2]

◆ setACoeffs() [2/2]

◆ setBCoeffs() [1/4]

◆ setBCoeffs() [2/4]

◆ setBCoeffs() [3/4]

◆ setBCoeffs() [4/4]

◆ setEBDirichlet() [1/3]

◆ setEBDirichlet() [2/3]

◆ setEBDirichlet() [3/3]

◆ setEBHomogDirichlet() [1/3]

◆ setEBHomogDirichlet() [2/3]

◆ setEBHomogDirichlet() [3/3]

◆ setPhiOnCentroid()

◆ setScalars()

◆ supportRobinBC()

◆ update()

Member Data Documentation

◆ m_a_coeffs

◆ m_a_scalar

◆ m_acoef_set

◆ m_b_coeffs

◆ m_b_scalar

◆ m_beta_loc

◆ m_cc_mask

◆ m_eb_b_coeffs

◆ m_eb_phi

◆ m_is_eb_inhomog

◆ m_is_singular

◆ m_ncomp

◆ m_needs_update

◆ m_phi_loc

◆ m_scalars_set