Chombo + EB + MF  3.2
Public Member Functions | Static Public Member Functions | Static Public Attributes | Protected Member Functions | Protected Attributes | Static Protected Attributes | Private Member Functions | List of all members
EBViscousTensorOp Class Reference

#include <EBViscousTensorOp.H>

Inheritance diagram for EBViscousTensorOp:
Inheritance graph
[legend]

Public Member Functions

 EBViscousTensorOp (const EBLevelGrid &a_eblgFine, const EBLevelGrid &a_eblg, const EBLevelGrid &a_eblgCoar, const EBLevelGrid &a_eblgCoarMG, const Real &a_alpha, const Real &a_beta, const RefCountedPtr< LevelData< EBCellFAB > > &a_acoef, const RefCountedPtr< LevelData< EBFluxFAB > > &a_eta, const RefCountedPtr< LevelData< EBFluxFAB > > &a_lambda, const RefCountedPtr< LevelData< BaseIVFAB< Real > > > &a_etaIrreg, const RefCountedPtr< LevelData< BaseIVFAB< Real > > > &a_lambdaIrreg, const Real &a_dx, const Real &a_dxCoar, const int &a_refToFine, const int &a_refToCoar, const RefCountedPtr< ViscousBaseDomainBC > &a_domainBC, const RefCountedPtr< ViscousBaseEBBC > &a_ebBC, const bool &a_hasMGObjects, const IntVect &a_ghostCellsPhi, const IntVect &a_ghostCellsRHS)
 
virtual ~EBViscousTensorOp ()
 Destructor. More...
 
virtual void setAlphaAndBeta (const Real &a_alpha, const Real &a_beta)
 for tga to reset stuff More...
 
void resetACoefficient (RefCountedPtr< LevelData< EBCellFAB > > &a_acoef)
 
virtual void outputLevel (LevelData< EBCellFAB > &a_rhs, string &a_name)
 
virtual void outputAMR (Vector< LevelData< EBCellFAB > * > &a_rhs, string &a_name)
 
virtual void kappaScale (LevelData< EBCellFAB > &a_rhs)
 it's tga's world—we just live in it. More...
 
void getShearStressDotGradU (LevelData< EBCellFAB > &a_source, const LevelData< EBCellFAB > &a_gradU, int a_level)
 compute (tau dot grad u) (for energy equation) More...
 
void getKappaDivSigmaU (LevelData< EBCellFAB > &a_divSigmaU, const LevelData< EBCellFAB > &a_velocity, const LevelData< EBCellFAB > *a_veloCoar, int a_level)
 compute volfrac(sigma dot grad u) (for energy equation) More...
 
void getCellCenteredCoefficients (LevelData< EBCellFAB > &a_etaCell, LevelData< EBCellFAB > &a_lambdaCell)
 
void getCCSigma (LevelData< EBCellFAB > &a_source, const LevelData< EBCellFAB > &a_gradU, const LevelData< EBCellFAB > &a_eta, const LevelData< EBCellFAB > &a_lambda)
 
virtual void diagonalScale (LevelData< EBCellFAB > &a_rhs, bool a_kappaWeighted)
 another tgaism More...
 
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)
 another leveltgaism More...
 
void getVelDotSigma (LevelData< EBFluxFAB > &a_velDotSigma, const LevelData< EBFluxFAB > &a_vel, const LevelData< EBFluxFAB > &a_sigma)
 
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)
 
void setTime (Real a_oldTime, Real a_mu, Real a_dt)
 
virtual void applyOpNoBoundary (LevelData< EBCellFAB > &a_opPhi, const LevelData< EBCellFAB > &a_phi)
 apply operator without any boundary or coarse-fine boundary conditions and no finer level More...
 
void finerOperatorChanged (const MGLevelOp< LevelData< EBCellFAB > > &a_operator, int a_coarseningFactor)
 
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)
 
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, bool a_homogeneousPhysBC)
 
void incrOpWithExternalFlux (EBCellFAB &a_lhs, const DataIndex &a_dit, const BaseIVFAB< Real > &a_ebflux)
 
void applyOp (LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_phi, bool a_homogeneous, const LevelData< BaseIVFAB< Real > > *a_ebflux)
 if the flux is null, it is not used More...
 
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 Real localMaxNorm (const LevelData< EBCellFAB > &a_rhs)
 
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 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 defineStencils ()
 (Re)define the stencils for the given coefficients. More...
 
- 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 setTime (Real a_time)
 
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 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 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...
 
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 ()
 

Static Public Member Functions

static Real staticMaxNorm (const LevelData< EBCellFAB > &a_rhs, const EBLevelGrid &a_eblg)
 
