The Belos::BlockCGSolMgr provides a powerful and fully-featured solver manager over the CG and BlockCG linear solver. More...

#include <BelosBlockCGSolMgr.hpp>

Inheritance diagram for Belos::BlockCGSolMgr< ScalarType, MV, OP >:

[legend]

List of all members.

Public Member Functions
Constructors/Destructor
	BlockCGSolMgr ()
	Empty constructor for BlockCGSolMgr. This constructor takes no arguments and sets the default values for the solver. The linear problem must be passed in using setProblem() before solve() is called on this object. The solver values can be changed using setParameters().
	BlockCGSolMgr (const Teuchos::RCP< LinearProblem< ScalarType, MV, OP > > &problem, const Teuchos::RCP< Teuchos::ParameterList > &pl)
	Basic constructor for BlockCGSolMgr.
virtual	~BlockCGSolMgr ()
	Destructor.
Accessor methods
const LinearProblem < ScalarType, MV, OP > &	getProblem () const
	Return a reference to the linear problem being solved by this solver manager.
Teuchos::RCP< const Teuchos::ParameterList >	getValidParameters () const
	Get a parameter list containing the valid parameters for this object.
Teuchos::RCP< const Teuchos::ParameterList >	getCurrentParameters () const
	Get a parameter list containing the current parameters for this object.
Teuchos::Array< Teuchos::RCP < Teuchos::Time > >	getTimers () const
	Return the timers for this object.
MagnitudeType	achievedTol () const
	Tolerance achieved by the last `solve()` invocation.
int	getNumIters () const
	Get the iteration count for the most recent call to `solve()`.
bool	isLOADetected () const
	Return whether a loss of accuracy was detected by this solver during the most current solve.
Set methods
void	setProblem (const Teuchos::RCP< LinearProblem< ScalarType, MV, OP > > &problem)
	Set the linear problem that needs to be solved.
void	setParameters (const Teuchos::RCP< Teuchos::ParameterList > &params)
	Set the parameters the solver manager should use to solve the linear problem.
Reset methods
void	reset (const ResetType type)
	Performs a reset of the solver manager specified by the `ResetType`. This informs the solver manager that the solver should prepare for the next call to solve by resetting certain elements of the iterative solver strategy.
Solver application methods
ReturnType	solve ()
	This method performs possibly repeated calls to the underlying linear solver's iterate() routine until the problem has been solved (as decided by the solver manager) or the solver manager decides to quit.
Overridden from Teuchos::Describable
std::string	description () const
	Method to return description of the block CG solver manager.

Detailed Description

template<class ScalarType, class MV, class OP>
class Belos::BlockCGSolMgr< ScalarType, MV, OP >

The Belos::BlockCGSolMgr provides a powerful and fully-featured solver manager over the CG and BlockCG linear solver.

Author:: Heidi Thornquist, Chris Baker, and Teri Barth

Examples:: BlockCG/BlockCGEpetraExFile.cpp, and BlockCG/BlockPrecCGEpetraExFile.cpp.

Definition at line 114 of file BelosBlockCGSolMgr.hpp.

Constructor & Destructor Documentation

template<class ScalarType , class MV , class OP >

Belos::BlockCGSolMgr< ScalarType, MV, OP >::BlockCGSolMgr ( )

Empty constructor for BlockCGSolMgr. This constructor takes no arguments and sets the default values for the solver. The linear problem must be passed in using setProblem() before solve() is called on this object. The solver values can be changed using setParameters().

Definition at line 391 of file BelosBlockCGSolMgr.hpp.

template<class ScalarType , class MV , class OP >

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

Basic constructor for BlockCGSolMgr.

This constructor accepts the LinearProblem to be solved in addition to a parameter list of options for the solver manager. These options include the following:

"Block Size" - an int specifying the block size to be used by the underlying block conjugate-gradient solver. Default: 1
"Adaptive Block Size" - a bool specifying whether the block size can be modified throughout the solve. Default: true
"Maximum Iterations" - an int specifying the maximum number of iterations the underlying solver is allowed to perform. Default: 1000
"Convergence Tolerance" - a MagnitudeType specifying the level that residual norms must reach to decide convergence. Default: 1e-8.
"Orthogonalization" - a std::string specifying the desired orthogonalization: DGKS ,ICGS, and IMGS. Default: "DGKS"
"Orthogonalization Constant" - a MagnitudeType used by DGKS orthogonalization to determine whether another step of classical Gram-Schmidt is necessary. Default: -1 (use DGKS default)
"Verbosity" - a sum of MsgType specifying the verbosity. Default: Belos::Errors
"Output Style" - a OutputType specifying the style of output. Default: Belos::General
"Output Stream" - a reference-counted pointer to the output stream where all solver output is sent. Default: Teuchos::rcp(&std::cout,false)
"Output Frequency" - an int specifying how often convergence information should be outputted. Default: -1 (never)
"Show Maximum Residual Norm Only" - a bool specifying whether that only the maximum relative residual norm is printed if convergence information is printed. Default: false
"Timer Label" - a std::string to use as a prefix for the timer labels. Default: "Belos"

Definition at line 413 of file BelosBlockCGSolMgr.hpp.

template<class ScalarType , class MV , class OP >

virtual Belos::BlockCGSolMgr< ScalarType, MV, OP >::~BlockCGSolMgr ( ) [inline, virtual]

Destructor.

Definition at line 167 of file BelosBlockCGSolMgr.hpp.

Member Function Documentation

template<class ScalarType , class MV , class OP >

const LinearProblem<ScalarType,MV,OP>& Belos::BlockCGSolMgr< ScalarType, MV, OP >::getProblem ( ) const [inline, virtual]

