template<class ScalarType, class MV, class OP, const bool supportsScalarType = Belos::Details::LapackSupportsScalar<ScalarType>::value && ! Teuchos::ScalarTraits<ScalarType>::isComplex>
class Belos::RCGSolMgr< ScalarType, MV, OP, supportsScalarType >
Implementation of the RCG (Recycling Conjugate Gradient) iterative linear solver.
Author
Michael Parks and Heidi Thornquist
Summary
This class implements the GCRODR (Recycling GMRES) iterative linear solver. This solver is suited for solving sequences of related linear systems , where each matrix is symmetric positive definite.
ScalarType must be real
This RCG implementation currently only supports real-valued (not complex-valued) ScalarType types. You may check whether ScalarType is complex using the following code:
if (Teuchos::ScalarTraits<ScalarType>::isComplex) {
// ScalarType is complex valued.
} else {
// ScalarType is real valued.
}
This is not a limitation of the RCG method itself, just of the current implementation. If there is sufficient interest, we can remedy this deficiency. For now, if you attempt to invoke the constructor when ScalarType is complex, the constructor will throw an exception. This is why this class inherits from Details::RealSolverManager. RCGSolMgr can still compile if ScalarType is complex, but you will not be able to construct a RCGSolMgr instance in that case, due to the aforementioned run-time error that the constructor raises. We do this so that the class will still compile, whether ScalarType is real or complex. This helps make SolverFactory valid to compile, whether ScalarType is real or complex.