static void getFaceCenteredFluxStencil (VoFStencil &a_fluxStencil, const RefCountedPtr< LevelData< EBFluxFAB > > &a_eta, const RefCountedPtr< LevelData< EBFluxFAB > > &a_lambda, const Real &a_dx, const EBLevelGrid &a_eblg, const FaceIndex &a_face, const DataIndex &a_dit, int a_ivar)
 static for code reuse in bcs More...
 
static void getFluxStencil (VoFStencil &a_fluxStencil, const RefCountedPtr< LevelData< EBFluxFAB > > &a_eta, const RefCountedPtr< LevelData< EBFluxFAB > > &a_lambda, const Real &a_dx, const EBLevelGrid &a_eblg, const FaceIndex &a_face, const DataIndex &a_dit, int a_ivar)
 static for code reuse in bcs More...
 
static void getDivergenceStencil (VoFStencil &a_divStencil, const FaceIndex &a_face, const DataIndex &a_dit, const Real &a_dx, const EBLevelGrid &a_eblg)
 static for code reuse in bcs More...
 
static void getGradientStencil (VoFStencil &a_gradStencil, int a_ivar, int a_diffDir, const FaceIndex &a_face, const DataIndex &a_dit, const Real &a_dx, const EBLevelGrid &a_eblg)
 static for code reuse in bcs More...
 
static void doLazyRelax (bool a_doLazyRelax)
 
static void setForceNoEBCF (bool a_forceNoEBCF)
 

Static Public Attributes

static int s_whichLev
 
static int s_step
 

Protected Member Functions

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)
 
void getFlux (FArrayBox &a_flux, const FArrayBox &a_phi, const FArrayBox &a_gradPhi, const Box &a_faceBox, const int &a_idir, const DataIndex &a_datInd)
 
void cfinterp (const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiCoarse)
 
void reflux (const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi, LevelData< EBCellFAB > &residual, AMRLevelOp< LevelData< EBCellFAB > > *a_finerOp)
 
void homogeneousCFInterp (LevelData< EBCellFAB > &a_phi)
 
void getVoFStencil (VoFStencil &a_vofStencil, const VolIndex &a_vof, const DataIndex &a_dit, int a_ivar)
 
void getDivergenceStencil (VoFStencil &a_gradStencil, const FaceIndex &a_face, const DataIndex &a_dit)
 
void getGradientStencil (VoFStencil &a_gradStencil, int a_ivar, int a_diffDir, const FaceIndex &a_face, const DataIndex &a_dit)
 
void gsrbColor (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_lph, const LevelData< EBCellFAB > &a_rhs, const IntVect &a_color)
 
void getFaceCenteredFluxStencil (VoFStencil &a_fluxStencil, const FaceIndex &a_face, const DataIndex &a_dit, int a_ivar)
 
void cellGrad (EBCellFAB &a_gradPhi, const EBCellFAB &a_phi, const Box &a_grid, const DataIndex &a_datInd)
 
void incrementFRCoar (EBFastFR &a_fr, const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi)
 
void averageToCells (LevelData< EBCellFAB > &a_cellCoef, const LevelData< EBFluxFAB > &a_faceCoef, const LevelData< BaseIVFAB< Real > > &a_irregCoef)
 
void incrementFRFine (EBFastFR &a_fr, const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi, AMRLevelOp< LevelData< EBCellFAB > > *a_finerOp)
 
Real getSafety ()
 
virtual void calculateAlphaWeight ()
 
virtual void calculateRelaxationCoefficient ()
 
void averageCellToFace (EBFaceFAB &a_fluxData, const EBCellFAB &a_cellData, const Box &a_grid, const EBISBox &a_ebisBox, const ProblemDomain &a_domain, const DataIndex &a_dit, int isrc, int idst, int inco, bool a_interpolateToCentroid)
 
void averageCellToFace (LevelData< EBFluxFAB > &a_fluxData, const LevelData< EBCellFAB > &a_cellData, const DisjointBoxLayout &a_grids, const EBISLayout &a_ebisl, const ProblemDomain &a_domain, int isrc, int idst, int inco, bool a_interpolateToCentroid)
 
void faceCenteredAverageCellsToFaces (EBFaceFAB &a_faceData, const EBCellFAB &a_cellData, const Box &ccFluxBox, const EBISBox &a_ebisBox, const ProblemDomain &a_domain, const DataIndex &a_dit, int isrc, int idst, int inco)
 

Protected Attributes

LayoutData< CFIVSm_loCFIVS [CH_SPACEDIM]
 
LayoutData< CFIVSm_hiCFIVS [CH_SPACEDIM]
 
LayoutData< RefCountedPtr< DivergenceStencil > > m_divergenceStencil
 
LayoutData< TensorFineStencilSetm_hiTanStencilSets [SpaceDim]
 
LayoutData< TensorFineStencilSetm_loTanStencilSets [SpaceDim]
 
IntVect m_ghostPhi
 
IntVect m_ghostRHS
 
EBLevelGrid m_eblgFine
 
EBLevelGrid m_eblg
 
EBLevelGrid m_eblgCoar
 
EBLevelGrid m_eblgCoarMG
 
Real m_alpha
 
Real m_beta
 
LayoutData< BaseIVFAB< Real > > m_alphaDiagWeight
 
LayoutData< BaseIVFAB< Real > > m_betaDiagWeight
 
LevelData< EBCellFABm_zeroCoarse
 
RefCountedPtr< LevelData< EBCellFAB > > m_acoef
 
RefCountedPtr< LevelData< EBFluxFAB > > m_eta
 
RefCountedPtr< LevelData< EBFluxFAB > > m_lambda
 
RefCountedPtr< LevelData< BaseIVFAB< Real > > > m_etaIrreg
 
RefCountedPtr< LevelData< BaseIVFAB< Real > > > m_lambdaIrreg
 
EBFastFR m_fastFR
 
Real m_dx
 
Real m_dxCoar
 
bool m_hasFine
 
bool m_hasCoar
 
int m_refToFine
 
int m_refToCoar
 
bool m_hasMGObjects
 
IntVect m_ghostCellsPhi
 
IntVect m_ghostCellsRHS
 
LayoutData< RefCountedPtr< EBStencil > > m_opEBStencil [CH_SPACEDIM]
 
LevelData< EBCellFABm_relCoef
 
LevelData< EBCellFABm_grad
 
LayoutData< VoFIteratorm_vofIterIrreg
 
LayoutData< VoFIteratorm_vofIterMulti
 
LayoutData< VoFIteratorm_vofIterDomLo [CH_SPACEDIM]
 
LayoutData< VoFIteratorm_vofIterDomHi [CH_SPACEDIM]
 
RefCountedPtr< EBTensorCFInterpm_interpWithCoarser
 
EBMGAverage m_ebAverage
 
EBMGAverage m_ebAverageMG
 
EBMGInterp m_ebInterp
 
EBMGInterp m_ebInterpMG
 
RefCountedPtr< ViscousBaseDomainBCm_domainBC
 
RefCountedPtr< ViscousBaseEBBCm_ebBC
 
Vector< IntVectm_colors
 
Copier m_exchangeCopier
 
Copier m_exchangeCopierGrad
 
LayoutData< IntVectSetm_ivsIrregCCFlux
 

Static Protected Attributes

static bool s_turnOffBCs
 
static bool s_forceNoEBCF
 
static bool s_doLazyRelax
 

Private Member Functions

 EBViscousTensorOp ()
 weak construction is bad More...
 
 EBViscousTensorOp (const EBViscousTensorOp &a_opin)
 
void operator= (const EBViscousTensorOp &a_opin)
 

Detailed Description

This class implements an operator that solves the equation alpha a I + (divF) = alpha a I + beta*div(eta(grad B + grad B^T) + lambda* I div B ) = rhs where beta, lambda, and eta are incorporated into the flux F. It uses the AMRLevelOp interface.

Constructor & Destructor Documentation

◆ EBViscousTensorOp() [1/3]

EBViscousTensorOp::EBViscousTensorOp ( const EBLevelGrid a_eblgFine,
const EBLevelGrid a_eblg,
const EBLevelGrid a_eblgCoar,
const EBLevelGrid a_eblgCoarMG,
const Real a_alpha,
const Real a_beta,
const RefCountedPtr< LevelData< EBCellFAB > > &  a_acoef,
const RefCountedPtr< LevelData< EBFluxFAB > > &  a_eta,
const RefCountedPtr< LevelData< EBFluxFAB > > &  a_lambda,
const RefCountedPtr< LevelData< BaseIVFAB< Real > > > &  a_etaIrreg,
const RefCountedPtr< LevelData< BaseIVFAB< Real > > > &  a_lambdaIrreg,
const Real a_dx,
const Real a_dxCoar,
const int &  a_refToFine,
const int &  a_refToCoar,
const RefCountedPtr< ViscousBaseDomainBC > &  a_domainBC,
const RefCountedPtr< ViscousBaseEBBC > &  a_ebBC,
const bool &  a_hasMGObjects,
const IntVect a_ghostCellsPhi,
const IntVect a_ghostCellsRHS 
)

Constructs a viscous tensor operator using the given data. This constructor is for time-independent a and b coefficients.

◆ ~EBViscousTensorOp()

virtual EBViscousTensorOp::~EBViscousTensorOp ( )
virtual

Destructor.

