slowEBCO Class Reference

#include <slowEBCO.H>

Inheritance diagram for slowEBCO:

Public Member Functions
Real	getSafety ()
virtual void	setAlphaAndBeta (const Real &a_alpha, const Real &a_beta)
virtual void	diagonalScale (LevelData< EBCellFAB > &a_rhs, bool a_kappaWeighted)
virtual void	divideByIdentityCoef (LevelData< EBCellFAB > &a_rhs)
virtual void	fillGrad (const LevelData< EBCellFAB > &a_phi)
	a leveltgaism More...
virtual void	getFlux (EBFluxFAB &a_flux, const LevelData< EBCellFAB > &a_data, const Box &a_grid, const DataIndex &a_dit, Real a_scale)
	a leveltgaism More...
void	getFlux (EBFaceFAB &a_fluxCentroid, const EBCellFAB &a_phi, const Box &a_ghostedBox, const Box &a_fabBox, const ProblemDomain &a_domain, const EBISBox &a_ebisBox, const Real &a_dx, const DataIndex &a_datInd, const int &a_idir)
virtual	~slowEBCO ()
	slowEBCO ()
void	AMRResidualNC (LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_rhs, bool a_homogeneousBC, AMRLevelOp< LevelData< EBCellFAB > > *a_finerOp)
void	AMROperatorNC (LevelData< EBCellFAB > &a_LofPhi, const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi, bool a_homogeneousBC, AMRLevelOp< LevelData< EBCellFAB > > *a_finerOp)
	slowEBCO (const EBLevelGrid &a_eblgFine, const EBLevelGrid &a_eblg, const EBLevelGrid &a_eblgCoar, const EBLevelGrid &a_eblgCoarMG, const RefCountedPtr< EBQuadCFInterp > &a_quadCFI, const RefCountedPtr< ConductivityBaseDomainBC > &a_domainBC, const RefCountedPtr< ConductivityBaseEBBC > &a_ebBC, const Real &a_dx, const Real &a_dxCoar, const int &a_refToFine, const int &a_refToCoar, const bool &a_hasFine, const bool &a_hasCoar, const bool &a_hasMGObjects, const bool &a_layoutChanged, const Real &a_alpha, const Real &a_beta, const RefCountedPtr< LevelData< EBCellFAB > > &a_acoef, const RefCountedPtr< LevelData< EBFluxFAB > > &a_bcoef, const RefCountedPtr< LevelData< BaseIVFAB< Real > > > &a_bcoIrreg, const IntVect &a_ghostCellsPhi, const IntVect &a_ghostCellsRHS, const int &a_relaxType)
virtual void	residual (LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_rhs, bool a_homogeneousPhysBC=false)
virtual void	preCond (LevelData< EBCellFAB > &a_opPhi, const LevelData< EBCellFAB > &a_phi)
virtual void	applyOp (LevelData< EBCellFAB > &a_opPhi, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > *const a_phiCoarse, const bool &a_homogeneousPhysBC, const bool &a_homogeneousCFBC)
virtual void	applyOpNoBoundary (LevelData< EBCellFAB > &a_opPhi, const LevelData< EBCellFAB > &a_phi)
	virtual function called by LevelTGA More...
virtual void	applyOp (LevelData< EBCellFAB > &a_opPhi, const LevelData< EBCellFAB > &a_phi, bool a_homogeneousPhysBC)
virtual void	create (LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs)
virtual void	createCoarsened (LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs, const int &a_refRat)
Real	AMRNorm (const LevelData< EBCellFAB > &a_coarResid, const LevelData< EBCellFAB > &a_fineResid, const int &a_refRat, const int &a_ord)
virtual void	assign (LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs)
virtual Real	dotProduct (const LevelData< EBCellFAB > &a_1, const LevelData< EBCellFAB > &a_2)
virtual void	incr (LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_x, Real a_scale)
virtual void	axby (LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_x, const LevelData< EBCellFAB > &a_y, Real a_a, Real a_b)
virtual void	scale (LevelData< EBCellFAB > &a_lhs, const Real &a_scale)
virtual Real	norm (const LevelData< EBCellFAB > &a_rhs, int a_ord)
virtual void	setToZero (LevelData< EBCellFAB > &a_lhs)
virtual void	setVal (LevelData< EBCellFAB > &a_lhs, const Real &a_value)
virtual void	createCoarser (LevelData< EBCellFAB > &a_coarse, const LevelData< EBCellFAB > &a_fine, bool a_ghosted)
virtual void	relax (LevelData< EBCellFAB > &a_e, const LevelData< EBCellFAB > &a_residual, int a_iterations)
virtual void	relaxGauSai (LevelData< EBCellFAB > &a_e, const LevelData< EBCellFAB > &a_residual, int a_iterations)
virtual void	relaxPoiJac (LevelData< EBCellFAB > &a_e, const LevelData< EBCellFAB > &a_residual, int a_iterations)
virtual void	restrictResidual (LevelData< EBCellFAB > &a_resCoarse, LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_rhsFine)
virtual void	prolongIncrement (LevelData< EBCellFAB > &a_phiThisLevel, const LevelData< EBCellFAB > &a_correctCoarse)
virtual int	refToCoarser ()
virtual int	refToFiner ()
virtual void	AMRResidual (LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiCoarse, const LevelData< EBCellFAB > &a_rhs, bool a_homogeneousBC, AMRLevelOp< LevelData< EBCellFAB > > *a_finerOp)
virtual void	AMRResidualNF (LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiCoarse, const LevelData< EBCellFAB > &a_rhs, bool a_homogeneousBC)
virtual void	AMROperator (LevelData< EBCellFAB > &a_LofPhi, const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiCoarse, bool a_homogeneousBC, AMRLevelOp< LevelData< EBCellFAB > > *a_finerOp)
virtual void	AMROperatorNF (LevelData< EBCellFAB > &a_LofPhi, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiCoarse, bool a_homogeneousBC)
virtual void	AMRRestrict (LevelData< EBCellFAB > &a_resCoarse, const LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_coarseCorrection, bool a_skip_res=false)
virtual void	AMRProlong (LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_coarseCorrection)
virtual void	AMRUpdateResidual (LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_coarseCorrection)
void	reflux (LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi, AMRLevelOp< LevelData< EBCellFAB > > *a_finerOp)
void	gsrbColor (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_lph, const LevelData< EBCellFAB > &a_rhs, const IntVect &a_color)
void	getDivFStencil (VoFStencil &a_vofStencil, const VolIndex &a_vof, const DataIndex &a_dit)
void	getFluxStencil (VoFStencil &a_fluxStencil, const FaceIndex &a_face, const DataIndex &a_dit)
void	getFaceCenteredFluxStencil (VoFStencil &a_fluxStencil, const FaceIndex &a_face, const DataIndex &a_dit)
void	incrOpRegularDir (EBCellFAB &a_lhs, const EBCellFAB &a_phi, const bool &a_homogeneous, const int &a_dir, const DataIndex &a_datInd)
void	applyOpIrregular (EBCellFAB &a_lhs, const EBCellFAB &a_phi, const bool &a_homogeneous, const DataIndex &a_datInd)
Public Member Functions inherited from LevelTGAHelmOp< LevelData< EBCellFAB >, EBFluxFAB >
	LevelTGAHelmOp ()
	LevelTGAHelmOp (bool a_isTimeIndependent)
