BISICLES AMR ice sheet model
0.9
|
#include <CGSolver.H>
Public Member Functions | |
CGSolver () | |
virtual | ~CGSolver () |
virtual void | setHomogeneous (bool a_homogeneous) |
virtual void | define (LinearOp< T > *a_op, bool a_homogeneous) |
virtual void | solve (T &a_phi, const T &a_rhs) |
solve the equation. More... | |
virtual void | setConvergenceMetrics (Real a_metric, Real a_tolerance) |
Public Attributes | |
bool | m_homogeneous |
LinearOp< T > * | m_op |
int | m_imax |
int | m_verbosity |
Real | m_eps |
Real | m_reps |
Real | m_convergenceMetric |
Real | m_hang |
int | m_exitStatus |
Real | m_small |
int | m_numRestarts |
int | m_normType |
Elliptic solver using the preconditioned conjugate gradient algorithm.
References CGSolver< T >::m_op.
|
virtual |
define the solver. a_op is the linear operator. a_homogeneous is whether the solver uses homogeneous boundary conditions.
References CGSolver< T >::m_homogeneous, and CGSolver< T >::m_op.
Referenced by CGSolver< T >::setHomogeneous(), and JFNKSolver::solve().
|
virtual |
References CGSolver< T >::m_convergenceMetric, and CGSolver< T >::m_eps.
Referenced by CGSolver< T >::setHomogeneous().
|
inlinevirtual |
|
virtual |
solve the equation.
References CGSolver< T >::m_convergenceMetric, CGSolver< T >::m_eps, CGSolver< T >::m_hang, CGSolver< T >::m_homogeneous, CGSolver< T >::m_imax, CGSolver< T >::m_normType, CGSolver< T >::m_op, and CGSolver< T >::m_verbosity.
Referenced by CGSolver< T >::setHomogeneous().
Real CGSolver< T >::m_convergenceMetric |
public member data: solver convergence metric – if negative, use initial residual; if positive, then use m_convergenceMetric
Referenced by CGSolver< T >::setConvergenceMetrics(), and CGSolver< T >::solve().
Real CGSolver< T >::m_eps |
public member data: solver tolerance
Referenced by CGSolver< T >::setConvergenceMetrics(), CGSolver< T >::solve(), and JFNKSolver::solve().
int CGSolver< T >::m_exitStatus |
public member data: set = -1 if solver exited for an unknown reason set = 1 if solver converged to tolerance set = 2 if rho = 0 set = 3 if max number of restarts was reached
Real CGSolver< T >::m_hang |
public member data: minium norm of solution should change per iterations
Referenced by CGSolver< T >::solve(), and JFNKSolver::solve().
bool CGSolver< T >::m_homogeneous |
public member data: whether the solver is restricted to homogeneous boundary conditions
Referenced by CGSolver< T >::define(), CGSolver< T >::setHomogeneous(), and CGSolver< T >::solve().
int CGSolver< T >::m_imax |
public member data: maximum number of iterations
Referenced by CGSolver< T >::solve(), and JFNKSolver::solve().
int CGSolver< T >::m_normType |
public member data: norm to be used when evaluating convergence. 0 is max norm, 1 is L(1), 2 is L(2) and so on.
Referenced by CGSolver< T >::solve(), and JFNKSolver::solve().
int CGSolver< T >::m_numRestarts |
public member data: number of times the algorithm can restart
Referenced by JFNKSolver::solve().
LinearOp<T>* CGSolver< T >::m_op |
public member data: operator to solve.
Referenced by CGSolver< T >::define(), CGSolver< T >::solve(), and CGSolver< T >::~CGSolver().
Real CGSolver< T >::m_reps |
public member data: relative solver tolerance
Real CGSolver< T >::m_small |
public member data: what the algorithm should consider "close to zero"
int CGSolver< T >::m_verbosity |
public member data: how much screen out put the user wants. set = 0 for no output.
Referenced by CGSolver< T >::solve(), and JFNKSolver::solve().