◆ EBViscousTensorOp() [2/3]

EBViscousTensorOp::EBViscousTensorOp ( )
inlineprivate

weak construction is bad

References MayDay::Error().

◆ EBViscousTensorOp() [3/3]

EBViscousTensorOp::EBViscousTensorOp ( const EBViscousTensorOp a_opin)
inlineprivate

References MayDay::Error().

Member Function Documentation

◆ setAlphaAndBeta()

virtual void EBViscousTensorOp::setAlphaAndBeta ( const Real a_alpha,
const Real a_beta 
)
virtual

for tga to reset stuff

Implements TGAHelmOp< LevelData< EBCellFAB > >.

◆ resetACoefficient()

void EBViscousTensorOp::resetACoefficient ( RefCountedPtr< LevelData< EBCellFAB > > &  a_acoef)
inline

This sets the data storage for the a coefficient to a different object and recalculates the stuff that depends on it. Use this only if you know what you're doing.

References calculateAlphaWeight(), calculateRelaxationCoefficient(), and m_acoef.

◆ outputLevel()

virtual void EBViscousTensorOp::outputLevel ( LevelData< EBCellFAB > &  a_rhs,
string &  a_name 
)
inlinevirtual

◆ outputAMR()

virtual void EBViscousTensorOp::outputAMR ( Vector< LevelData< EBCellFAB > * > &  a_rhs,
string &  a_name 
)
inlinevirtual

◆ kappaScale()

virtual void EBViscousTensorOp::kappaScale ( LevelData< EBCellFAB > &  a_rhs)
inlinevirtual

◆ getShearStressDotGradU()

void EBViscousTensorOp::getShearStressDotGradU ( LevelData< EBCellFAB > &  a_source,
const LevelData< EBCellFAB > &  a_gradU,
int  a_level 
)

compute (tau dot grad u) (for energy equation)

Referenced by kappaScale().

◆ getKappaDivSigmaU()

void EBViscousTensorOp::getKappaDivSigmaU ( LevelData< EBCellFAB > &  a_divSigmaU,
const LevelData< EBCellFAB > &  a_velocity,
const LevelData< EBCellFAB > *  a_veloCoar,
int  a_level 
)

compute volfrac(sigma dot grad u) (for energy equation)

Referenced by kappaScale().

◆ getCellCenteredCoefficients()

void EBViscousTensorOp::getCellCenteredCoefficients ( LevelData< EBCellFAB > &  a_etaCell,
LevelData< EBCellFAB > &  a_lambdaCell 
)

Referenced by kappaScale().

◆ getCCSigma()

void EBViscousTensorOp::getCCSigma ( LevelData< EBCellFAB > &  a_source,
const LevelData< EBCellFAB > &  a_gradU,
const LevelData< EBCellFAB > &  a_eta,
const LevelData< EBCellFAB > &  a_lambda 
)

Referenced by kappaScale().

◆ diagonalScale()

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

◆ divideByIdentityCoef()

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

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_rhsThe right hand side of the equation to be scaled.

Implements TGAHelmOp< LevelData< EBCellFAB > >.

References BoxLayout::dataIterator(), fillGrad(), EBLevelGrid::getDBL(), getFlux(), getVelDotSigma(), m_acoef, m_eblg, LayoutIterator::ok(), setTime(), and SpaceDim.

◆ fillGrad()

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

◆ getFlux() [1/3]

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

◆ getVelDotSigma()

void EBViscousTensorOp::getVelDotSigma ( LevelData< EBFluxFAB > &  a_velDotSigma,
const LevelData< EBFluxFAB > &  a_vel,
const LevelData< EBFluxFAB > &  a_sigma 
)

Referenced by divideByIdentityCoef().

◆ getFlux() [2/3]

void EBViscousTensorOp::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 
)

◆ setTime()

void EBViscousTensorOp::setTime ( Real  a_oldTime,
Real  a_mu,
Real  a_dt 
)
virtual

Sets the time-dependent state of the operator. This version of setTime allows one to linearly interpolate coefficients across an integration step, since it accepts arguments that define where in the step it is to be updated. The default implementation calls setTime(a_oldTime + a_mu * a_dt).

Parameters
a_oldTimeThe time at the beginning of the current step.
a_muThe fraction of the current step that has elapsed.
a_dtThe size of the current step.

Reimplemented from TGAHelmOp< LevelData< EBCellFAB > >.

Referenced by divideByIdentityCoef().

◆ applyOpNoBoundary()

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

◆ finerOperatorChanged()

void EBViscousTensorOp::finerOperatorChanged ( const MGLevelOp< LevelData< EBCellFAB > > &  a_operator,
int  a_coarseningFactor 
)
virtual

