dune-istl
2.4.1
|
Generalized preconditioned conjugate gradient solver. More...
#include <dune/istl/solvers.hh>
Public Types | |
typedef X | domain_type |
The domain type of the operator to be inverted. More... | |
typedef X | range_type |
The range type of the operator to be inverted. More... | |
typedef X::field_type | field_type |
The field type of the operator to be inverted. More... | |
typedef FieldTraits < field_type >::real_type | real_type |
The real type of the field type (is the same if using real numbers, but differs for std::complex) More... | |
Public Member Functions | |
template<class L , class P > | |
GeneralizedPCGSolver (L &op, P &prec, real_type reduction, int maxit, int verbose, int restart=10) | |
Set up nonlinear preconditioned conjugate gradient solver. More... | |
template<class L , class P , class S > | |
GeneralizedPCGSolver (L &op, S &sp, P &prec, real_type reduction, int maxit, int verbose, int restart=10) | |
Set up nonlinear preconditioned conjugate gradient solver. More... | |
virtual void | apply (X &x, X &b, InverseOperatorResult &res) |
Apply inverse operator. More... | |
virtual void | apply (X &x, X &b, double reduction, InverseOperatorResult &res) |
Apply inverse operator with given reduction factor. More... | |
Protected Types | |
enum | |
Protected Member Functions | |
void | printHeader (std::ostream &s) const |
helper function for printing header of solver output More... | |
void | printOutput (std::ostream &s, const CountType &iter, const DataType &norm, const DataType &norm_old) const |
helper function for printing solver output More... | |
void | printOutput (std::ostream &s, const CountType &iter, const DataType &norm) const |
helper function for printing solver output More... | |
Generalized preconditioned conjugate gradient solver.
A preconditioned conjugate gradient that allows the preconditioner to change between iterations.
One example for such preconditioner is AMG when used without a direct coarse solver. In this case the number of iterations performed on the coarsest level might change between applications.
In contrast to CGSolver the search directions are stored and the orthogonalization is done explicitly.
typedef X Dune::GeneralizedPCGSolver< X >::domain_type |
The domain type of the operator to be inverted.
typedef X::field_type Dune::GeneralizedPCGSolver< X >::field_type |
The field type of the operator to be inverted.
typedef X Dune::GeneralizedPCGSolver< X >::range_type |
The range type of the operator to be inverted.
typedef FieldTraits<field_type>::real_type Dune::GeneralizedPCGSolver< X >::real_type |
The real type of the field type (is the same if using real numbers, but differs for std::complex)
|
protectedinherited |
|
inline |
Set up nonlinear preconditioned conjugate gradient solver.
restart | When to restart the construction of the Krylov search space. |
References Dune::SolverCategory::sequential.
|
inline |
Set up nonlinear preconditioned conjugate gradient solver.
restart | When to restart the construction of the Krylov search space. |
|
inlinevirtual |
Apply inverse operator.
Apply inverse operator,.
x | The left hand side to store the result in. |
b | The right hand side |
res | Object to store the statistics about applying the operator. |
Implements Dune::InverseOperator< X, X >.
References Dune::Preconditioner< X, Y >::apply(), Dune::LinearOperator< X, Y >::apply(), Dune::LinearOperator< X, Y >::applyscaleadd(), Dune::InverseOperatorResult::clear(), Dune::InverseOperatorResult::conv_rate, Dune::InverseOperatorResult::converged, Dune::InverseOperatorResult::elapsed, Dune::InverseOperatorResult::iterations, Dune::Preconditioner< X, Y >::post(), Dune::Preconditioner< X, Y >::pre(), Dune::InverseOperator< X, X >::printHeader(), Dune::InverseOperator< X, X >::printOutput(), and Dune::InverseOperatorResult::reduction.
|
inlinevirtual |
Apply inverse operator with given reduction factor.
apply inverse operator, with given convergence criteria.
x | The left hand side to store the result in. |
b | The right hand side |
reduction | The minimum defect reduction to achieve. |
res | Object to store the statistics about applying the operator. |
Implements Dune::InverseOperator< X, X >.
|
inlineprotectedinherited |
helper function for printing header of solver output
References Dune::InverseOperator< X, Y >::iterationSpacing, and Dune::InverseOperator< X, Y >::normSpacing.
Referenced by Dune::LoopSolver< X >::apply(), Dune::GradientSolver< X >::apply(), Dune::CGSolver< X >::apply(), Dune::BiCGSTABSolver< X >::apply(), Dune::MINRESSolver< X >::apply(), and Dune::GeneralizedPCGSolver< X >::apply().
|
inlineprotectedinherited |
helper function for printing solver output
References Dune::InverseOperator< X, Y >::iterationSpacing, and Dune::InverseOperator< X, Y >::normSpacing.
Referenced by Dune::LoopSolver< X >::apply(), Dune::GradientSolver< X >::apply(), Dune::CGSolver< X >::apply(), Dune::BiCGSTABSolver< X >::apply(), Dune::MINRESSolver< X >::apply(), and Dune::GeneralizedPCGSolver< X >::apply().
|
inlineprotectedinherited |
helper function for printing solver output
References Dune::InverseOperator< X, Y >::iterationSpacing, and Dune::InverseOperator< X, Y >::normSpacing.