Chombo + EB
3.0
|
#include <EBViscousTensorOp.H>
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) | |
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_acoef0, const RefCountedPtr< LevelData< EBCellFAB > > &a_acoef1, 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... | |
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) |
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 | 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) |
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 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 | AMRRestrictS (LevelData< EBCellFAB > &a_resCoarse, const LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_coarseCorrection, LevelData< EBCellFAB > &scratch) |
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... | |
void | addObserver (MGLevelOpObserver< LevelData< EBCellFAB > > *a_observer) |
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 void | getDivergenceStencil (VoFStencil &a_divStencil, const FaceIndex &a_face, const DataIndex &a_dit, const Real &a_dx, const EBLevelGrid &a_eblg) |
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 void | doLazyRelax (bool a_doLazyRelax) |
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 | getFluxStencil (VoFStencil &a_fluxStencil, const FaceIndex &a_face, 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 () |
void | calculateAlphaWeight () |
void | calculateRelaxationCoefficient () |
Static Protected Attributes | |
static bool | s_turnOffBCs |
static bool | s_doLazyRelax |
Private Member Functions | |
EBViscousTensorOp () | |
weak construction is bad More... | |
EBViscousTensorOp (const EBViscousTensorOp &a_opin) | |
void | operator= (const EBViscousTensorOp &a_opin) |
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.
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::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_acoef0, | ||
const RefCountedPtr< LevelData< EBCellFAB > > & | a_acoef1, | ||
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 a time-dependent a coefficient.
|
virtual |
Destructor.
|
inlineprivate |
weak construction is bad
References MayDay::Error().
|
inlineprivate |
References MayDay::Error().
|
virtual |
for tga to reset stuff
Implements TGAHelmOp< LevelData< EBCellFAB > >.
|
inlinevirtual |
Reimplemented from AMRLevelOp< LevelData< EBCellFAB > >.
References s_step, s_whichLev, and writeEBLevelname().
|
inlinevirtual |
Reimplemented from AMRLevelOp< LevelData< EBCellFAB > >.
References s_step, s_whichLev, and writeEBAMRname().
it's tga's world—we just live in it.
Reimplemented from TGAHelmOp< LevelData< EBCellFAB > >.
References getCCSigma(), getCellCenteredCoefficients(), getKappaDivSigmaU(), getShearStressDotGradU(), and EBLevelDataOps::kappaWeight().
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().
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().
void EBViscousTensorOp::getCellCenteredCoefficients | ( | LevelData< EBCellFAB > & | a_etaCell, |
LevelData< EBCellFAB > & | a_lambdaCell | ||
) |
Referenced by kappaScale().
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().
|
inlinevirtual |
another tgaism
Reimplemented from TGAHelmOp< LevelData< EBCellFAB > >.
References BoxLayout::dataIterator(), EBLevelGrid::getDBL(), EBLevelDataOps::kappaWeight(), m_acoef, m_eblg, LayoutIterator::ok(), and SpaceDim.
|
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).
a_rhs | The 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.
a leveltgaism
Implements LevelTGAHelmOp< LevelData< EBCellFAB >, EBFluxFAB >.
Referenced by divideByIdentityCoef().
|
virtual |
another leveltgaism
Implements LevelTGAHelmOp< LevelData< EBCellFAB >, EBFluxFAB >.
Referenced by divideByIdentityCoef(), and doLazyRelax().
void EBViscousTensorOp::getVelDotSigma | ( | LevelData< EBFluxFAB > & | a_velDotSigma, |
const LevelData< EBFluxFAB > & | a_vel, | ||
const LevelData< EBFluxFAB > & | a_sigma | ||
) |
Referenced by divideByIdentityCoef().
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 | ||
) |
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).
a_oldTime | The time at the beginning of the current step. |
a_mu | The fraction of the current step that has elapsed. |
a_dt | The size of the current step. |
Reimplemented from TGAHelmOp< LevelData< EBCellFAB > >.
Referenced by divideByIdentityCoef().
|
inlinevirtual |
apply operator without any boundary or coarse-fine boundary conditions and no finer level
Implements TGAHelmOp< LevelData< EBCellFAB > >.
References AMRNorm(), AMROperator(), AMROperatorNC(), AMROperatorNF(), AMRProlong(), AMRResidual(), AMRResidualNC(), AMRResidualNF(), AMRRestrict(), AMRUpdateResidual(), applyOp(), assign(), axby(), create(), createCoarsened(), createCoarser(), dotProduct(), finerOperatorChanged(), getDivergenceStencil(), getGradientStencil(), incr(), norm(), preCond(), prolongIncrement(), refToCoarser(), refToFiner(), relax(), residual(), restrictResidual(), s_turnOffBCs, scale(), setToZero(), and setVal().
|
virtual |
This is called on multigrid operators when their AMR operators are altered.
Reimplemented from MGLevelOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
virtual |
a_residual = a_rhs - L(a_phiFine, a_phi) no coaser AMR level
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
virtual |
apply AMR operator no coaser AMR level
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
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().
|
virtual |
Given the current state of the residual the correction, apply your preconditioner to a_cor.
Implements LinearOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
virtual |
this is the linearop function. CFBC is set to homogeneous. phic is null
Implements LinearOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
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().
|
virtual |
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
virtual |
Reimplemented from AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
virtual |
Set a_lhs equal to a_rhs.
Implements LinearOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
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().
|
virtual |
Increment by scaled amount (a_lhs += a_scale*a_x).
Implements LinearOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
virtual |
Set input to a scaled sum (a_lhs = a_a*a_x + a_b*a_y).
Implements LinearOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
virtual |
Multiply the input by a given scale (a_lhs *= a_scale).
Implements LinearOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
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().
Set a_lhs to zero.
Implements LinearOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
virtual |
Referenced by applyOpNoBoundary().
|
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().
|
virtual |
Use your relaxtion 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().
|
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().
|
virtual |
Correct the fine solution based on coarse correction: a_phiThisLevel += I[2h->h] (a_correctCoarse)
Implements MGLevelOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
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().
|
virtual |
Refinement ratio between this level and coarser level. Returns 1 when there are no coarser AMRLevelOp objects
Referenced by applyOpNoBoundary().
|
virtual |
a_residual = a_rhs - L(a_phi, a_phiFine, a_phiCoarse)
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
virtual |
a_residual = a_rhs - L(a_phi, a_phiCoarse)
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
virtual |
a_residual = a_rhs - L(a_phi, a_phiFine, a_phiCoarse)
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
virtual |
a_residual = a_rhs - L(a_phi, a_phiCoarse)
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
virtual |
a_resCoarse = I[h-2h] (a_residual - L(a_correction, a_coarseCorrection))
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
virtual |
a_correction += I[2h->h](a_coarseCorrection)
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
virtual |
a_residual = a_residual - L(a_correction, a_coarseCorrection)
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by applyOpNoBoundary().
|
static |
a_residual = a_residual - L(a_correction, a_coarseCorrection)
Referenced by applyOpNoBoundary(), and doLazyRelax().
|
static |
a_residual = a_residual - L(a_correction, a_coarseCorrection)
Referenced by applyOpNoBoundary(), and doLazyRelax().
|
inlinestatic |
References applyOpIrregular(), averageToCells(), cellGrad(), cfinterp(), defineStencils(), getDivergenceStencil(), getFaceCenteredFluxStencil(), getFlux(), getFluxStencil(), getGradientStencil(), getVoFStencil(), gsrbColor(), homogeneousCFInterp(), incrementFRCoar(), incrementFRFine(), incrOpRegularDir(), reflux(), and s_doLazyRelax.
void EBViscousTensorOp::defineStencils | ( | ) |
(Re)define the stencils for the given coefficients.
Referenced by doLazyRelax().
|
protected |
Referenced by doLazyRelax().
|
protected |
Referenced by doLazyRelax().
|
protected |
|
protected |
Referenced by doLazyRelax().
|
protected |
Referenced by doLazyRelax().
Referenced by doLazyRelax().
|
protected |
Referenced by doLazyRelax().
|
protected |
Referenced by doLazyRelax().
|
protected |
|
protected |
|
protected |
Referenced by doLazyRelax().
|
protected |
Referenced by doLazyRelax().
|
protected |
Referenced by doLazyRelax().
|
protected |
Referenced by doLazyRelax().
|
protected |
Referenced by doLazyRelax().
|
protected |
Referenced by doLazyRelax().
|
protected |
|
protected |
|
protected |
|
inlineprivate |
References MayDay::Error().
|
static |
Referenced by outputAMR(), and outputLevel().
|
static |
Referenced by outputAMR(), and outputLevel().
|
staticprotected |
Referenced by applyOpNoBoundary().
|
staticprotected |
Referenced by doLazyRelax().
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
Referenced by diagonalScale(), and divideByIdentityCoef().
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
"Current" (time-interpolated) value of the a coefficient. For a time-independent a coefficient, this is where the coefficient lives.
Referenced by diagonalScale(), and divideByIdentityCoef().
|
protected |
The a coefficient at the beginning of the time step. This is NULL if the conductivity operator is time independent.
|
protected |
The a coefficient at the end of the time step. This is NULL if the conductivity operator is time independent.
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |