Chombo + EB
3.2
|
#include <AMRMultiGrid.H>
Public Member Functions | |
AMRMultiGrid () | |
virtual | ~AMRMultiGrid () |
void | outputAMR (Vector< T * > &a_data, string &a_name, int a_lmax, int a_lbase) |
virtual void | define (const ProblemDomain &a_coarseDomain, AMRLevelOpFactory< T > &a_factory, LinearSolver< T > *a_bottomSolver, int a_numLevels) |
void | addInspector (RefCountedPtr< AMRMultiGridInspector< T > > &a_inspector) |
virtual void | solve (Vector< T * > &a_phi, const Vector< T * > &a_rhs, int l_max, int l_base, bool a_zeroPhi=true, bool forceHomogeneous=false) |
virtual void | solveNoInit (Vector< T * > &a_phi, const Vector< T * > &a_rhs, int l_max, int l_base, bool a_zeroPhi=true, bool forceHomogeneous=false) |
virtual void | solveNoInitResid (Vector< T * > &a_phi, Vector< T * > &a_finalResid, const Vector< T * > &a_rhs, int l_max, int l_base, bool a_zeroPhi=true, bool forceHomogeneous=false) |
use if you want final residual More... | |
void | relaxOnlyHomogeneous (Vector< T * > &a_phi, const Vector< T * > &a_rhs, int l_max, int l_base) |
virtual void | AMRVCycle (Vector< T * > &a_correction, Vector< T * > &a_residual, int l, int l_max, int l_base) |
void | setMGCycle (int a_numMG) |
virtual void | init (const Vector< T * > &a_phi, const Vector< T * > &a_rhs, int l_max, int l_base) |
void | revert (const Vector< T * > &a_phi, const Vector< T * > &a_rhs, int l_max, int l_base) |
AMRLevelOp< T > & | levelOp (int level) |
Real | computeAMRResidual (Vector< T * > &a_resid, Vector< T * > &a_phi, const Vector< T * > &a_rhs, int l_max, int l_base, bool a_homogeneousBC=false, bool a_computeNorm=true) |
Real | computeAMRResidual (Vector< T * > &a_phi, const Vector< T * > &a_rhs, int l_max, int l_min) |
void | computeAMROperator (Vector< T * > &a_lph, Vector< T * > &a_phi, int l_max, int l_base, bool a_homogeneousBC=false) |
Vector< MGLevelOp< T > * > | getAllOperators () |
Vector< MGLevelOp< T > * > | getOperatorsOp () |
Vector< Vector< MGLevelOp< T > * > > | getOperatorsMG () |
Vector< AMRLevelOp< T > * > & | getAMROperators () |
void | setSolverParameters (const int &a_pre, const int &a_post, const int &a_bottom, const int &a_numMG, const int &a_iterMax, const Real &a_eps, const Real &a_hang, const Real &a_normThresh) |
void | setBottomSolver (int l_max, int l_base) |
set up bottom solver for internal MG solver More... | |
void | setBottomSolverEpsCushion (Real a_bottomSolverEpsCushion) |
void | getInfo () const |
Dump out the state of the solver. AMR levels, MG levels, bottom solver, box counts, etc. More... | |
Public Attributes | |
Real | m_eps |
Real | m_hang |
Real | m_normThresh |
bool | m_solverParamsSet |
int | m_imin |
int | m_iterMax |
int | m_iterMin |
int | m_verbosity |
int | m_exitStatus |
int | m_pre |
int | m_post |
int | m_bottom |
int | m_numMG |
int | m_maxDepth |
max no. of coarsenings – -1 (default) means coarsen as far as possible More... | |
Real | m_convergenceMetric |
Real | m_bottomSolverEpsCushion |
Protected Member Functions | |
void | relax (T &phi, T &R, int depth, int nRelax=2) |
void | computeAMRResidualLevel (Vector< T * > &a_resid, Vector< T * > &a_phi, const Vector< T * > &a_rhs, int l_max, int l_base, int ilev, bool a_homogeneousBC) |
void | clear () |
Protected Attributes | |
Vector< AMRLevelOp< T > * > | m_op |
Vector< MultiGrid< T > * > | m_mg |
Vector< T * > | m_correction |
Vector< T * > | m_residual |
Vector< T * > | m_resC |
Vector< Copier > | m_resCopier |
Vector< Copier > | m_reverseCopier |
NoOpSolver< T > | m_nosolve |
LinearSolver< T > * | m_bottomSolver |
Vector< char > | m_hasInitBeenCalled |
Private Member Functions | |
AMRMultiGrid (const AMRMultiGrid< T > &) | |
AMRMultiGrid & | operator= (const AMRMultiGrid< T > &) |
Private Attributes | |
Vector< RefCountedPtr < AMRMultiGridInspector< T > > > | m_inspectors |
Class to solve elliptic equations using the Martin and Cartwright algorithm.
AMRMultiGrid< T >::AMRMultiGrid | ( | ) |
References AMRMultiGrid< T >::m_solverParamsSet.
|
virtual |
References CH_TIME.
|
private |
void AMRMultiGrid< T >::outputAMR | ( | Vector< T * > & | a_data, |
string & | a_name, | ||
int | a_lmax, | ||
int | a_lbase | ||
) |
References Vector< T >::push_back().
|
virtual |
Define the solver. a_coarseDomain is the index space on the coarsest AMR level. a_factory is the operator factory through which all special information is conveyed. a_bottomSolver is the solver to be used at the termination of multigrid coarsening. It is the client's responsibility to free up the dynamically-allocated memory. a_numLevels is the number of AMR levels.
Reimplemented in AMRFASMultiGrid< T >.
References AMRLevelOpFactory< T >::AMRnewOp(), CH_TIME, ProblemDomain::refine(), and AMRLevelOpFactory< T >::refToFiner().
Referenced by AMRFASMultiGrid< T >::define().
|
inline |
Add an inspector to the list of inspectors maintained by this AMRMultiGrid instance. It will be given the opportunity to record intermediate data.
|
virtual |
Solve L(phi) = rho from l_base to l_max. To solve over all levels, l_base = 0 and l_max = max_level = numLevels-1.
References CH_TIME.
Referenced by BaseLevelHeatSolver< LevelData< FArrayBox >, FluxBox, LevelFluxRegister >::solveHelm().
|
virtual |
same as "solve" except user has taken the reponsibility of having previously called "init" so solver can allocate temporary holders.
Reimplemented in AMRFASMultiGrid< T >.
References CH_TIME, and Vector< T >::size().
|
virtual |
void AMRMultiGrid< T >::relaxOnlyHomogeneous | ( | Vector< T * > & | a_phi, |
const Vector< T * > & | a_rhs, | ||
int | l_max, | ||
int | l_base | ||
) |
References CH_assert, CH_TIME, pout(), and Vector< T >::size().
|
virtual |
References CH_assert.
void AMRMultiGrid< T >::setMGCycle | ( | int | a_numMG | ) |
|
virtual |
Reimplemented in AMRFASMultiGrid< T >.
References AMRLevelOp< T >::buildCopier(), CH_TIME, LinearOp< T >::create(), AMRLevelOp< T >::createCoarsened(), and AMRLevelOp< T >::refToCoarser().
Referenced by AMRFASMultiGrid< T >::init().
void AMRMultiGrid< T >::revert | ( | const Vector< T * > & | a_phi, |
const Vector< T * > & | a_rhs, | ||
int | l_max, | ||
int | l_base | ||
) |
References CH_TIME.
AMRLevelOp< T > & AMRMultiGrid< T >::levelOp | ( | int | level | ) |
Real AMRMultiGrid< T >::computeAMRResidual | ( | Vector< T * > & | a_resid, |
Vector< T * > & | a_phi, | ||
const Vector< T * > & | a_rhs, | ||
int | l_max, | ||
int | l_base, | ||
bool | a_homogeneousBC = false , |
||
bool | a_computeNorm = true |
||
) |
resid = L(phi) - rhs
References CH_TIME, MayDay::Error(), Max(), and MPI_CH_REAL.
Real AMRMultiGrid< T >::computeAMRResidual | ( | Vector< T * > & | a_phi, |
const Vector< T * > & | a_rhs, | ||
int | l_max, | ||
int | l_min | ||
) |
just return the normed value of computeAMRResidual. used for benchmarking
void AMRMultiGrid< T >::computeAMROperator | ( | Vector< T * > & | a_lph, |
Vector< T * > & | a_phi, | ||
int | l_max, | ||
int | l_base, | ||
bool | a_homogeneousBC = false |
||
) |
Vector< MGLevelOp< T > * > AMRMultiGrid< T >::getAllOperators | ( | ) |
For changing coefficients. Use at thy own peril.
References Vector< T >::append(), and Vector< T >::push_back().
Referenced by BaseLevelHeatSolver< LevelData< FArrayBox >, FluxBox, LevelFluxRegister >::resetSolverAlphaAndBeta().
Vector< MGLevelOp< T > * > AMRMultiGrid< T >::getOperatorsOp | ( | ) |
References Vector< T >::push_back().
Vector< Vector< MGLevelOp< T > * > > AMRMultiGrid< T >::getOperatorsMG | ( | ) |
|
inline |
void AMRMultiGrid< T >::setSolverParameters | ( | const int & | a_pre, |
const int & | a_post, | ||
const int & | a_bottom, | ||
const int & | a_numMG, | ||
const int & | a_iterMax, | ||
const Real & | a_eps, | ||
const Real & | a_hang, | ||
const Real & | a_normThresh | ||
) |
Set parameters of the solve. a_pre is the number of smoothings before averaging. a_post is the number of smoothings after averaging. a_bottom is the number of smoothings at the bottom level. a_numMG = 1 for vcycle, =2 for wcycle (use 1). a_itermax is the max number of v cycles. a_hang is the minimum amount of change per vcycle. a_eps is the solution tolerance. a_normThresh is how close to zero eps*resid is allowed to get.
void AMRMultiGrid< T >::setBottomSolver | ( | int | l_max, |
int | l_base | ||
) |
set up bottom solver for internal MG solver
This function is normally called by the solve(...) function. However, it must be called if solve will not be called (in particular, if only the V-cycle is being used)
void AMRMultiGrid< T >::setBottomSolverEpsCushion | ( | Real | a_bottomSolverEpsCushion | ) |
References CH_assert.
void AMRMultiGrid< T >::getInfo | ( | ) | const |
|
protected |
References CH_TIME.
|
protected |
|
protected |
|
private |
Real AMRMultiGrid< T >::m_eps |
Real AMRMultiGrid< T >::m_hang |
Real AMRMultiGrid< T >::m_normThresh |
bool AMRMultiGrid< T >::m_solverParamsSet |
Referenced by AMRMultiGrid< T >::AMRMultiGrid().
int AMRMultiGrid< T >::m_imin |
int AMRMultiGrid< T >::m_iterMax |
int AMRMultiGrid< T >::m_iterMin |
int AMRMultiGrid< T >::m_verbosity |
int AMRMultiGrid< T >::m_exitStatus |
int AMRMultiGrid< T >::m_pre |
int AMRMultiGrid< T >::m_post |
int AMRMultiGrid< T >::m_bottom |
int AMRMultiGrid< T >::m_numMG |
int AMRMultiGrid< T >::m_maxDepth |
max no. of coarsenings – -1 (default) means coarsen as far as possible
If using a value besides the default, need to set it before define function is called
Real AMRMultiGrid< T >::m_convergenceMetric |
Real AMRMultiGrid< T >::m_bottomSolverEpsCushion |
|
protected |
Referenced by AMRMultiGrid< LevelData< T > >::getAMROperators().
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
private |
Referenced by AMRMultiGrid< LevelData< T > >::addInspector().