Belos's class for applying the GMRES polynomial operator that is used by the hybrid-GMRES linear solver. More...

#include <BelosGmresPolyOp.hpp>

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

Constructor/Destructor
	GmresPolyOp (const Teuchos::RCP< LinearProblem< ScalarType, MV, OP > > &problem_in, const Teuchos::RCP< Teuchos::ParameterList > &params_in)
	Basic contstructor. More...

	GmresPolyOp (const Teuchos::RCP< LinearProblem< ScalarType, MV, OP > > &problem_in)
	Given no ParameterList, constructor creates no polynomial and only applies the given operator. More...

virtual	~GmresPolyOp ()
	Destructor. More...

Parameter processing method
void	setParameters (const Teuchos::RCP< Teuchos::ParameterList > &params_in)
	Process the passed in parameters. More...

Polynomial creation method
void	generateArnoldiPoly ()
	This routine takes the matrix, preconditioner, and vectors from the linear problem as well as the parameters to generate the Arnoldi polynomial. More...

void	generateGmresPoly ()
	This routine takes the matrix, preconditioner, and vectors from the linear problem as well as the parameters to generate the Gmres polynomial. More...

Operator application method
typedef int	OT

typedef MultiVecTraits < ScalarType, MV >	MVT

typedef Teuchos::ScalarTraits < ScalarType >	SCT

typedef Teuchos::ScalarTraits < ScalarType >::magnitudeType	MagnitudeType

typedef Teuchos::ScalarTraits < MagnitudeType >	MCT

std::string	polyUpdateLabel_

Teuchos::RCP< LinearProblem < ScalarType, MV, OP > >	problem_

Teuchos::RCP < Teuchos::ParameterList >	params_

Teuchos::RCP< const OP >	LP_

Teuchos::RCP< const OP >	RP_

Teuchos::RCP< OutputManager < ScalarType > >	printer_

Teuchos::RCP< std::ostream >	outputStream_ = Teuchos::rcp(outputStream_default_,false)

Teuchos::RCP< MatOrthoManager < ScalarType, MV, OP > >	ortho_

MagnitudeType	polyTol_ = DefaultSolverParameters::polyTol

int	maxDegree_ = maxDegree_default_

int	verbosity_ = verbosity_default_

bool	randomRHS_ = randomRHS_default_

std::string	label_ = label_default_

std::string	polyType_ = polyType_default_

std::string	orthoType_ = orthoType_default_

int	dim_ = 0

bool	damp_ = damp_default_

bool	addRoots_ = addRoots_default_

Teuchos::RCP< MV >	V_

Teuchos::RCP< MV >	wL_

Teuchos::RCP< MV >	wR_

Teuchos::SerialDenseMatrix< OT, ScalarType >	H_

Teuchos::SerialDenseMatrix< OT, ScalarType >	y_

Teuchos::SerialDenseVector< OT, ScalarType >	r0_

bool	autoDeg = false

Teuchos::SerialDenseMatrix< OT, ScalarType >	pCoeff_

Teuchos::SerialDenseMatrix< OT, MagnitudeType >	theta_

static constexpr int	maxDegree_default_ = 25

static constexpr int	verbosity_default_ = Belos::Errors

static constexpr bool	randomRHS_default_ = true

static constexpr const char *	label_default_ = "Belos"

static constexpr const char *	polyType_default_ = "Roots"

static constexpr const char *	orthoType_default_ = "DGKS"

static constexpr std::ostream *	outputStream_default_ = &std::cout

static constexpr bool	damp_default_ = false

static constexpr bool	addRoots_default_ = true

void	ApplyPoly (const MV &x, MV &y) const
	This routine takes the MV `x` and applies the polynomial operator `phi(OP)` to it resulting in the MV `y`, which is returned. More...

void	ApplyArnoldiPoly (const MV &x, MV &y) const

void	ApplyGmresPoly (const MV &x, MV &y) const

void	ApplyRootsPoly (const MV &x, MV &y) const

void	Apply (const MultiVec< ScalarType > &x, MultiVec< ScalarType > &y, ETrans=NOTRANS) const
	This routine casts the MultiVec to GmresPolyMv to retrieve the MV. Then the above apply method is called. More...

int	polyDegree () const

void	SortModLeja (Teuchos::SerialDenseMatrix< OT, MagnitudeType > &thetaN, std::vector< int > &index) const

void	ComputeAddedRoots ()

Detailed Description

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

Belos's class for applying the GMRES polynomial operator that is used by the hybrid-GMRES linear solver.

This operator is used as the interface to the matrix polynomial (phi(A)), solution (X), and right-hand side (B) of the linear system phi(A)X = B. Furthermore, it is also the interface to left/right preconditioning of the linear system.

Author: Heidi Thornquist and Jennifer Loe

Definition at line 198 of file BelosGmresPolyOp.hpp.