Chombo + EB + MF
3.2
|
#include <EBAMRPoissonOp.H>
Classes | |
class | StencilIndex |
class | StencilIndexComparator |
Public Member Functions | |
void | getVoFStencil (VoFStencil &a_stenc, const VolIndex &a_vof, const DataIndex &a_dit) const |
virtual void | dumpAMR (Vector< LevelData< EBCellFAB > *> &a_data, string name) |
virtual void | dumpLevel (LevelData< EBCellFAB > &a_data, string name) |
int | getPetscMatrix (Mat &a_petsc_mat) |
Fill a_petsc_mat with a matrix that describes the operator stencil. More... | |
int | getLevelDataFromPetscVector (LevelData< EBCellFAB > &a_data, const Vec &a_petsc_vec_real, const Vec &a_petsc_vec_imag) |
put the real components into comp 0, the imaginary ones into comp 1 More... | |
virtual void | setAlphaAndBeta (const Real &a_alpha, const Real &a_beta) |
virtual void | diagonalScale (LevelData< EBCellFAB > &a_rhs, bool a_kappaWeighted=true) |
virtual void | divideByIdentityCoef (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... | |
Real | dx () const |
returns m_dx, such function is required by some LinearSolvers More... | |
virtual void | fillGrad (const LevelData< EBCellFAB > &a_phi) |
a leveltgaism More... | |
void | dumpStencilMatrix () |
void | getDomainFluxStencil (VoFStencil &a_stencil, const VolIndex &a_vof, const int a_comp, const DataIndex &a_dit) |
void | dumpReferenceStencilMatrix () |
virtual void | getFlux (EBFluxFAB &a_flux, const LevelData< EBCellFAB > &a_data, const Box &a_grid, const DataIndex &a_dit, Real a_scale) |
EBLevelGrid | getEBLG () |
EBLevelGrid | getEBLGCoarMG () |
virtual void | setTime (Real a_time) |
virtual | ~EBAMRPoissonOp () |
EBAMRPoissonOp () | |
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) |
EBAMRPoissonOp (const EBLevelGrid &a_eblgFine, const EBLevelGrid &a_eblg, const EBLevelGrid &a_eblgCoar, const EBLevelGrid &a_eblgCoarMG, const RefCountedPtr< EBQuadCFInterp > &a_quadCFI, const RefCountedPtr< BaseDomainBC > &a_domainBC, const RefCountedPtr< BaseEBBC > &a_ebBC, const RealVect &a_dx, const RealVect &a_dxCoar, const RealVect &a_origin, 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 int &a_numPreCondIters, const int &a_relaxType, const Real &a_alpha, const Real &a_beta, const IntVect &a_ghostCellsPhi, const IntVect &a_ghostCellsRHS, int a_testRef=2) | |
virtual void | residual (LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_rhs, bool a_homogeneousPhysBC=false) |
virtual void | getOpMatrix (const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_rhs) |
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, const LevelData< EBCellFAB > *const a_phiCoar, const bool &a_homogeneousPhysBC, const bool &a_homogeneousCFBC, const LevelData< BaseIVFAB< Real > > *const a_ebFluxBCLD) |
virtual void | applyOp (LevelData< EBCellFAB > &a_opPhi, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > *const a_phiCoar, DataIterator &a_dit, const bool &a_homogeneousPhysBC, const bool &a_homogeneousCFBC, const LevelData< BaseIVFAB< Real > > *const a_ebFluxBCLD) |
virtual void | GSColorAllRegular (BaseFab< Real > &a_phi, const BaseFab< Real > &a_rhs, const int &a_icolor, const Real &a_weight, const bool &a_homogeneousPhysBC, const DataIndex &a_dit) |
virtual void | GSColorAllIrregular (EBCellFAB &a_phi, const EBCellFAB &a_rhs, const int &a_icolor, const bool &a_homogeneousPhysBC, const DataIndex &a_dit) |
virtual void | GSColorAllRegularClone (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiOld, const LevelData< EBCellFAB > &a_rhs, const int &a_icolor, const Real &a_weight, const bool &a_homogeneousPhysBC) |
virtual void | GSColorAllIrregularClone (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiOld, const LevelData< EBCellFAB > &a_rhs, const int &a_icolor, const bool &a_homogeneousPhysBC) |
void | mvApplyDomainFlux (BaseFab< Real > &a_phiFAB, const Box &a_grid, const DataIndex &a_dit) |
apply domainflux in multivariable mode More... | |
void | mvBetaLaplacianGrid (EBCellFAB &a_lph, const EBCellFAB &a_phi, const DataIndex &a_dit) |
virtual void | applyOp (LevelData< EBCellFAB > &a_opPhi, const LevelData< EBCellFAB > &a_phi, bool a_homogeneousPhysBC) |
void | applyOpNoCFBCs (LevelData< EBCellFAB > &a_opPhi, const LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > *const a_phiCoar, DataIterator &a_dit, const bool &a_homogeneousPhysBC, const bool &a_homogeneousCFBC, const LevelData< BaseIVFAB< Real > > *const a_ebFluxBCLD) |
no exchange of cf interp 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) |
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) |
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 void | assignLocal (LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs) |
copier definition was killing us. More... | |
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 | AMRUpdateResidual (LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_correction, const LevelData< EBCellFAB > &a_coarseCorrection, const LevelData< BaseIVFAB< Real > > *const a_ebFluxBCLD) |
void | reflux (LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi, AMRLevelOp< LevelData< EBCellFAB > > *a_finerOp) |
void | fast_reflux (LevelData< EBCellFAB > &a_residual, const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi, AMRLevelOp< LevelData< EBCellFAB > > *a_finerOp) |
void | setEBBC (const RefCountedPtr< BaseEBBC > &a_ebBC) |
void | slowGSRBColor (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_res, int a_iredblack) |
void | levelSlowRelax (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_rhs) |
void | levelJacobi (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_rhs, int a_iterations) |
void | levelMultiColorGS (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_rhs) |
void | levelMultiColorGS (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_resid, const IntVect &color) |
void | colorGS (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_rhs, const int &a_icolor) |
void | levelGSRB (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_rhs) |
void | levelGSRB (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_rhs, const int a_color) |
void | levelMultiColorGSClone (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_rhs) |
void | colorGSClone (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > &a_phiOld, const LevelData< EBCellFAB > &a_rhs, const int &a_icolor) |
void | getDivFStencil (VoFStencil &a_vofStencil, const VolIndex &a_vof, const DataIndex &a_dit, bool doFaceInterp) const |
void | getVoFStencil (LayoutData< BaseIVFAB< VoFStencil > > const *&a_vofStencil) |
void | getAlphaDiagWeight (LayoutData< BaseIVFAB< Real > > const *&a_alphaDiagWeight) |
void | getAlphaBeta (Real &a_alpha, Real &a_beta) |
const RefCountedPtr< BaseDomainBC > | getDomainBC () |
void | setListValue (const LevelData< EBCellFAB > &a_data, Real a_value) |
void | setListValue (EBCellFAB &a_data, const Vector< VolIndex > &a_setList, Real a_value) |
void | setRhsSetList (const LayoutData< Vector< VolIndex > > &a_list) |
LayoutData< Vector< VolIndex > > & | getRhsSetList () |
const LayoutData< Vector< VolIndex > > & | getListValue () |
void | applyCFBCs (LevelData< EBCellFAB > &a_phi, const LevelData< EBCellFAB > *const a_phiCoarse, bool a_homogeneousCFBC, bool a_doOnlyRegularInterp=false) |
void | getInvDiagRHS (LevelData< EBCellFAB > &a_lhs, const LevelData< EBCellFAB > &a_rhs) |
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_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 > > | |
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 | 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... | |
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 | mDotProduct (const LevelData< EBCellFAB > &a_1, const int a_sz, const LevelData< EBCellFAB > a_2[], Real a_mdots[]) |
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 int | getMatrixIndexingLD (LevelData< BaseEBCellFAB< int > > &a_gids, int &a_data, const EBLevelGrid &a_eblg, const IntVect &a_ghostCellsPhi, const bool &a_hasCoar) |
static void | getAggregatedLayout (DisjointBoxLayout &a_dblCoar, const ProblemDomain &a_domainCoar, const EBIndexSpace *const a_ebisPtr, const int &a_maxBoxSize) |
static bool | getCoarserLayouts (DisjointBoxLayout &a_dblCoar, ProblemDomain &a_domainCoar, const DisjointBoxLayout &a_dblFine, const EBISLayout &a_ebislFine, const ProblemDomain &a_domainFine, int a_refToCoar, const EBIndexSpace *a_ebisPtr, int a_maxBoxSize, bool &a_layoutChanged, int a_testRef=2) |
static Real | staticMaxNorm (const LevelData< EBCellFAB > &a_rhs, const EBLevelGrid &a_eblg) |
static void | setOperatorTime (Real a_time) |
static void | doLazyRelax (bool a_doLazyRelax) |
static void | doEBEllipticLoadBalance (bool a_doEBEllipticLoadBalance) |
static void | areaFracWeighted (bool a_areaFracWeighted) |
static void | getDivFStencil (VoFStencil &a_vofStencil, const VolIndex &a_vof, const EBISBox &a_ebisBox, const IntVectSet &a_cfivs, const RealVect &a_dx, bool doFaceInterp=true) |
static void | getFluxStencil (VoFStencil &a_fluxStencil, const FaceIndex &a_face, const EBISBox &a_ebisBox, const IntVectSet &a_cfivs, const RealVect &a_dx, bool a_doFaceInterp) |
static void | getFaceCenteredFluxStencil (VoFStencil &a_fluxStencil, const FaceIndex &a_face, const RealVect &a_dx) |
Static Public Attributes | |
static int | s_numComps |
static int | s_whichComp |
Protected Member Functions | |
void | defineMGObjects (const EBLevelGrid &a_eblgCoarMG) |
void | defineWithCoarser (const EBLevelGrid &a_eblgCoar, const int &a_refToCoar) |
void | defineWithFiner (const EBLevelGrid &a_eblgFine, const int &a_refToFine) |
void | defineStencils () |
void | defineEBCFStencils () |
Static Protected Attributes | |
static bool | s_turnOffBCs |
static bool | s_doEBEllipticLoadBalance |
static bool | s_areaFracWeighted |
static Real | s_time |
static bool | s_doLazyRelax |
static bool | s_doInconsistentRelax |
static bool | s_doTrimEdges |
static bool | s_doSetListValueResid |
Private Member Functions | |
void | fast_incrementFRCoar (const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi) |
void | fast_incrementFRFine (const LevelData< EBCellFAB > &a_phiFine, const LevelData< EBCellFAB > &a_phi, AMRLevelOp< LevelData< EBCellFAB > > *a_finerOp) |
void | getFlux (EBFaceFAB &a_flux, const EBCellFAB &a_phi, const Box &a_ghostedBox, const Box &a_fabBox, const ProblemDomain &a_domainBox, const EBISBox &a_ebisBox, const RealVect &a_dx, const int &a_idir) |
this one gets called by base level tga (called by the public getFlux) More... | |
void | getFluxEBCF (EBFaceFAB &a_flux, const EBCellFAB &a_phi, const Box &a_ghostedBox, Vector< FaceIndex > &a_faceitEBCF, Vector< VoFStencil > &a_ebcfsten, const RealVect &a_dx) |
this one is internal (called by refluxing) More... | |
void | getFluxRegO (EBFaceFAB &a_flux, const EBCellFAB &a_phi, const Box &a_ghostedBox, const RealVect &a_dx) |
this one is internal (called by refluxing) More... | |
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) |
Real | getRadius (const FaceIndex &a_face, const RealVect &a_centroid) |
void | applyOpRegularAllDirs (Box *a_loBox, Box *a_hiBox, int *a_hasLo, int *a_hasHi, Box &a_curOpPhiBox, Box &a_curPhiBox, int a_nComps, BaseFab< Real > &a_curOpPhiFAB, const BaseFab< Real > &a_curPhiFAB, bool a_homogeneousPhysBC, const DataIndex &a_dit, const Real &a_beta) |
void | applyDomainFlux (Box *a_loBox, Box *a_hiBox, int *a_hasLo, int *a_hasHi, Box &a_curPhiBox, int a_nComps, BaseFab< Real > &a_phiFAB, bool a_homogeneousPhysBC, const DataIndex &a_dit, const Real &a_beta) |
EBAMRPoissonOp (const EBAMRPoissonOp &a_opin) | |
void | operator= (const EBAMRPoissonOp &a_opin) |
Operator to solve (alpha + beta lapl)phi = rhs. This follows the AMRLevelOp interface.
|
virtual |
Referenced by setTime().
EBAMRPoissonOp::EBAMRPoissonOp | ( | ) |
Referenced by setTime().
EBAMRPoissonOp::EBAMRPoissonOp | ( | const EBLevelGrid & | a_eblgFine, |
const EBLevelGrid & | a_eblg, | ||
const EBLevelGrid & | a_eblgCoar, | ||
const EBLevelGrid & | a_eblgCoarMG, | ||
const RefCountedPtr< EBQuadCFInterp > & | a_quadCFI, | ||
const RefCountedPtr< BaseDomainBC > & | a_domainBC, | ||
const RefCountedPtr< BaseEBBC > & | a_ebBC, | ||
const RealVect & | a_dx, | ||
const RealVect & | a_dxCoar, | ||
const RealVect & | a_origin, | ||
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 int & | a_numPreCondIters, | ||
const int & | a_relaxType, | ||
const Real & | a_alpha, | ||
const Real & | a_beta, | ||
const IntVect & | a_ghostCellsPhi, | ||
const IntVect & | a_ghostCellsRHS, | ||
int | a_testRef = 2 |
||
) |
If you are approaching this operator from this interface, consider backing away and using EBAMRPoissonOpFactory 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_alpha: coefficent of identity \ a_beta: coefficient of laplacian.\ 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.
|
inlineprivate |
References MayDay::Error().
void EBAMRPoissonOp::getVoFStencil | ( | VoFStencil & | a_stenc, |
const VolIndex & | a_vof, | ||
const DataIndex & | a_dit | ||
) | const |
|
inlinevirtual |
Reimplemented from AMRLevelOp< LevelData< EBCellFAB > >.
References writeEBAMRname().
|
inlinevirtual |
Reimplemented from AMRLevelOp< LevelData< EBCellFAB > >.
References diagonalScale(), getAggregatedLayout(), getCoarserLayouts(), getLevelDataFromPetscVector(), getMatrixIndexingLD(), getPetscMatrix(), setAlphaAndBeta(), staticMaxNorm(), and writeEBLevelname().
int EBAMRPoissonOp::getPetscMatrix | ( | Mat & | a_petsc_mat | ) |
Fill a_petsc_mat with a matrix that describes the operator stencil.
Referenced by dumpLevel().
int EBAMRPoissonOp::getLevelDataFromPetscVector | ( | LevelData< EBCellFAB > & | a_data, |
const Vec & | a_petsc_vec_real, | ||
const Vec & | a_petsc_vec_imag | ||
) |
put the real components into comp 0, the imaginary ones into comp 1
Referenced by dumpLevel().
|
static |
get the matrix indicies associated with every point on the level Not really going to work if level does not cover the domain
Referenced by dumpLevel().
|
static |
Referenced by dumpLevel().
|
static |
version that does not fill ebislcoar
Referenced by dumpLevel().
|
static |
Referenced by dumpLevel().
Sets the scaling constants in the Helmholtz equation.
a_alpha | The scaling constant that multiplies the identity term. |
a_beta | The scaling constant that multiplies the derivative term. |
Implements TGAHelmOp< LevelData< EBCellFAB > >.
Referenced by dumpLevel().
|
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.
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 > >.
Referenced by dumpLevel(), and kappaScale().
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 > >.
for eb only. kappa weight the rhs but do not multiply by the identity coefficient
Reimplemented from TGAHelmOp< LevelData< EBCellFAB > >.
References diagonalScale().
|
inlinevirtual |
returns m_dx, such function is required by some LinearSolvers
Reimplemented from LinearOp< LevelData< EBCellFAB > >.
References m_dx.
a leveltgaism
Implements LevelTGAHelmOp< LevelData< EBCellFAB >, EBFluxFAB >.
void EBAMRPoissonOp::dumpStencilMatrix | ( | ) |
Referenced by EBAMRPoissonOp::StencilIndexComparator::operator()().
void EBAMRPoissonOp::getDomainFluxStencil | ( | VoFStencil & | a_stencil, |
const VolIndex & | a_vof, | ||
const int | a_comp, | ||
const DataIndex & | a_dit | ||
) |
Referenced by EBAMRPoissonOp::StencilIndexComparator::operator()().
void EBAMRPoissonOp::dumpReferenceStencilMatrix | ( | ) |
Referenced by EBAMRPoissonOp::StencilIndexComparator::operator()().
|
inlinevirtual |
Implements LevelTGAHelmOp< LevelData< EBCellFAB >, EBFluxFAB >.
References EBLevelGrid::getDomain(), EBLevelGrid::getEBISL(), Box::grow(), m_dx, m_eblg, and SpaceDim.
|
inline |
References m_eblg.
|
inline |
References m_eblgCoarMG.
|
inlinestatic |
|
inlinevirtual |
Sets the time-dependent state of the operator. The default implementation does nothing and is appropriate for time-independent operators.
a_time | The time to be used to update the time-dependent operator. |
Reimplemented from TGAHelmOp< LevelData< EBCellFAB > >.
References AMRNorm(), AMROperator(), AMROperatorNC(), AMROperatorNF(), AMRProlong(), AMRResidual(), AMRResidualNC(), AMRResidualNF(), AMRRestrict(), AMRUpdateResidual(), applyOp(), applyOpNoBoundary(), applyOpNoCFBCs(), areaFracWeighted(), assign(), assignCopier(), assignLocal(), axby(), buildCopier(), colorGS(), colorGSClone(), create(), createCoarsened(), createCoarser(), doEBEllipticLoadBalance(), doLazyRelax(), dotProduct(), EBAMRPoissonOp(), fast_reflux(), getDivFStencil(), getOpMatrix(), GSColorAllIrregular(), GSColorAllIrregularClone(), GSColorAllRegular(), GSColorAllRegularClone(), incr(), levelGSRB(), levelJacobi(), levelMultiColorGS(), levelMultiColorGSClone(), levelSlowRelax(), localMaxNorm(), mvApplyDomainFlux(), mvBetaLaplacianGrid(), norm(), preCond(), prolongIncrement(), reflux(), refToCoarser(), refToFiner(), relax(), residual(), restrictResidual(), scale(), setEBBC(), setOperatorTime(), setToZero(), setVal(), slowGSRBColor(), and ~EBAMRPoissonOp().
|
virtual |
a_residual = a_rhs - L(a_phiFine, a_phi) no coaser AMR level
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
apply AMR operator no coaser AMR level
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
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 setTime().
|
virtual |
For debugging purpose, get the matrix of the op (could be change from time to time due to EBBC update) by calculating 0 - L(phi(i)) for i=0,1,...,n, each time phi(i)=1 for i-th cell and 0 for all other cells
Referenced by setTime().
|
virtual |
Given the current state of the residual the correction, apply your preconditioner to a_cor.
Implements LinearOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
This function assumes that coarse-fine boundary condtions have been dealt with.
Referenced by setTime().
|
virtual |
virtual function called by LevelTGA
Implements TGAHelmOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
get EB flux from leveldata if we have one. otherwise use m_ebBC ignore input data in the case of homogeneous Phys BC
|
virtual |
|
virtual |
Referenced by setTime().
|
virtual |
Referenced by setTime().
|
virtual |
Referenced by setTime().
|
virtual |
Referenced by setTime().
void EBAMRPoissonOp::mvApplyDomainFlux | ( | BaseFab< Real > & | a_phiFAB, |
const Box & | a_grid, | ||
const DataIndex & | a_dit | ||
) |
apply domainflux in multivariable mode
Referenced by setTime().
void EBAMRPoissonOp::mvBetaLaplacianGrid | ( | EBCellFAB & | a_lph, |
const EBCellFAB & | a_phi, | ||
const DataIndex & | a_dit | ||
) |
for taking the multi-variable laplacian (think source term in viscous flow) evaluates beta*laplacian (with inhomogeneous bcs).
Referenced by setTime().
|
virtual |
this is the linearop function. CFBC is set to homogeneous. phic is null
Implements LinearOp< LevelData< EBCellFAB > >.
void EBAMRPoissonOp::applyOpNoCFBCs | ( | LevelData< EBCellFAB > & | a_opPhi, |
const LevelData< EBCellFAB > & | a_phi, | ||
const LevelData< EBCellFAB > *const | a_phiCoar, | ||
DataIterator & | a_dit, | ||
const bool & | a_homogeneousPhysBC, | ||
const bool & | a_homogeneousCFBC, | ||
const LevelData< BaseIVFAB< Real > > *const | a_ebFluxBCLD | ||
) |
no exchange of cf interp
Referenced by setTime().
|
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 setTime().
|
virtual |
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
Reimplemented from AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
Reimplemented from AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
Reimplemented from AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
copier definition was killing us.
Reimplemented from LinearOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
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 setTime().
|
virtual |
Increment by scaled amount (a_lhs += a_scale*a_x).
Implements LinearOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
Set input to a scaled sum (a_lhs = a_a*a_x + a_b*a_y).
Implements LinearOp< LevelData< EBCellFAB > >.
Referenced by setTime().
Multiply the input by a given scale (a_lhs *= a_scale).
Implements LinearOp< LevelData< EBCellFAB > >.
Referenced by setTime().
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 setTime().
Reimplemented from AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
Referenced by setTime().
|
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 setTime().
|
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 setTime().
|
virtual |
Calculate restricted residual: a_resCoarse[2h] = I[h->2h] (a_rhsFine[h] - L[h](a_phiFine[h]))
Implements MGLevelOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
Correct the fine solution based on coarse correction: a_phiThisLevel += I[2h->h] (a_correctCoarse)
Implements MGLevelOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
Refinement ratio between this level and coarser level. Returns 1 when there are no coarser AMRLevelOp objects
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
Refinement ratio between this level and coarser level. Returns 1 when there are no coarser AMRLevelOp objects
Referenced by setTime().
|
virtual |
a_residual = a_rhs - L(a_phi, a_phiFine, a_phiCoarse)
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
a_residual = a_rhs - L(a_phi, a_phiCoarse)
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
a_residual = a_rhs - L(a_phi, a_phiFine, a_phiCoarse)
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
a_residual = a_rhs - L(a_phi, a_phiCoarse)
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
a_resCoarse = I[h-2h] (a_residual - L(a_correction, a_coarseCorrection))
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
a_correction += I[2h->h](a_coarseCorrection)
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by setTime().
|
virtual |
a_residual = a_residual - L(a_correction, a_coarseCorrection)
Implements AMRLevelOp< LevelData< EBCellFAB > >.
Referenced by setTime().
void EBAMRPoissonOp::AMRUpdateResidual | ( | LevelData< EBCellFAB > & | a_residual, |
const LevelData< EBCellFAB > & | a_correction, | ||
const LevelData< EBCellFAB > & | a_coarseCorrection, | ||
const LevelData< BaseIVFAB< Real > > *const | a_ebFluxBCLD | ||
) |
a_residual = a_residual - L(a_correction, a_coarseCorrection) used in multifluid
void EBAMRPoissonOp::reflux | ( | LevelData< EBCellFAB > & | a_residual, |
const LevelData< EBCellFAB > & | a_phiFine, | ||
const LevelData< EBCellFAB > & | a_phi, | ||
AMRLevelOp< LevelData< EBCellFAB > > * | a_finerOp | ||
) |
Referenced by setTime().
void EBAMRPoissonOp::fast_reflux | ( | LevelData< EBCellFAB > & | a_residual, |
const LevelData< EBCellFAB > & | a_phiFine, | ||
const LevelData< EBCellFAB > & | a_phi, | ||
AMRLevelOp< LevelData< EBCellFAB > > * | a_finerOp | ||
) |
Referenced by setTime().
void EBAMRPoissonOp::setEBBC | ( | const RefCountedPtr< BaseEBBC > & | a_ebBC | ) |
Referenced by setTime().
void EBAMRPoissonOp::slowGSRBColor | ( | LevelData< EBCellFAB > & | a_phi, |
const LevelData< EBCellFAB > & | a_res, | ||
int | a_iredblack | ||
) |
Referenced by setTime().
void EBAMRPoissonOp::levelSlowRelax | ( | LevelData< EBCellFAB > & | a_phi, |
const LevelData< EBCellFAB > & | a_rhs | ||
) |
Referenced by setTime().
void EBAMRPoissonOp::levelJacobi | ( | LevelData< EBCellFAB > & | a_phi, |
const LevelData< EBCellFAB > & | a_rhs, | ||
int | a_iterations | ||
) |
Referenced by setTime().
void EBAMRPoissonOp::levelMultiColorGS | ( | LevelData< EBCellFAB > & | a_phi, |
const LevelData< EBCellFAB > & | a_rhs | ||
) |
Referenced by setTime().
void EBAMRPoissonOp::levelMultiColorGS | ( | LevelData< EBCellFAB > & | a_phi, |
const LevelData< EBCellFAB > & | a_resid, | ||
const IntVect & | color | ||
) |
void EBAMRPoissonOp::colorGS | ( | LevelData< EBCellFAB > & | a_phi, |
const LevelData< EBCellFAB > & | a_rhs, | ||
const int & | a_icolor | ||
) |
Referenced by setTime().
void EBAMRPoissonOp::levelGSRB | ( | LevelData< EBCellFAB > & | a_phi, |
const LevelData< EBCellFAB > & | a_rhs | ||
) |
Referenced by setTime().
void EBAMRPoissonOp::levelGSRB | ( | LevelData< EBCellFAB > & | a_phi, |
const LevelData< EBCellFAB > & | a_rhs, | ||
const int | a_color | ||
) |
void EBAMRPoissonOp::levelMultiColorGSClone | ( | LevelData< EBCellFAB > & | a_phi, |
const LevelData< EBCellFAB > & | a_rhs | ||
) |
Referenced by setTime().
void EBAMRPoissonOp::colorGSClone | ( | LevelData< EBCellFAB > & | a_phi, |
const LevelData< EBCellFAB > & | a_phiOld, | ||
const LevelData< EBCellFAB > & | a_rhs, | ||
const int & | a_icolor | ||
) |
Referenced by setTime().
|
static |
Referenced by setTime().
|
static |
Referenced by setTime().
|
static |
Referenced by setTime().
|
static |
Referenced by setTime().
void EBAMRPoissonOp::getDivFStencil | ( | VoFStencil & | a_vofStencil, |
const VolIndex & | a_vof, | ||
const DataIndex & | a_dit, | ||
bool | doFaceInterp | ||
) | const |
|
inline |
References m_opStencil.
|
inline |
References m_alphaDiagWeight.
|
inline |
References m_domainBC, setListValue(), and setRhsSetList().
Referenced by getDomainBC().
void EBAMRPoissonOp::setListValue | ( | EBCellFAB & | a_data, |
const Vector< VolIndex > & | a_setList, | ||
Real | a_value | ||
) |
void EBAMRPoissonOp::setRhsSetList | ( | const LayoutData< Vector< VolIndex > > & | a_list | ) |
Referenced by getDomainBC().
|
inline |
References m_rhsSetList.
|
inline |
References applyCFBCs(), getFaceCenteredFluxStencil(), getFluxStencil(), getInvDiagRHS(), and m_rhsSetList.
|
static |
Referenced by getListValue().
|
static |
Referenced by getListValue().
void EBAMRPoissonOp::applyCFBCs | ( | LevelData< EBCellFAB > & | a_phi, |
const LevelData< EBCellFAB > *const | a_phiCoarse, | ||
bool | a_homogeneousCFBC, | ||
bool | a_doOnlyRegularInterp = false |
||
) |
Referenced by getListValue().
void EBAMRPoissonOp::getInvDiagRHS | ( | LevelData< EBCellFAB > & | a_lhs, |
const LevelData< EBCellFAB > & | a_rhs | ||
) |
This function computes: a_lhs = (1/diagonal)*a_rhs It is used to initialize the preconditioner, and by MFPoissonOp::levelJacobi. Consider using one of the level Gauss-Seidel methods instead of monkeying with this.
Referenced by getListValue().
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
private |
|
private |
|
private |
this one gets called by base level tga (called by the public getFlux)
|
private |
this one is internal (called by refluxing)
|
private |
this one is internal (called by refluxing)
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
applyOp() on the regular cells for all directions opPhi comes in holding alpha*phi. this adds on beta*lapl phi
|
private |
|
inlineprivate |
References MayDay::Error().
|
static |
|
static |
|
staticprotected |
|
staticprotected |
|
staticprotected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
Referenced by getEBLGCoarMG().
|
protected |
|
protected |
Referenced by getDomainBC().
|
protected |
|
protected |
Referenced by getListValue(), and getRhsSetList().
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
Referenced by getAlphaBeta().
|
protected |
|
protected |
Referenced by getAlphaBeta().
|
protected |
|
staticprotected |
Referenced by setOperatorTime().
|
staticprotected |
|
staticprotected |
|
staticprotected |
|
staticprotected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
Referenced by getVoFStencil().
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
Referenced by getAlphaDiagWeight().
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |