Implementation of the RCG (Recycling Conjugate Gradient) iterative linear solver. More...

#include <BelosRCGSolMgr.hpp>

Inheritance diagram for Belos::RCGSolMgr< ScalarType, MV, OP, supportsScalarType >:

Public Member Functions
	RCGSolMgr ()

	RCGSolMgr (const Teuchos::RCP< LinearProblem< ScalarType, MV, OP > > &problem, const Teuchos::RCP< Teuchos::ParameterList > &pl)

virtual	~RCGSolMgr ()

Teuchos::RCP< SolverManager < ScalarType, MV, OP > >	clone () const override
	clone for Inverted Injection (DII) More...

Private Types
typedef Details::SolverManagerRequiresRealLapack < ScalarType, MV, OP, scalarTypeIsSupported >	base_type

Static Private Attributes
static const bool	scalarTypeIsSupported

Detailed Description

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 $A_i x_i = b_i$ , where each matrix $A_i$ 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.

Definition at line 115 of file BelosRCGSolMgr.hpp.