Main Page   Namespace List   Class Hierarchy   Alphabetical List   Compound List   File List   Compound Members   File Members  

LevelSolver Class Reference

Solves elliptic equations on a level. More...

#include <LevelSolver.H>

Collaboration diagram for LevelSolver:

Collaboration graph
[legend]
List of all members.

Public Methods

 LevelSolver ()
virtual ~LevelSolver ()
 LevelSolver (const DisjointBoxLayout &a_grids, const DisjointBoxLayout *a_baseGrids, const Box &a_domain, Real a_dxLevel, int a_nRefCrse, const LevelOp *const a_opin, int a_ncomp=1)
 LevelSolver (const DisjointBoxLayout &a_grids, const DisjointBoxLayout *a_baseGrids, const ProblemDomain &a_domain, Real a_dxLevel, int a_nRefCrse, const LevelOp *const a_opin, int a_ncomp)
virtual void define (const DisjointBoxLayout &a_grids, const DisjointBoxLayout *a_baseGrids, const Box &a_domain, Real a_dxLevel, int a_nRefCrse, const LevelOp *const a_opin, int a_ncomp=1)
 Full define function.

virtual void define (const DisjointBoxLayout &a_grids, const DisjointBoxLayout *a_baseGrids, const ProblemDomain &a_domain, Real a_dxLevel, int a_nRefCrse, const LevelOp *const a_opin, int a_ncomp=1)
 Full define function.

bool isDefined () const
void clearMemory ()
void setDefaultValues ()
void levelSolveH (LevelData< FArrayBox > &phi, const LevelData< FArrayBox > &rhs, bool initializePhiToZero=true)
void levelSolve (LevelData< FArrayBox > &phi, const LevelData< FArrayBox > *phic, const LevelData< FArrayBox > &rhs, bool initializePhiToZero=true)
void setTolerance (Real a_tolerance)
 set solution tolerance

void setOperatorTolerance (Real a_operatorTolerance)
void setVerbose (bool a_verbose)
void setMaxIter (int a_maxIter)
 set max number multigrid vcycles to use in solution

void setMinIter (int a_minIter)

Protected Attributes

LevelMG m_levelMG
LevelOpm_levelOpPtr
bool m_isDefined
DisjointBoxLayout m_grids
ProblemDomain m_domain
LevelData< FArrayBoxm_resid
LevelData< FArrayBoxm_scratch
LevelData< FArrayBoxm_corr
int m_nRefCrse
Real m_dxLevel
int m_maxIter
int m_minIter
Real m_tolerance
Real m_operatorTolerance
bool m_verbose

Private Methods

void operator= (const LevelSolver &)
 LevelSolver (const LevelSolver &)

Detailed Description

Solves elliptic equations on a level.

LevelSolver Solves elliptic equations on a level using multigrid.


Constructor & Destructor Documentation

LevelSolver::LevelSolver  
 

Default constructor leaves LevelSolver undefined.

virtual LevelSolver::~LevelSolver   [virtual]
 

LevelSolver::LevelSolver const DisjointBoxLayout   a_grids,
const DisjointBoxLayout   a_baseGrids,
const Box   a_domain,
Real    a_dxLevel,
int    a_nRefCrse,
const LevelOp *const    a_opin,
int    a_ncomp = 1
 

Full constructor. Calls full define function.\ {\bf Inputs:} \ a_grids: Grid on which solution will be computed. \ a_baseGrids: pointer to next coarser grids. Should be == NULL if there is no coarser level. \ a_domain: domain of solution level \ a_dxLevel: grid spacing at solution level \ a_nRefCrse: refinement ratio between two levels. Not used if there is no finer level \ a_opin: operator L used in solution of Lphi=rho

LevelSolver::LevelSolver const DisjointBoxLayout   a_grids,
const DisjointBoxLayout   a_baseGrids,
const ProblemDomain   a_domain,
Real    a_dxLevel,
int    a_nRefCrse,
const LevelOp *const    a_opin,
int    a_ncomp
 

Full constructor. Calls full define function.\ {\bf Inputs:} \ a_grids: Grid on which solution will be computed. \ a_baseGrids: pointer to next coarser grids. Should be == NULL if there is no coarser level. \ a_domain: domain of solution level \ a_dxLevel: grid spacing at solution level \ a_nRefCrse: refinement ratio between two levels. Not used if there is no finer level \ a_opin: operator L used in solution of Lphi=rho