This is called on multigrid operators when their AMR operators are altered.

Reimplemented from MGLevelOp< LevelData< EBCellFAB > >.

Referenced by applyOpNoBoundary().

◆ AMRResidualNC()

void EBViscousTensorOp::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 applyOpNoBoundary().

◆ AMROperatorNC()

void EBViscousTensorOp::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 applyOpNoBoundary().

◆ residual()

virtual void EBViscousTensorOp::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 applyOpNoBoundary().

◆ preCond()

virtual void EBViscousTensorOp::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 applyOpNoBoundary().

◆ applyOp() [1/2]

virtual void EBViscousTensorOp::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 > >.

Referenced by applyOpNoBoundary().

◆ incrOpWithExternalFlux()

void EBViscousTensorOp::incrOpWithExternalFlux ( EBCellFAB a_lhs,
const DataIndex a_dit,
const BaseIVFAB< Real > &  a_ebflux 
)

Referenced by applyOpNoBoundary().

◆ applyOp() [2/2]

void EBViscousTensorOp::applyOp ( LevelData< EBCellFAB > &  a_lhs,
const LevelData< EBCellFAB > &  a_phi,
bool  a_homogeneous,
const LevelData< BaseIVFAB< Real > > *  a_ebflux 
)

if the flux is null, it is not used

◆ create()

virtual void EBViscousTensorOp::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 applyOpNoBoundary().

◆ createCoarsened()

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

◆ AMRNorm()

Real EBViscousTensorOp::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 applyOpNoBoundary().

◆ assign()

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

Set a_lhs equal to a_rhs.

Implements LinearOp< LevelData< EBCellFAB > >.

Referenced by applyOpNoBoundary().

◆ dotProduct()

virtual Real EBViscousTensorOp::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 applyOpNoBoundary().

◆ incr()

virtual void EBViscousTensorOp::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 applyOpNoBoundary().

◆ axby()

virtual void EBViscousTensorOp::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 applyOpNoBoundary().

◆ scale()

virtual void EBViscousTensorOp::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 applyOpNoBoundary().

◆ norm()

virtual Real EBViscousTensorOp::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 applyOpNoBoundary().

◆ staticMaxNorm()

static Real EBViscousTensorOp::staticMaxNorm ( const LevelData< EBCellFAB > &  a_rhs,
const EBLevelGrid a_eblg 
)
static

Referenced by applyOpNoBoundary().

◆ localMaxNorm()

virtual Real EBViscousTensorOp::localMaxNorm ( const LevelData< EBCellFAB > &  a_rhs)
virtual

Reimplemented from AMRLevelOp< LevelData< EBCellFAB > >.

Referenced by applyOpNoBoundary().

◆ setToZero()

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

Set a_lhs to zero.

Implements LinearOp< LevelData< EBCellFAB > >.

Referenced by applyOpNoBoundary().

◆ setVal()

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

Referenced by applyOpNoBoundary().

◆ createCoarser()

virtual void EBViscousTensorOp::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 applyOpNoBoundary().

◆ relax()

virtual void EBViscousTensorOp::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 applyOpNoBoundary().

◆ restrictResidual()

virtual void EBViscousTensorOp::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 applyOpNoBoundary().

◆ prolongIncrement()

virtual void EBViscousTensorOp::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 applyOpNoBoundary().

◆ refToCoarser()

virtual int EBViscousTensorOp::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 applyOpNoBoundary().

◆ refToFiner()

virtual int EBViscousTensorOp::refToFiner ( )
virtual

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

Referenced by applyOpNoBoundary().

◆ AMRResidual()

virtual void EBViscousTensorOp::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 applyOpNoBoundary().

◆ AMRResidualNF()

virtual void EBViscousTensorOp::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 applyOpNoBoundary().

◆ AMROperator()

virtual void EBViscousTensorOp::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 applyOpNoBoundary().

◆ AMROperatorNF()

virtual void EBViscousTensorOp::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 applyOpNoBoundary().

◆ AMRRestrict()

virtual void EBViscousTensorOp::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 applyOpNoBoundary().

◆ AMRProlong()

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

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

Implements AMRLevelOp< LevelData< EBCellFAB > >.

Referenced by applyOpNoBoundary().

◆ AMRUpdateResidual()

virtual void EBViscousTensorOp::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 applyOpNoBoundary().

◆ getFaceCenteredFluxStencil() [1/2]

static void EBViscousTensorOp::getFaceCenteredFluxStencil ( VoFStencil a_fluxStencil,
const RefCountedPtr< LevelData< EBFluxFAB > > &  a_eta,
const RefCountedPtr< LevelData< EBFluxFAB > > &  a_lambda,
const Real a_dx,
const EBLevelGrid a_eblg,
const FaceIndex a_face,
const DataIndex a_dit,
int  a_ivar 
)
static