virtual	~LevelTGAHelmOp ()
	Destructor. More...
Public Member Functions inherited from TGAHelmOp< LevelData< EBCellFAB > >
	TGAHelmOp ()
	TGAHelmOp (bool a_isTimeDependent)
virtual	~TGAHelmOp ()
	Destructor. More...
virtual void	diagonalScale (LevelData< EBCellFAB > &a_rhs)
virtual void	kappaScale (LevelData< EBCellFAB > &a_rhs)
	for eb only. kappa weight the rhs but do not multiply by the identity coefficient More...
virtual void	setTime (Real a_time)
virtual void	setTime (Real a_oldTime, Real a_mu, Real a_dt)
bool	isTimeDependent () const
	Returns true if the operator is time-dependent, false otherwise. More...
virtual LevelData< EBCellFAB > &	identityCoef ()
	Allows access to the identity coefficient data for the operator. More...
Public Member Functions inherited from AMRLevelOp< LevelData< EBCellFAB > >
virtual void	dumpAMR (Vector< LevelData< EBCellFAB > * > &a_data, string name)
virtual void	dumpLevel (LevelData< EBCellFAB > &a_data, string name)
	AMRLevelOp ()
	Constructor. More...
virtual void	dumpStuff (Vector< LevelData< EBCellFAB > * > data, string filename)
virtual	~AMRLevelOp ()
	Destructor. More...
virtual void	outputLevel (LevelData< EBCellFAB > &a_rhs, string &a_name)
virtual void	outputAMR (Vector< LevelData< EBCellFAB > * > &a_rhs, string &a_name)
virtual void	buildCopier (Copier &a_copier, const LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs)
virtual void	assignCopier (LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs, const Copier &a_copier)
virtual void	zeroCovered (LevelData< EBCellFAB > &a_lhs, LevelData< EBCellFAB > &a_rhs, const Copier &a_copier)
virtual Real	localMaxNorm (const LevelData< EBCellFAB > &a_phi)
virtual void	AMRProlongS (LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_coarseCorrection, LevelData< EBCellFAB > &a_temp, const Copier &a_copier)
virtual void	AMRProlongS_2 (LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_coarseCorrection, LevelData< EBCellFAB > &a_temp, const Copier &a_copier, const Copier &a_cornerCopier, const AMRLevelOp< LevelData< FArrayBox > > *a_crsOp)
virtual void	AMRRestrictS (LevelData< EBCellFAB > &a_resCoarse, const LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_coarseCorrection, LevelData< EBCellFAB > &scratch, bool a_skip_res=false)
virtual unsigned int	orderOfAccuracy (void) const
virtual void	enforceCFConsistency (LevelData< EBCellFAB > &a_coarseCorrection, const LevelData< EBCellFAB > &a_correction)
	This routine is for operators with orderOfAccuracy()>2. More...
Public Member Functions inherited from MGLevelOp< LevelData< EBCellFAB > >
	MGLevelOp ()
	Constructor. More...
virtual	~MGLevelOp ()
	Destructor. More...
virtual void	relaxNF (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > *a_phiCoarse, const LevelData< EBCellFAB > &a_rhs, int a_iterations)
	specialized no-fine relax function, useful for full-multigrid schemes, defaults to regular relax More...
virtual void	restrictResidual (LevelData< EBCellFAB > &a_resCoarse, LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > *a_phiCoarse, const LevelData< EBCellFAB > &a_rhsFine, bool homogeneous)
	full-multigrid version of restrictResidual, useful for FAS-type schemes. Defaults to standard restriction More...
virtual void	restrictR (LevelData< EBCellFAB > &a_phiCoarse, const LevelData< EBCellFAB > &a_phiFine)
	simple restriction function More...
void	addObserver (MGLevelOpObserver< LevelData< EBCellFAB > > *a_observer)
virtual void	applyOpMg (LevelData< EBCellFAB > &a_lhs, LevelData< EBCellFAB > &a_phi, LevelData< EBCellFAB > *a_phiCoarse, bool a_homogeneous)
	Apply an operator. More...
virtual void	residualNF (LevelData< EBCellFAB > &a_lhs, LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > *a_phiCoarse, const LevelData< EBCellFAB > &a_rhs, bool a_homogeneous=false)
void	removeObserver (MGLevelOpObserver< LevelData< EBCellFAB > > *a_observer)
void	addCoarserObserver (MGLevelOp< LevelData< EBCellFAB > > *a_operator, int a_coarseningFactor)
void	notifyObserversOfChange ()
	This should be called whenever the operator's data is updated. More...
virtual void	finerOperatorChanged (const MGLevelOp< LevelData< EBCellFAB > > &a_operator, int a_coarseningFactor)
int	numObservers () const
	Returns the number of objects observing this operator. More...
Public Member Functions inherited from LinearOp< LevelData< EBCellFAB > >
virtual	~LinearOp ()
virtual void	clear (LevelData< EBCellFAB > &a_lhs)
virtual void	assignLocal (LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs)
virtual void	mDotProduct (const LevelData< EBCellFAB > &a_1, const int a_sz, const LevelData< EBCellFAB > a_2[], Real a_mdots[])
virtual Real	dx () const
virtual void	write (const LevelData< EBCellFAB > *a, const char *filename)
Public Member Functions inherited from MGLevelOpObserver< LevelData< EBCellFAB > >
	MGLevelOpObserver ()
	Base level Constructor. Called by all subclasses. More...
virtual	~MGLevelOpObserver ()
	Destructor. More...
virtual void	operatorChanged (const MGLevelOp< LevelData< EBCellFAB > > &a_operator)
void	setObservee (MGLevelOp< LevelData< EBCellFAB > > *a_observee)
void	clearObservee ()

Protected Member Functions
void	defineStencils ()
void	calculateRelaxationCoefficient ()

Protected Attributes
int	m_relaxType
const IntVect	m_ghostCellsPhi
const IntVect	m_ghostCellsRHS
RefCountedPtr< EBQuadCFInterp >	m_quadCFIWithCoar
EBLevelGrid	m_eblg
EBLevelGrid	m_eblgFine
EBLevelGrid	m_eblgCoar
EBLevelGrid	m_eblgCoarMG
EBLevelGrid	m_eblgCoarsenedFine
RefCountedPtr< ConductivityBaseDomainBC >	m_domainBC
RefCountedPtr< ConductivityBaseEBBC >	m_ebBC
Real	m_dxFine
Real	m_dx
Real	m_dxCoar
RefCountedPtr< LevelData< EBCellFAB > >	m_acoef
RefCountedPtr< LevelData< EBFluxFAB > >	m_bcoef
RefCountedPtr< LevelData< BaseIVFAB< Real > > >	m_bcoIrreg
Real	m_alpha
Real	m_beta
LayoutData< BaseIVFAB< Real > >	m_alphaDiagWeight
LayoutData< BaseIVFAB< Real > >	m_betaDiagWeight
int	m_refToFine
int	m_refToCoar
bool	m_hasFine
bool	m_hasInterpAve
bool	m_hasCoar
EBMGAverage	m_ebAverage
EBMGInterp	m_ebInterp
LayoutData< BaseIVFAB< VoFStencil > >	m_opStencil
LevelData< EBCellFAB >	m_relCoef
LayoutData< VoFIterator >	m_vofIterIrreg
LayoutData< VoFIterator >	m_vofIterMulti
LayoutData< VoFIterator >	m_vofIterDomLo [CH_SPACEDIM]
LayoutData< VoFIterator >	m_vofIterDomHi [CH_SPACEDIM]
LayoutData< CFIVS >	m_loCFIVS [SpaceDim]
LayoutData< CFIVS >	m_hiCFIVS [SpaceDim]
EBFastFR	m_fastFR
bool	m_hasMGObjects
bool	m_layoutChanged
EBMGAverage	m_ebAverageMG
EBMGInterp	m_ebInterpMG
DisjointBoxLayout	m_dblCoarMG
EBISLayout	m_ebislCoarMG
ProblemDomain	m_domainCoarMG
Vector< IntVect >	m_colors

Static Protected Attributes
static bool	s_turnOffBCs

Private Member Functions
void	incrementFRCoar (EBFastFR &a_fluxReg, const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi)
void	incrementFRFine (EBFastFR &a_fluxReg, const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi, AMRLevelOp< LevelData< EBCellFAB > > *a_finerOp)
void	getFlux (FArrayBox &a_flux, const FArrayBox &a_phi, const Box &a_faceBox, const int &a_idir, const Real &a_dx, const DataIndex &a_datInd)
void	applyCFBCs (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > *const a_phiCoarse, bool a_homogeneousCFBC)
void	getOpVoFStencil (VoFStencil &a_stencil, const EBISBox &a_ebisbox, const VolIndex &a_vof)
void	getOpVoFStencil (VoFStencil &a_stencil, const int &a_dir, const Vector< VolIndex > &a_allMonotoneVoFs, const EBISBox &a_ebisbox, const VolIndex &a_vof, const bool &a_lowOrder)
void	getOpFaceStencil (VoFStencil &a_stencil, const Vector< VolIndex > &a_allMonotoneVofs, const EBISBox &a_ebisbox, const VolIndex &a_vof, int a_dir, const Side::LoHiSide &a_side, const FaceIndex &a_face, const bool &a_lowOrder)
void	levelJacobi (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_rhs)
void	applyHomogeneousCFBCs (LevelData< EBCellFAB > &a_phi)
void	applyHomogeneousCFBCs (EBCellFAB &a_phi, const DataIndex &a_datInd, int a_idir, Side::LoHiSide a_hiorlo)
	slowEBCO (const slowEBCO &a_opin)
