Hybrid block GMRES iterative linear solver. More...

#include <BelosGmresPolySolMgr.hpp>

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

Public Member Functions
Constructors/Destructor
	GmresPolySolMgr ()
	Empty constructor for GmresPolySolMgr. 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().
	GmresPolySolMgr (const Teuchos::RCP< LinearProblem< ScalarType, MV, OP > > &problem, const Teuchos::RCP< Teuchos::ParameterList > &pl)
	Basic constructor for GmresPolySolMgr.
virtual	~GmresPolySolMgr ()
	Destructor.
Accessor methods
const LinearProblem < ScalarType, MV, OP > &	getProblem () const
	Get current linear problem being solved for in this object.
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.
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)
	Reset the solver.
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 hybrid block GMRES solver manager.

Detailed Description

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

Hybrid block GMRES iterative linear solver.

Author:: Heidi Thornquist

Examples:: BlockGmres/BlockGmresPolyEpetraExFile.cpp.

Definition at line 156 of file BelosGmresPolySolMgr.hpp.

Constructor & Destructor Documentation

template<class ScalarType , class MV , class OP >

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

Empty constructor for GmresPolySolMgr. 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 434 of file BelosGmresPolySolMgr.hpp.

template<class ScalarType , class MV , class OP >

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

Basic constructor for GmresPolySolMgr.

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" - a int specifying the block size to be used by the underlying block GMRES solver. Default: 1
"Num Blocks" - a int specifying the number of blocks allocated for the Krylov basis. Default: 300
"Maximum Iterations" - a int specifying the maximum number of iterations the underlying solver is allowed to perform. Default: 1000
"Maximum Restarts" - a int specifying the maximum number of restarts the underlying solver is allowed to perform. Default: 20
"Orthogonalization" - a std::string specifying the desired orthogonalization: DGKS, ICGS, and IMGS. Default: "DGKS"
"Verbosity" - a sum of MsgType specifying the verbosity. Default: Belos::Errors
"Output Style" - a OutputType specifying the style of output. Default: Belos::General
"Convergence Tolerance" - a MagnitudeType specifying the level that residual norms must reach to decide convergence. Default: 1e-8

Definition at line 465 of file BelosGmresPolySolMgr.hpp.

template<class ScalarType , class MV , class OP >

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

Destructor.

Definition at line 194 of file BelosGmresPolySolMgr.hpp.

Member Function Documentation

template<class ScalarType , class MV , class OP >

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

Get current linear problem being solved for in this object.

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

Definition at line 202 of file BelosGmresPolySolMgr.hpp.

template<class ScalarType , class MV , class OP >

Teuchos::RCP< const Teuchos::ParameterList > Belos::GmresPolySolMgr< 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 510 of file BelosGmresPolySolMgr.hpp.

template<class ScalarType , class MV , class OP >

Teuchos::RCP<const Teuchos::ParameterList> Belos::GmresPolySolMgr< 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 212 of file BelosGmresPolySolMgr.hpp.

template<class ScalarType , class MV , class OP >

Teuchos::Array<Teuchos::RCP<Teuchos::Time> > Belos::GmresPolySolMgr< 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 219 of file BelosGmresPolySolMgr.hpp.

template<class ScalarType , class MV , class OP >

int Belos::GmresPolySolMgr< 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 224 of file BelosGmresPolySolMgr.hpp.

template<class ScalarType , class MV , class OP >

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

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

Note:: This flag will be reset the next time solve() is called.

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

Definition at line 231 of file BelosGmresPolySolMgr.hpp.

template<class ScalarType , class MV , class OP >

void Belos::GmresPolySolMgr< 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 239 of file BelosGmresPolySolMgr.hpp.

template<class ScalarType , class MV , class OP >

void Belos::GmresPolySolMgr< 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 571 of file BelosGmresPolySolMgr.hpp.

template<class ScalarType , class MV , class OP >

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

Reset the solver.

Parameters:

type	[in] How to reset the solver.

If type includes Belos::Problem, then reset the solver's state. This clears out the stored coefficients, so that the next call to solve() actually computes a full block GMRES solve, instead of just reusing the coefficients from the first solve.

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

Definition at line 256 of file BelosGmresPolySolMgr.hpp.

template<class ScalarType , class MV , class OP >

ReturnType Belos::GmresPolySolMgr< 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 BlockGmresIter::iterate(), which will return either because a specially constructed status test evaluates to Passed or an std::exception is thrown.

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

the maximum number of restarts 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 1075 of file BelosGmresPolySolMgr.hpp.

template<class ScalarType , class MV , class OP >

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

Method to return description of the hybrid block GMRES solver manager.

Reimplemented from Teuchos::Describable.

Definition at line 1432 of file BelosGmresPolySolMgr.hpp.

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

BelosGmresPolySolMgr.hpp

Public Member Functions

Detailed Description

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

Constructor & Destructor Documentation

Member Function Documentation

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