LevelSolver::LevelSolver const LevelSolver &    [inline, private]
 


Member Function Documentation

void LevelSolver::clearMemory  
 

virtual void LevelSolver::define const DisjointBoxLayout   a_grids,
const DisjointBoxLayout   a_baseGrids,
const ProblemDomain   a_domain,
Real    a_dxLevel,
int    a_nRefCrse,
const LevelOp *const    a_opin,
int    a_ncomp = 1
[virtual]
 

Full define function.

Full define function.\ {\bf Inputs:} \ a_grids: Grid on which solution will be computed. \ a_baseGrids: pointer to next coarser grids. Should be == NULL if there is no coarser level. \ a_domain: domain of solution level \ a_dxLevel: grid spacing at solution level \ a_nRefCrse: refinement ratio between two levels. Not used if there is no finer level \ a_opin: operator L used in solution of Lphi=rho

virtual void LevelSolver::define const DisjointBoxLayout   a_grids,
const DisjointBoxLayout   a_baseGrids,
const Box   a_domain,
Real    a_dxLevel,
int    a_nRefCrse,
const LevelOp *const    a_opin,
int    a_ncomp = 1
[virtual]
 

Full define function.

Full define function.\ {\bf Inputs:} \ a_grids: Grid on which solution will be computed. \ a_baseGrids: pointer to next coarser grids. Should be == NULL if there is no coarser level. \ a_domain: domain of solution level \ a_dxLevel: grid spacing at solution level \ a_nRefCrse: refinement ratio between two levels. Not used if there is no finer level \ a_opin: operator L used in solution of Lphi=rho

bool LevelSolver::isDefined   const
 

void LevelSolver::levelSolve LevelData< FArrayBox > &    phi,
const LevelData< FArrayBox > *    phic,
const LevelData< FArrayBox > &    rhs,
bool    initializePhiToZero = true
 

do a level solve on this level using multigrid uses inhomogeneous bcs at interfaces\ {\bf Inputs:} \ phi: initial guess at solution. phic: solution at next coarser level. should be == NULL if there is no coarser level. rhs: right-hand side. \ {\bf Outputs:} \ phi: solution.

void LevelSolver::levelSolveH LevelData< FArrayBox > &    phi,
const LevelData< FArrayBox > &    rhs,
bool    initializePhiToZero = true
 

do a level solve on this level using multigrid uses homogeneous bcs at coarse/fine interfaces \ {\bf Inputs:} \ phi: initial guess at solution. rhs: right-hand side. \ {\bf Outputs:} \ phi: solution.

void LevelSolver::operator= const LevelSolver &    [inline, private]
 

void LevelSolver::setDefaultValues  
 

void LevelSolver::setMaxIter int    a_maxIter
 

set max number multigrid vcycles to use in solution

void LevelSolver::setMinIter int    a_minIter
 

Set min number of iterations. Only relevant when residual is not decreasing fast enough to satisfy the "operator tolerance". Default is 4.

void LevelSolver::setOperatorTolerance Real    a_operatorTolerance
 

Set "operator tolerance" of iterative solution. Iteration will stop if (new_residual/old_residual) > 1-operatorTolerance (and at least minIter iterations have been performed). Default is 1.0e-5.

void LevelSolver::setTolerance Real    a_tolerance
 

set solution tolerance

void LevelSolver::setVerbose bool    a_verbose
 


Member Data Documentation

LevelData<FArrayBox> LevelSolver::m_corr [protected]
 

ProblemDomain LevelSolver::m_domain [protected]
 

Real LevelSolver::m_dxLevel [protected]
 

DisjointBoxLayout LevelSolver::m_grids [protected]
 

bool LevelSolver::m_isDefined [protected]
 

LevelMG LevelSolver::m_levelMG [protected]
 

LevelOp* LevelSolver::m_levelOpPtr [protected]
 

int LevelSolver::m_maxIter [protected]
 

int LevelSolver::m_minIter [protected]
 

int LevelSolver::m_nRefCrse [protected]
 

Real LevelSolver::m_operatorTolerance [protected]
 

LevelData<FArrayBox> LevelSolver::m_resid [protected]
 

LevelData<FArrayBox> LevelSolver::m_scratch [protected]
 

Real LevelSolver::m_tolerance [protected]
 

bool LevelSolver::m_verbose [protected]
 


The documentation for this class was generated from the following file:
Generated on Wed Apr 16 14:35:24 2003 for EBChombo by doxygen1.2.16