void	operator= (const slowEBCO &a_opin)

Detailed Description

Operator to solve (alpha a + beta div (b grad) )phi = rhs. This follows the AMRLevelOp interface.

Constructor & Destructor Documentation

~slowEBCO()

virtual slowEBCO::~slowEBCO ( )

inlinevirtual

slowEBCO() [1/3]

slowEBCO::slowEBCO ( )

inline

References AMRNorm(), AMROperator(), AMROperatorNC(), AMROperatorNF(), AMRProlong(), AMRResidual(), AMRResidualNC(), AMRResidualNF(), AMRRestrict(), AMRUpdateResidual(), applyOp(), applyOpIrregular(), applyOpNoBoundary(), assign(), axby(), create(), createCoarsened(), createCoarser(), dotProduct(), getDivFStencil(), getFaceCenteredFluxStencil(), getFluxStencil(), gsrbColor(), incr(), incrOpRegularDir(), norm(), preCond(), prolongIncrement(), reflux(), refToCoarser(), refToFiner(), relax(), relaxGauSai(), relaxPoiJac(), residual(), restrictResidual(), scale(), setToZero(), and setVal().

slowEBCO() [2/3]

slowEBCO::slowEBCO	(	const EBLevelGrid &	a_eblgFine,
		const EBLevelGrid &	a_eblg,
		const EBLevelGrid &	a_eblgCoar,
		const EBLevelGrid &	a_eblgCoarMG,
		const RefCountedPtr< EBQuadCFInterp > &	a_quadCFI,
		const RefCountedPtr< ConductivityBaseDomainBC > &	a_domainBC,
		const RefCountedPtr< ConductivityBaseEBBC > &	a_ebBC,
		const Real &	a_dx,
		const Real &	a_dxCoar,
		const int &	a_refToFine,
		const int &	a_refToCoar,
		const bool &	a_hasFine,
		const bool &	a_hasCoar,
		const bool &	a_hasMGObjects,
		const bool &	a_layoutChanged,
		const Real &	a_alpha,
		const Real &	a_beta,
		const RefCountedPtr< LevelData< EBCellFAB > > &	a_acoef,
		const RefCountedPtr< LevelData< EBFluxFAB > > &	a_bcoef,
		const RefCountedPtr< LevelData< BaseIVFAB< Real > > > &	a_bcoIrreg,
		const IntVect &	a_ghostCellsPhi,
		const IntVect &	a_ghostCellsRHS,
		const int &	a_relaxType
	)

If you are approaching this operator from this interface, consider backing away and using slowEBCOFactory to generate these objects. Really. a_eblgFine, : grid at finer level \ a_eblg, : grid at this level \ a_eblgCoar, : grid at coarser level \ a_eblgCoarMG, : grid at intermediate multigrid level \ a_domainBC, : domain boundary conditions at this level \ a_ebBC: eb boundary conditions at this level \ a_dx: grid spacing at this level \ a_origin: offset to lowest corner of the domain \ a_refToFine: refinement ratio to finer level \ a_refToCoar: refinement ratio to coarser level \ a_hasFiner: true if there is a finer AMR level, false otherwise. \ a_hasCoarser: true if there is a coarser AMR level. \ a_hasCoarserMG: true if there is a coarser MultiGrid level. \ a_preCondIters: number of iterations to do for pre-conditioning \ a_relaxType: 0 means point Jacobi, 1 is Gauss-Seidel. \ a_acoef: coefficent of identity \ a_bcoef: coefficient of gradient.\ a_ghostCellsPhi: Number of ghost cells in phi, correction\ a_ghostCellsRhs: Number of ghost cells in RHS, residual, lphi\ Ghost cell arguments are there for caching reasons. Once you set them, an error is thrown if you send in data that does not match.

slowEBCO() [3/3]

slowEBCO::slowEBCO ( const slowEBCO &  a_opin )

inlineprivate

References MayDay::Error().

Member Function Documentation

getSafety()

Real slowEBCO::getSafety

(

)

setAlphaAndBeta()

virtual void slowEBCO::setAlphaAndBeta ( const Real &  a_alpha, const Real &  a_beta  )

virtual

Sets the scaling constants in the Helmholtz equation.

Parameters

a_alpha	The scaling constant that multiplies the identity term.
a_beta	The scaling constant that multiplies the derivative term.

Implements TGAHelmOp< LevelData< EBCellFAB > >.

diagonalScale()

virtual void slowEBCO::diagonalScale ( LevelData< EBCellFAB > &  a_rhs, bool  a_kappaWeighted  )

virtual