static for code reuse in bcs

Referenced by applyOpNoBoundary(), and setForceNoEBCF().

◆ getFluxStencil()

static void EBViscousTensorOp::getFluxStencil ( VoFStencil a_fluxStencil,
const RefCountedPtr< LevelData< EBFluxFAB > > &  a_eta,
const RefCountedPtr< LevelData< EBFluxFAB > > &  a_lambda,
const Real a_dx,
const EBLevelGrid a_eblg,
const FaceIndex a_face,
const DataIndex a_dit,
int  a_ivar 
)
static

static for code reuse in bcs

Referenced by applyOpNoBoundary().

◆ getDivergenceStencil() [1/2]

static void EBViscousTensorOp::getDivergenceStencil ( VoFStencil a_divStencil,
const FaceIndex a_face,
const DataIndex a_dit,
const Real a_dx,
const EBLevelGrid a_eblg 
)
static

static for code reuse in bcs

Referenced by applyOpNoBoundary(), and setForceNoEBCF().

◆ getGradientStencil() [1/2]

static void EBViscousTensorOp::getGradientStencil ( VoFStencil a_gradStencil,
int  a_ivar,
int  a_diffDir,
const FaceIndex a_face,
const DataIndex a_dit,
const Real a_dx,
const EBLevelGrid a_eblg 
)
static

static for code reuse in bcs

Referenced by applyOpNoBoundary(), and setForceNoEBCF().

◆ doLazyRelax()

static void EBViscousTensorOp::doLazyRelax ( bool  a_doLazyRelax)
inlinestatic

References s_doLazyRelax.

◆ setForceNoEBCF()

static void EBViscousTensorOp::setForceNoEBCF ( bool  a_forceNoEBCF)
inlinestatic

◆ defineStencils()

void EBViscousTensorOp::defineStencils ( )

(Re)define the stencils for the given coefficients.

Referenced by setForceNoEBCF().

◆ incrOpRegularDir()

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

Referenced by setForceNoEBCF().

◆ applyOpIrregular()

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

Referenced by setForceNoEBCF().

◆ getFlux() [3/3]

void EBViscousTensorOp::getFlux ( FArrayBox a_flux,
const FArrayBox a_phi,
const FArrayBox a_gradPhi,
const Box a_faceBox,
const int &  a_idir,
const DataIndex a_datInd 
)
protected

◆ cfinterp()

void EBViscousTensorOp::cfinterp ( const LevelData< EBCellFAB > &  a_phi,
const LevelData< EBCellFAB > &  a_phiCoarse 
)
protected

Referenced by setForceNoEBCF().

◆ reflux()

void EBViscousTensorOp::reflux ( const LevelData< EBCellFAB > &  a_phiFine,
const LevelData< EBCellFAB > &  a_phi,
LevelData< EBCellFAB > &  residual,
AMRLevelOp< LevelData< EBCellFAB > > *  a_finerOp 
)
protected

Referenced by setForceNoEBCF().

◆ homogeneousCFInterp()

void EBViscousTensorOp::homogeneousCFInterp ( LevelData< EBCellFAB > &  a_phi)
protected

Referenced by setForceNoEBCF().

◆ getVoFStencil()

void EBViscousTensorOp::getVoFStencil ( VoFStencil a_vofStencil,
const VolIndex a_vof,
const DataIndex a_dit,
int  a_ivar 
)
protected

Referenced by setForceNoEBCF().

◆ getDivergenceStencil() [2/2]

void EBViscousTensorOp::getDivergenceStencil ( VoFStencil a_gradStencil,
const FaceIndex a_face,
const DataIndex a_dit 
)
protected

◆ getGradientStencil() [2/2]

void EBViscousTensorOp::getGradientStencil ( VoFStencil a_gradStencil,
int  a_ivar,
int  a_diffDir,
const FaceIndex a_face,
const DataIndex a_dit 
)
protected

◆ gsrbColor()

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

Referenced by setForceNoEBCF().

◆ getFaceCenteredFluxStencil() [2/2]

void EBViscousTensorOp::getFaceCenteredFluxStencil ( VoFStencil a_fluxStencil,
const FaceIndex a_face,
const DataIndex a_dit,
int  a_ivar 
)
protected

◆ cellGrad()

void EBViscousTensorOp::cellGrad ( EBCellFAB a_gradPhi,
const EBCellFAB a_phi,
const Box a_grid,
const DataIndex a_datInd 
)
protected

Referenced by setForceNoEBCF().

◆ incrementFRCoar()

void EBViscousTensorOp::incrementFRCoar ( EBFastFR a_fr,
const LevelData< EBCellFAB > &  a_phiFine,
const LevelData< EBCellFAB > &  a_phi 
)
protected