Return a reference to the linear problem being solved by this solver manager.

Implements Belos::SolverManager< ScalarType, MV, OP >.

Definition at line 173 of file BelosBlockCGSolMgr.hpp.

template<class ScalarType , class MV , class OP >

Teuchos::RCP< const Teuchos::ParameterList > Belos::BlockCGSolMgr< ScalarType, MV, OP >::getValidParameters ( ) const [virtual]

Get a parameter list containing the valid parameters for this object.

Implements Belos::SolverManager< ScalarType, MV, OP >.

Definition at line 683 of file BelosBlockCGSolMgr.hpp.

template<class ScalarType , class MV , class OP >

Teuchos::RCP<const Teuchos::ParameterList> Belos::BlockCGSolMgr< ScalarType, MV, OP >::getCurrentParameters ( ) const [inline, virtual]

Get a parameter list containing the current parameters for this object.

Implements Belos::SolverManager< ScalarType, MV, OP >.

Definition at line 183 of file BelosBlockCGSolMgr.hpp.

template<class ScalarType , class MV , class OP >

Teuchos::Array<Teuchos::RCP<Teuchos::Time> > Belos::BlockCGSolMgr< ScalarType, MV, OP >::getTimers ( ) const [inline]

Return the timers for this object.

The timers are ordered as follows:

time spent in solve() routine

Definition at line 190 of file BelosBlockCGSolMgr.hpp.

template<class ScalarType , class MV , class OP >

MagnitudeType Belos::BlockCGSolMgr< ScalarType, MV, OP >::achievedTol ( ) const [inline, virtual]

Tolerance achieved by the last solve() invocation.

This is the maximum over all right-hand sides' achieved convergence tolerances, and is set whether or not the solve actually managed to achieve the desired convergence tolerance.

Reimplemented from Belos::SolverManager< ScalarType, MV, OP >.

Definition at line 199 of file BelosBlockCGSolMgr.hpp.

template<class ScalarType , class MV , class OP >

int Belos::BlockCGSolMgr< ScalarType, MV, OP >::getNumIters ( ) const [inline, virtual]

Get the iteration count for the most recent call to solve().

Implements Belos::SolverManager< ScalarType, MV, OP >.

Definition at line 204 of file BelosBlockCGSolMgr.hpp.

template<class ScalarType , class MV , class OP >

bool Belos::BlockCGSolMgr< ScalarType, MV, OP >::isLOADetected ( ) const [inline, virtual]

Return whether a loss of accuracy was detected by this solver during the most current solve.

Implements Belos::SolverManager< ScalarType, MV, OP >.

Definition at line 210 of file BelosBlockCGSolMgr.hpp.

template<class ScalarType , class MV , class OP >

void Belos::BlockCGSolMgr< ScalarType, MV, OP >::setProblem ( const Teuchos::RCP< LinearProblem< ScalarType, MV, OP > > & problem ) [inline]

Set the linear problem that needs to be solved.

Definition at line 218 of file BelosBlockCGSolMgr.hpp.

template<class ScalarType , class MV , class OP >

void Belos::BlockCGSolMgr< ScalarType, MV, OP >::setParameters ( const Teuchos::RCP< Teuchos::ParameterList > & params )

Set the parameters the solver manager should use to solve the linear problem.

Definition at line 446 of file BelosBlockCGSolMgr.hpp.

template<class ScalarType , class MV , class OP >

void Belos::BlockCGSolMgr< ScalarType, MV, OP >::reset ( const ResetType type ) [inline, virtual]

Performs a reset of the solver manager specified by the ResetType. This informs the solver manager that the solver should prepare for the next call to solve by resetting certain elements of the iterative solver strategy.

Implements Belos::SolverManager< ScalarType, MV, OP >.

Definition at line 231 of file BelosBlockCGSolMgr.hpp.

template<class ScalarType , class MV , class OP >

ReturnType Belos::BlockCGSolMgr< ScalarType, MV, OP >::solve ( ) [virtual]

This method performs possibly repeated calls to the underlying linear solver's iterate() routine until the problem has been solved (as decided by the solver manager) or the solver manager decides to quit.

This method calls BlockCGIter::iterate() or CGIter::iterate(), which will return either because a specially constructed status test evaluates to Passed or an std::exception is thrown.

A return from BlockCGIter::iterate() signifies one of the following scenarios:

the maximum number of iterations has been exceeded. In this scenario, the current solutions to the linear system will be placed in the linear problem and return Unconverged.
global convergence has been met. In this case, the current solutions to the linear system will be placed in the linear problem and the solver manager will return Converged

Returns:

ReturnType specifying:

Converged: the linear problem was solved to the specification required by the solver manager.
Unconverged: the linear problem was not solved to the specification desired by the solver manager.

Implements Belos::SolverManager< ScalarType, MV, OP >.

Definition at line 732 of file BelosBlockCGSolMgr.hpp.

template<class ScalarType , class MV , class OP >

std::string Belos::BlockCGSolMgr< ScalarType, MV, OP >::description ( ) const [virtual]

Method to return description of the block CG solver manager.

Reimplemented from Teuchos::Describable.

Definition at line 1018 of file BelosBlockCGSolMgr.hpp.

The documentation for this class was generated from the following file:

BelosBlockCGSolMgr.hpp

Public Member Functions

Detailed Description

template<class ScalarType, class MV, class OP> class Belos::BlockCGSolMgr< ScalarType, MV, OP >

Constructor & Destructor Documentation

Member Function Documentation

template<class ScalarType, class MV, class OP>
class Belos::BlockCGSolMgr< ScalarType, MV, OP >