Scales the right hand side of the Helmholtz equation by the identity term in the operator. If you are solving rho(x) dphi/dt = L(phi), this means multiply by rho (or kappa * rho in embedded boundary calculations.

Parameters

a_rhs	The right hand side of the equation to be scaled.
a_kappaWeighted	If set to true, a_rhs will be scaled by the volume fraction in addition to the identity term.

Reimplemented from TGAHelmOp< LevelData< EBCellFAB > >.

divideByIdentityCoef()

virtual void slowEBCO::divideByIdentityCoef ( LevelData< EBCellFAB > &  a_rhs )

virtual

Divides the right hand side of the Helmholtz equation by the identity coefficient rho(x) in the equation rho(x) dphi/dt = L(phi).

Parameters

a_rhs

The right hand side of the equation to be scaled.

Implements TGAHelmOp< LevelData< EBCellFAB > >.

fillGrad()

virtual void slowEBCO::fillGrad ( const LevelData< EBCellFAB > &  a_phi )

inlinevirtual

a leveltgaism

Implements LevelTGAHelmOp< LevelData< EBCellFAB >, EBFluxFAB >.

References getFlux().

getFlux() [1/3]

virtual void slowEBCO::getFlux ( EBFluxFAB &  a_flux, const LevelData< EBCellFAB > &  a_data, const Box &  a_grid, const DataIndex &  a_dit, Real  a_scale  )

virtual

a leveltgaism

Implements LevelTGAHelmOp< LevelData< EBCellFAB >, EBFluxFAB >.

Referenced by fillGrad().

getFlux() [2/3]

void slowEBCO::getFlux	(	EBFaceFAB &	a_fluxCentroid,
		const EBCellFAB &	a_phi,
		const Box &	a_ghostedBox,
		const Box &	a_fabBox,
		const ProblemDomain &	a_domain,
		const EBISBox &	a_ebisBox,
		const Real &	a_dx,
		const DataIndex &	a_datInd,
		const int &	a_idir
	)

AMRResidualNC()

void slowEBCO::AMRResidualNC ( LevelData< EBCellFAB > &  a_residual, const LevelData< EBCellFAB > &  a_phiFine, const LevelData< EBCellFAB > &  a_phi, const LevelData< EBCellFAB > &  a_rhs, bool  a_homogeneousBC, AMRLevelOp< LevelData< EBCellFAB > > *  a_finerOp  )

virtual

a_residual = a_rhs - L(a_phiFine, a_phi) no coaser AMR level

Implements AMRLevelOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

AMROperatorNC()

void slowEBCO::AMROperatorNC ( LevelData< EBCellFAB > &  a_LofPhi, const LevelData< EBCellFAB > &  a_phiFine, const LevelData< EBCellFAB > &  a_phi, bool  a_homogeneousBC, AMRLevelOp< LevelData< EBCellFAB > > *  a_finerOp  )

virtual

apply AMR operator no coaser AMR level

Implements AMRLevelOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

residual()

virtual void slowEBCO::residual ( LevelData< EBCellFAB > &  a_lhs, const LevelData< EBCellFAB > &  a_phi, const LevelData< EBCellFAB > &  a_rhs, bool  a_homogeneous = false  )

virtual

Say you are solving L(phi) = rhs. Make a_lhs = L(a_phi) - a_rhs. If a_homogeneous is true, evaluate the operator using homogeneous boundary conditions.

Implements LinearOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

preCond()

virtual void slowEBCO::preCond ( LevelData< EBCellFAB > &  a_cor, const LevelData< EBCellFAB > &  a_residual  )

virtual

Given the current state of the residual the correction, apply your preconditioner to a_cor.

Implements LinearOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

applyOp() [1/2]

virtual void slowEBCO::applyOp ( LevelData< EBCellFAB > &  a_opPhi, const LevelData< EBCellFAB > &  a_phi, const LevelData< EBCellFAB > *const  a_phiCoarse, const bool &  a_homogeneousPhysBC, const bool &  a_homogeneousCFBC  )

virtual

This function assumes that coarse-fine boundary condtions have been dealt with.

Referenced by slowEBCO().

applyOpNoBoundary()

virtual void slowEBCO::applyOpNoBoundary ( LevelData< EBCellFAB > &  a_opPhi, const LevelData< EBCellFAB > &  a_phi  )

virtual

virtual function called by LevelTGA

Implements TGAHelmOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

applyOp() [2/2]

virtual void slowEBCO::applyOp ( LevelData< EBCellFAB > &  a_opPhi, const LevelData< EBCellFAB > &  a_phi, bool  a_homogeneousPhysBC  )

virtual

this is the linearop function. CFBC is set to homogeneous. phic is null

Implements LinearOp< LevelData< EBCellFAB > >.

create()

virtual void slowEBCO::create ( LevelData< EBCellFAB > &  a_lhs, const LevelData< EBCellFAB > &  a_rhs  )

virtual

Creat data holder a_lhs that mirrors a_rhs. You do not need to copy the data of a_rhs, just make a holder the same size.

Implements LinearOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

createCoarsened()

virtual void slowEBCO::createCoarsened ( LevelData< EBCellFAB > &  a_lhs, const LevelData< EBCellFAB > &  a_rhs, const int &  a_refRat  )

virtual

Implements AMRLevelOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

AMRNorm()

Real slowEBCO::AMRNorm ( const LevelData< EBCellFAB > &  a_coarResid, const LevelData< EBCellFAB > &  a_fineResid, const int &  a_refRat, const int &  a_ord  )

virtual

Reimplemented from AMRLevelOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

assign()

virtual void slowEBCO::assign ( LevelData< EBCellFAB > &  a_lhs, const LevelData< EBCellFAB > &  a_rhs  )

virtual

Set a_lhs equal to a_rhs.

Implements LinearOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

dotProduct()

virtual Real slowEBCO::dotProduct ( const LevelData< EBCellFAB > &  a_1, const LevelData< EBCellFAB > &  a_2  )

virtual

Compute and return the dot product of a_1 and a_2. In most contexts, this means return the sum over all data points of a_1*a_2.

Implements LinearOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

incr()

virtual void slowEBCO::incr ( LevelData< EBCellFAB > &  a_lhs, const LevelData< EBCellFAB > &  a_x, Real  a_scale  )

virtual

Increment by scaled amount (a_lhs += a_scale*a_x).

Implements LinearOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

axby()

virtual void slowEBCO::axby ( LevelData< EBCellFAB > &  a_lhs, const LevelData< EBCellFAB > &  a_x, const LevelData< EBCellFAB > &  a_y, Real  a_a, Real  a_b  )

virtual

Set input to a scaled sum (a_lhs = a_a*a_x + a_b*a_y).

Implements LinearOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

scale()

virtual void slowEBCO::scale ( LevelData< EBCellFAB > &  a_lhs, const Real &  a_scale  )

virtual

Multiply the input by a given scale (a_lhs *= a_scale).

Implements LinearOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

norm()

virtual Real slowEBCO::norm ( const LevelData< EBCellFAB > &  a_rhs, int  a_ord  )

virtual

Return the norm of a_rhs. a_ord == 0 max norm, a_ord == 1 sum(abs(a_rhs)), else, L(a_ord) norm.

Implements LinearOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

setToZero()

virtual void slowEBCO::setToZero ( LevelData< EBCellFAB > &  a_lhs )

virtual

Set a_lhs to zero.

Implements LinearOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

setVal()

virtual void slowEBCO::setVal ( LevelData< EBCellFAB > &  a_lhs, const Real &  a_value  )

virtual

Referenced by slowEBCO().

createCoarser()

virtual void slowEBCO::createCoarser ( LevelData< EBCellFAB > &  a_coarse, const LevelData< EBCellFAB > &  a_fine, bool  ghosted  )

virtual

Create a coarsened (by two) version of the input data. This does not include averaging the data. So if a_fine is over a Box of (0, 0, 0) (63, 63, 63), a_fine should be over a Box (0, 0, 0) (31, 31, 31).

Implements MGLevelOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

relax()

virtual void slowEBCO::relax ( LevelData< EBCellFAB > &  a_correction, const LevelData< EBCellFAB > &  a_residual, int  a_iterations  )

virtual

Use your relaxation operator to remove the high frequency wave numbers from the correction so that it may be averaged to a coarser refinement. A point relaxtion scheme, for example takes the form a_correction -= lambda*(L(a_correction) - a_residual).

Implements MGLevelOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

relaxGauSai()

virtual void slowEBCO::relaxGauSai ( LevelData< EBCellFAB > &  a_e, const LevelData< EBCellFAB > &  a_residual, int  a_iterations  )

virtual

Referenced by slowEBCO().

relaxPoiJac()

virtual void slowEBCO::relaxPoiJac ( LevelData< EBCellFAB > &  a_e, const LevelData< EBCellFAB > &  a_residual, int  a_iterations  )

virtual

Referenced by slowEBCO().

restrictResidual()

virtual void slowEBCO::restrictResidual ( LevelData< EBCellFAB > &  a_resCoarse, LevelData< EBCellFAB > &  a_phiFine, const LevelData< EBCellFAB > &  a_rhsFine  )

virtual

Calculate restricted residual: a_resCoarse[2h] = I[h->2h] (a_rhsFine[h] - L[h](a_phiFine[h]))

Implements MGLevelOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

prolongIncrement()

virtual void slowEBCO::prolongIncrement ( LevelData< EBCellFAB > &  a_phiThisLevel, const LevelData< EBCellFAB > &  a_correctCoarse  )

virtual

Correct the fine solution based on coarse correction: a_phiThisLevel += I[2h->h] (a_correctCoarse)

Implements MGLevelOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

refToCoarser()

virtual int slowEBCO::refToCoarser ( )

virtual

Refinement ratio between this level and coarser level. Returns 1 when there are no coarser AMRLevelOp objects

Implements AMRLevelOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

refToFiner()

virtual int slowEBCO::refToFiner ( )

virtual

Refinement ratio between this level and coarser level. Returns 1 when there are no coarser AMRLevelOp objects

Referenced by slowEBCO().

AMRResidual()

virtual void slowEBCO::AMRResidual ( LevelData< EBCellFAB > &  a_residual, const LevelData< EBCellFAB > &  a_phiFine, const LevelData< EBCellFAB > &  a_phi, const LevelData< EBCellFAB > &  a_phiCoarse, const LevelData< EBCellFAB > &  a_rhs, bool  a_homogeneousBC, AMRLevelOp< LevelData< EBCellFAB > > *  a_finerOp  )

virtual

a_residual = a_rhs - L(a_phi, a_phiFine, a_phiCoarse)

Implements AMRLevelOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

AMRResidualNF()

virtual void slowEBCO::AMRResidualNF ( LevelData< EBCellFAB > &  a_residual, const LevelData< EBCellFAB > &  a_phi, const LevelData< EBCellFAB > &  a_phiCoarse, const LevelData< EBCellFAB > &  a_rhs, bool  a_homogeneousBC  )

virtual

a_residual = a_rhs - L(a_phi, a_phiCoarse)

Implements AMRLevelOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

AMROperator()

virtual void slowEBCO::AMROperator ( LevelData< EBCellFAB > &  a_LofPhi, const LevelData< EBCellFAB > &  a_phiFine, const LevelData< EBCellFAB > &  a_phi, const LevelData< EBCellFAB > &  a_phiCoarse, bool  a_homogeneousBC, AMRLevelOp< LevelData< EBCellFAB > > *  a_finerOp  )

virtual

a_residual = a_rhs - L(a_phi, a_phiFine, a_phiCoarse)

Implements AMRLevelOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

AMROperatorNF()

virtual void slowEBCO::AMROperatorNF ( LevelData< EBCellFAB > &  a_LofPhi, const LevelData< EBCellFAB > &  a_phi, const LevelData< EBCellFAB > &  a_phiCoarse, bool  a_homogeneousBC  )

virtual

a_residual = a_rhs - L(a_phi, a_phiCoarse)

Implements AMRLevelOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

AMRRestrict()

virtual void slowEBCO::AMRRestrict ( LevelData< EBCellFAB > &  a_resCoarse, const LevelData< EBCellFAB > &  a_residual, const LevelData< EBCellFAB > &  a_correction, const LevelData< EBCellFAB > &  a_coarseCorrection, bool  a_skip_res = false  )

virtual

a_resCoarse = I[h-2h] (a_residual - L(a_correction, a_coarseCorrection))

Implements AMRLevelOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

AMRProlong()

virtual void slowEBCO::AMRProlong ( LevelData< EBCellFAB > &  a_correction, const LevelData< EBCellFAB > &  a_coarseCorrection  )

virtual

a_correction += I[2h->h](a_coarseCorrection)

Implements AMRLevelOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

AMRUpdateResidual()

virtual void slowEBCO::AMRUpdateResidual ( LevelData< EBCellFAB > &  a_residual, const LevelData< EBCellFAB > &  a_correction, const LevelData< EBCellFAB > &  a_coarseCorrection  )

virtual

a_residual = a_residual - L(a_correction, a_coarseCorrection)

Implements AMRLevelOp< LevelData< EBCellFAB > >.

Referenced by slowEBCO().

reflux()

void slowEBCO::reflux	(	LevelData< EBCellFAB > &	a_residual,
		const LevelData< EBCellFAB > &	a_phiFine,
		const LevelData< EBCellFAB > &	a_phi,
		AMRLevelOp< LevelData< EBCellFAB > > *	a_finerOp
	)

Referenced by slowEBCO().

gsrbColor()

void slowEBCO::gsrbColor	(	LevelData< EBCellFAB > &	a_phi,
		const LevelData< EBCellFAB > &	a_lph,
		const LevelData< EBCellFAB > &	a_rhs,
		const IntVect &	a_color
	)

Referenced by slowEBCO().

getDivFStencil()

void slowEBCO::getDivFStencil	(	VoFStencil &	a_vofStencil,
		const VolIndex &	a_vof,
		const DataIndex &	a_dit
	)

Referenced by slowEBCO().

getFluxStencil()

void slowEBCO::getFluxStencil	(	VoFStencil &	a_fluxStencil,
		const FaceIndex &	a_face,
		const DataIndex &	a_dit
	)

Referenced by slowEBCO().

getFaceCenteredFluxStencil()

void slowEBCO::getFaceCenteredFluxStencil	(	VoFStencil &	a_fluxStencil,
		const FaceIndex &	a_face,
		const DataIndex &	a_dit
	)

Referenced by slowEBCO().

incrOpRegularDir()

void slowEBCO::incrOpRegularDir	(	EBCellFAB &	a_lhs,
		const EBCellFAB &	a_phi,
		const bool &	a_homogeneous,
		const int &	a_dir,
		const DataIndex &	a_datInd
	)

Referenced by slowEBCO().

applyOpIrregular()

void slowEBCO::applyOpIrregular	(	EBCellFAB &	a_lhs,
		const EBCellFAB &	a_phi,
		const bool &	a_homogeneous,
		const DataIndex &	a_datInd
	)

Referenced by slowEBCO().

defineStencils()

void slowEBCO::defineStencils ( )

protected

calculateRelaxationCoefficient()

void slowEBCO::calculateRelaxationCoefficient ( )

protected

incrementFRCoar()

void slowEBCO::incrementFRCoar ( EBFastFR &  a_fluxReg, const LevelData< EBCellFAB > &  a_phiFine, const LevelData< EBCellFAB > &  a_phi  )

private

incrementFRFine()

void slowEBCO::incrementFRFine ( EBFastFR &  a_fluxReg, const LevelData< EBCellFAB > &  a_phiFine, const LevelData< EBCellFAB > &  a_phi, AMRLevelOp< LevelData< EBCellFAB > > *  a_finerOp  )

private

getFlux() [3/3]

void slowEBCO::getFlux ( FArrayBox &  a_flux, const FArrayBox &  a_phi, const Box &  a_faceBox, const int &  a_idir, const Real &  a_dx, const DataIndex &  a_datInd  )

private

applyCFBCs()

void slowEBCO::applyCFBCs ( LevelData< EBCellFAB > &  a_phi, const LevelData< EBCellFAB > *const  a_phiCoarse, bool  a_homogeneousCFBC  )

private

getOpVoFStencil() [1/2]

void slowEBCO::getOpVoFStencil ( VoFStencil &  a_stencil, const EBISBox &  a_ebisbox, const VolIndex &  a_vof  )

private

getOpVoFStencil() [2/2]

void slowEBCO::getOpVoFStencil ( VoFStencil &  a_stencil, const int &  a_dir, const Vector< VolIndex > &  a_allMonotoneVoFs, const EBISBox &  a_ebisbox, const VolIndex &  a_vof, const bool &  a_lowOrder  )

private

getOpFaceStencil()

void slowEBCO::getOpFaceStencil ( VoFStencil &  a_stencil, const Vector< VolIndex > &  a_allMonotoneVofs, const EBISBox &  a_ebisbox, const VolIndex &  a_vof, int  a_dir, const Side::LoHiSide &  a_side, const FaceIndex &  a_face, const bool &  a_lowOrder  )

private

levelJacobi()

void slowEBCO::levelJacobi ( LevelData< EBCellFAB > &  a_phi, const LevelData< EBCellFAB > &  a_rhs  )

private

applyHomogeneousCFBCs() [1/2]

void slowEBCO::applyHomogeneousCFBCs ( LevelData< EBCellFAB > &  a_phi )

private

applyHomogeneousCFBCs() [2/2]

void slowEBCO::applyHomogeneousCFBCs ( EBCellFAB &  a_phi, const DataIndex &  a_datInd, int  a_idir, Side::LoHiSide  a_hiorlo  )

private

operator=()

void slowEBCO::operator= ( const slowEBCO &  a_opin )

inlineprivate

References MayDay::Error().

Member Data Documentation

s_turnOffBCs

bool slowEBCO::s_turnOffBCs

staticprotected

m_relaxType

int slowEBCO::m_relaxType

protected

m_ghostCellsPhi

const IntVect slowEBCO::m_ghostCellsPhi

protected

m_ghostCellsRHS

const IntVect slowEBCO::m_ghostCellsRHS

protected

m_quadCFIWithCoar

RefCountedPtr<EBQuadCFInterp> slowEBCO::m_quadCFIWithCoar

protected

m_eblg

EBLevelGrid slowEBCO::m_eblg

protected

m_eblgFine

EBLevelGrid slowEBCO::m_eblgFine

protected

m_eblgCoar

EBLevelGrid slowEBCO::m_eblgCoar

protected

m_eblgCoarMG

EBLevelGrid slowEBCO::m_eblgCoarMG

protected

m_eblgCoarsenedFine

EBLevelGrid slowEBCO::m_eblgCoarsenedFine

protected

m_domainBC

RefCountedPtr<ConductivityBaseDomainBC> slowEBCO::m_domainBC

protected

m_ebBC

RefCountedPtr<ConductivityBaseEBBC> slowEBCO::m_ebBC

protected

m_dxFine

Real slowEBCO::m_dxFine

protected

m_dx

Real slowEBCO::m_dx

protected

m_dxCoar

Real slowEBCO::m_dxCoar

protected

m_acoef

RefCountedPtr<LevelData<EBCellFAB> > slowEBCO::m_acoef

protected

m_bcoef

RefCountedPtr<LevelData<EBFluxFAB> > slowEBCO::m_bcoef

protected

m_bcoIrreg

RefCountedPtr<LevelData<BaseIVFAB<Real> > > slowEBCO::m_bcoIrreg

protected

m_alpha

Real slowEBCO::m_alpha

protected

m_beta

Real slowEBCO::m_beta

protected

m_alphaDiagWeight

LayoutData<BaseIVFAB<Real> > slowEBCO::m_alphaDiagWeight

protected

m_betaDiagWeight

LayoutData<BaseIVFAB<Real> > slowEBCO::m_betaDiagWeight

protected

m_refToFine

int slowEBCO::m_refToFine

protected

m_refToCoar

int slowEBCO::m_refToCoar

protected

m_hasFine

bool slowEBCO::m_hasFine

protected

m_hasInterpAve

bool slowEBCO::m_hasInterpAve

protected

m_hasCoar

bool slowEBCO::m_hasCoar

protected

m_ebAverage

EBMGAverage slowEBCO::m_ebAverage

protected

m_ebInterp

EBMGInterp slowEBCO::m_ebInterp

protected

m_opStencil

LayoutData<BaseIVFAB<VoFStencil> > slowEBCO::m_opStencil

protected

m_relCoef

LevelData<EBCellFAB> slowEBCO::m_relCoef

protected

m_vofIterIrreg

LayoutData<VoFIterator > slowEBCO::m_vofIterIrreg

protected

m_vofIterMulti

LayoutData<VoFIterator > slowEBCO::m_vofIterMulti

protected

m_vofIterDomLo

LayoutData<VoFIterator > slowEBCO::m_vofIterDomLo[CH_SPACEDIM]

protected

m_vofIterDomHi

LayoutData<VoFIterator > slowEBCO::m_vofIterDomHi[CH_SPACEDIM]

protected

m_loCFIVS

LayoutData<CFIVS> slowEBCO::m_loCFIVS[SpaceDim]

protected

m_hiCFIVS

LayoutData<CFIVS> slowEBCO::m_hiCFIVS[SpaceDim]

protected

m_fastFR

EBFastFR slowEBCO::m_fastFR

protected

m_hasMGObjects

bool slowEBCO::m_hasMGObjects

protected

m_layoutChanged

bool slowEBCO::m_layoutChanged

protected

m_ebAverageMG

EBMGAverage slowEBCO::m_ebAverageMG

protected

m_ebInterpMG

EBMGInterp slowEBCO::m_ebInterpMG

protected

m_dblCoarMG

DisjointBoxLayout slowEBCO::m_dblCoarMG

protected

m_ebislCoarMG

EBISLayout slowEBCO::m_ebislCoarMG

protected

m_domainCoarMG

ProblemDomain slowEBCO::m_domainCoarMG

protected

m_colors

Vector<IntVect> slowEBCO::m_colors

protected

---

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

• slowEBCO.H