Referenced by setForceNoEBCF().

◆ averageToCells()

void EBViscousTensorOp::averageToCells ( LevelData< EBCellFAB > &  a_cellCoef,
const LevelData< EBFluxFAB > &  a_faceCoef,
const LevelData< BaseIVFAB< Real > > &  a_irregCoef 
)
protected

Referenced by setForceNoEBCF().

◆ incrementFRFine()

void EBViscousTensorOp::incrementFRFine ( EBFastFR a_fr,
const LevelData< EBCellFAB > &  a_phiFine,
const LevelData< EBCellFAB > &  a_phi,
AMRLevelOp< LevelData< EBCellFAB > > *  a_finerOp 
)
protected

Referenced by setForceNoEBCF().

◆ getSafety()

Real EBViscousTensorOp::getSafety ( )
protected

◆ calculateAlphaWeight()

virtual void EBViscousTensorOp::calculateAlphaWeight ( )
protectedvirtual

Referenced by resetACoefficient().

◆ calculateRelaxationCoefficient()

virtual void EBViscousTensorOp::calculateRelaxationCoefficient ( )
protectedvirtual

Referenced by resetACoefficient().

◆ averageCellToFace() [1/2]

void EBViscousTensorOp::averageCellToFace ( EBFaceFAB a_fluxData,
const EBCellFAB a_cellData,
const Box a_grid,
const EBISBox a_ebisBox,
const ProblemDomain a_domain,
const DataIndex a_dit,
int  isrc,
int  idst,
int  inco,
bool  a_interpolateToCentroid 
)
protected

◆ averageCellToFace() [2/2]

void EBViscousTensorOp::averageCellToFace ( LevelData< EBFluxFAB > &  a_fluxData,
const LevelData< EBCellFAB > &  a_cellData,
const DisjointBoxLayout a_grids,
const EBISLayout a_ebisl,
const ProblemDomain a_domain,
int  isrc,
int  idst,
int  inco,
bool  a_interpolateToCentroid 
)
protected

◆ faceCenteredAverageCellsToFaces()

void EBViscousTensorOp::faceCenteredAverageCellsToFaces ( EBFaceFAB a_faceData,
const EBCellFAB a_cellData,
const Box ccFluxBox,
const EBISBox a_ebisBox,
const ProblemDomain a_domain,
const DataIndex a_dit,
int  isrc,
int  idst,
int  inco 
)
protected

◆ operator=()

void EBViscousTensorOp::operator= ( const EBViscousTensorOp a_opin)
inlineprivate

References MayDay::Error().

Member Data Documentation

◆ s_whichLev

int EBViscousTensorOp::s_whichLev
static

Referenced by outputAMR(), and outputLevel().

◆ s_step

int EBViscousTensorOp::s_step
static

Referenced by outputAMR(), and outputLevel().

◆ s_turnOffBCs

bool EBViscousTensorOp::s_turnOffBCs
staticprotected

Referenced by applyOpNoBoundary().

◆ s_forceNoEBCF

bool EBViscousTensorOp::s_forceNoEBCF
staticprotected

Referenced by setForceNoEBCF().

◆ s_doLazyRelax

bool EBViscousTensorOp::s_doLazyRelax
staticprotected

Referenced by doLazyRelax().

◆ m_loCFIVS

LayoutData<CFIVS> EBViscousTensorOp::m_loCFIVS[CH_SPACEDIM]
protected

◆ m_hiCFIVS

LayoutData<CFIVS> EBViscousTensorOp::m_hiCFIVS[CH_SPACEDIM]
protected

◆ m_divergenceStencil

LayoutData<RefCountedPtr<DivergenceStencil> > EBViscousTensorOp::m_divergenceStencil
protected

◆ m_hiTanStencilSets

LayoutData<TensorFineStencilSet> EBViscousTensorOp::m_hiTanStencilSets[SpaceDim]
protected

◆ m_loTanStencilSets

LayoutData<TensorFineStencilSet> EBViscousTensorOp::m_loTanStencilSets[SpaceDim]
protected

◆ m_ghostPhi

IntVect EBViscousTensorOp::m_ghostPhi
protected

◆ m_ghostRHS

IntVect EBViscousTensorOp::m_ghostRHS
protected

◆ m_eblgFine

EBLevelGrid EBViscousTensorOp::m_eblgFine
protected

◆ m_eblg

EBLevelGrid EBViscousTensorOp::m_eblg
protected

◆ m_eblgCoar

EBLevelGrid EBViscousTensorOp::m_eblgCoar
protected

◆ m_eblgCoarMG

EBLevelGrid EBViscousTensorOp::m_eblgCoarMG
protected

◆ m_alpha

Real EBViscousTensorOp::m_alpha
protected

◆ m_beta

Real EBViscousTensorOp::m_beta
protected

◆ m_alphaDiagWeight

LayoutData<BaseIVFAB<Real> > EBViscousTensorOp::m_alphaDiagWeight
protected

◆ m_betaDiagWeight

LayoutData<BaseIVFAB<Real> > EBViscousTensorOp::m_betaDiagWeight
protected

◆ m_zeroCoarse

LevelData<EBCellFAB> EBViscousTensorOp::m_zeroCoarse
protected

◆ m_acoef

RefCountedPtr<LevelData<EBCellFAB> > EBViscousTensorOp::m_acoef
protected

"Current" (time-interpolated) value of the a coefficient. For a time-independent a coefficient, this is where the coefficient lives.

Referenced by diagonalScale(), divideByIdentityCoef(), and resetACoefficient().

◆ m_eta

RefCountedPtr<LevelData<EBFluxFAB> > EBViscousTensorOp::m_eta
protected

◆ m_lambda

RefCountedPtr<LevelData<EBFluxFAB> > EBViscousTensorOp::m_lambda
protected

◆ m_etaIrreg

RefCountedPtr<LevelData<BaseIVFAB<Real> > > EBViscousTensorOp::m_etaIrreg
protected

◆ m_lambdaIrreg

RefCountedPtr<LevelData<BaseIVFAB<Real> > > EBViscousTensorOp::m_lambdaIrreg
protected

◆ m_fastFR

EBFastFR EBViscousTensorOp::m_fastFR
protected

◆ m_dx

Real EBViscousTensorOp::m_dx
protected

◆ m_dxCoar

Real EBViscousTensorOp::m_dxCoar
protected

◆ m_hasFine

bool EBViscousTensorOp::m_hasFine
protected

◆ m_hasCoar

bool EBViscousTensorOp::m_hasCoar
protected

◆ m_refToFine

int EBViscousTensorOp::m_refToFine
protected

◆ m_refToCoar

int EBViscousTensorOp::m_refToCoar
protected

◆ m_hasMGObjects

bool EBViscousTensorOp::m_hasMGObjects
protected

◆ m_ghostCellsPhi

IntVect EBViscousTensorOp::m_ghostCellsPhi
protected

◆ m_ghostCellsRHS

IntVect EBViscousTensorOp::m_ghostCellsRHS
protected

◆ m_opEBStencil

LayoutData<RefCountedPtr<EBStencil> > EBViscousTensorOp::m_opEBStencil[CH_SPACEDIM]
protected

◆ m_relCoef

LevelData<EBCellFAB> EBViscousTensorOp::m_relCoef
protected

◆ m_grad

LevelData<EBCellFAB> EBViscousTensorOp::m_grad
protected

◆ m_vofIterIrreg

LayoutData<VoFIterator > EBViscousTensorOp::m_vofIterIrreg
protected

◆ m_vofIterMulti

LayoutData<VoFIterator > EBViscousTensorOp::m_vofIterMulti
protected

◆ m_vofIterDomLo

LayoutData<VoFIterator > EBViscousTensorOp::m_vofIterDomLo[CH_SPACEDIM]
protected

◆ m_vofIterDomHi

LayoutData<VoFIterator > EBViscousTensorOp::m_vofIterDomHi[CH_SPACEDIM]
protected

◆ m_interpWithCoarser

RefCountedPtr<EBTensorCFInterp> EBViscousTensorOp::m_interpWithCoarser
protected

◆ m_ebAverage

EBMGAverage EBViscousTensorOp::m_ebAverage
protected

◆ m_ebAverageMG

EBMGAverage EBViscousTensorOp::m_ebAverageMG
protected

◆ m_ebInterp

EBMGInterp EBViscousTensorOp::m_ebInterp
protected

◆ m_ebInterpMG

EBMGInterp EBViscousTensorOp::m_ebInterpMG
protected

◆ m_domainBC

RefCountedPtr<ViscousBaseDomainBC> EBViscousTensorOp::m_domainBC
protected

◆ m_ebBC

RefCountedPtr<ViscousBaseEBBC> EBViscousTensorOp::m_ebBC
protected

◆ m_colors

Vector<IntVect> EBViscousTensorOp::m_colors
protected

◆ m_exchangeCopier

Copier EBViscousTensorOp::m_exchangeCopier
protected

◆ m_exchangeCopierGrad

Copier EBViscousTensorOp::m_exchangeCopierGrad
protected

◆ m_ivsIrregCCFlux

LayoutData<IntVectSet> EBViscousTensorOp::m_ivsIrregCCFlux
protected

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