A stochastic preconditioner based on applying two iterations of approximate Jacobi. More...

#include <Stokhos_ApproxJacobiPreconditioner.hpp>

Inheritance diagram for Stokhos::ApproxJacobiPreconditioner:

Collaboration diagram for Stokhos::ApproxJacobiPreconditioner:

Public Member Functions
	ApproxJacobiPreconditioner (const Teuchos::RCP< const EpetraExt::MultiComm > &sg_comm, const Teuchos::RCP< const Stokhos::OrthogPolyBasis< int, double > > &sg_basis, const Teuchos::RCP< const Stokhos::EpetraSparse3Tensor > &epetraCijk, const Teuchos::RCP< const Epetra_Map > &base_map, const Teuchos::RCP< const Epetra_Map > &sg_map, const Teuchos::RCP< Stokhos::AbstractPreconditionerFactory > &prec_factory, const Teuchos::RCP< Teuchos::ParameterList > &params)
	Constructor.
virtual	~ApproxJacobiPreconditioner ()
	Destructor.
Stokhos::SGPreconditioner methods
virtual void	setupPreconditioner (const Teuchos::RCP< Stokhos::SGOperator > &sg_op, const Epetra_Vector &x)
	Setup preconditioner.
Public Member Functions inherited from Stokhos::SGPreconditioner
	SGPreconditioner ()
	Constructor.
virtual	~SGPreconditioner ()
	Destructor.

Epetra_Operator methods
std::string	label
	Label for operator.
Teuchos::RCP< const EpetraExt::MultiComm >	sg_comm
	Stores SG parallel communicator.
Teuchos::RCP< const Stokhos::OrthogPolyBasis< int, double > >	sg_basis
	Stochastic Galerking basis.
Teuchos::RCP< const Stokhos::EpetraSparse3Tensor >	epetraCijk
	Stores Epetra Cijk tensor.
Teuchos::RCP< const Epetra_Map >	base_map
	Stores base map.
Teuchos::RCP< const Epetra_Map >	sg_map
	Stores SG map.
Teuchos::RCP< Stokhos::AbstractPreconditionerFactory >	prec_factory
	Stores factory for building mean preconditioner.
Teuchos::RCP< Epetra_Operator >	mean_prec
	Stores mean preconditioner.
bool	useTranspose
	Flag indicating whether transpose was selected.
int	num_iter
	Number of Jacobi iterations.
Teuchos::RCP< Stokhos::SGOperator >	sg_op
	Pointer to the SG operator.
Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly >	sg_poly
	Pointer to the PCE expansion of Jacobian.
Teuchos::RCP< Stokhos::SGOperator >	mat_free_op
	SG operator to implement SG mat-vec.
Teuchos::RCP< EpetraExt::BlockMultiVector >	rhs_block
	Temporary vector for storing rhs in Gauss-Seidel loop.
virtual int	SetUseTranspose (bool UseTranspose)
	Set to true if the transpose of the operator is requested.
virtual int	Apply (const Epetra_MultiVector &Input, Epetra_MultiVector &Result) const
	Returns the result of a Epetra_Operator applied to a Epetra_MultiVector Input in Result as described above.
virtual int	ApplyInverse (const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
	Returns the result of the inverse of the operator applied to a Epetra_MultiVector Input in Result as described above.
virtual double	NormInf () const
	Returns an approximate infinity norm of the operator matrix.
virtual const char *	Label () const
	Returns a character string describing the operator.
virtual bool	UseTranspose () const
	Returns the current UseTranspose setting.
virtual bool	HasNormInf () const
	Returns true if the this object can provide an approximate Inf-norm, false otherwise.
virtual const Epetra_Comm &	Comm () const
	Returns a reference to the Epetra_Comm communicator associated with this operator.
virtual const Epetra_Map &	OperatorDomainMap () const
	Returns the Epetra_Map object associated with the domain of this matrix operator.
virtual const Epetra_Map &	OperatorRangeMap () const
	Returns the Epetra_Map object associated with the range of this matrix operator.

Detailed Description

A stochastic preconditioner based on applying two iterations of approximate Jacobi.

Constructor & Destructor Documentation

◆ ApproxJacobiPreconditioner()

Stokhos::ApproxJacobiPreconditioner::ApproxJacobiPreconditioner	(	const Teuchos::RCP< const EpetraExt::MultiComm > &	sg_comm,
		const Teuchos::RCP< const Stokhos::OrthogPolyBasis< int, double > > &	sg_basis,
		const Teuchos::RCP< const Stokhos::EpetraSparse3Tensor > &	epetraCijk,
		const Teuchos::RCP< const Epetra_Map > &	base_map,
		const Teuchos::RCP< const Epetra_Map > &	sg_map,
		const Teuchos::RCP< Stokhos::AbstractPreconditionerFactory > &	prec_factory,
		const Teuchos::RCP< Teuchos::ParameterList > &	params )

Constructor.

References base_map, Stokhos::SGOperatorFactory::build(), epetraCijk, Stokhos::EpetraSparse3Tensor::EpetraSparse3Tensor(), label, mat_free_op, mean_prec, num_iter, prec_factory, rhs_block, sg_basis, sg_comm, sg_map, sg_op, sg_poly, and useTranspose.

Member Function Documentation

◆ Apply()

int Stokhos::ApproxJacobiPreconditioner::Apply	(	const Epetra_MultiVector &	Input,
		Epetra_MultiVector &	Result ) const

virtual

Returns the result of a Epetra_Operator applied to a Epetra_MultiVector Input in Result as described above.

References sg_op.

◆ ApplyInverse()

int Stokhos::ApproxJacobiPreconditioner::ApplyInverse	(	const Epetra_MultiVector &	X,
		Epetra_MultiVector &	Y ) const

virtual

Returns the result of the inverse of the operator applied to a Epetra_MultiVector Input in Result as described above.

References base_map, epetraCijk, mat_free_op, mean_prec, num_iter, rhs_block, and sg_map.

◆ Comm()

const Epetra_Comm & Stokhos::ApproxJacobiPreconditioner::Comm ( ) const

virtual

Returns a reference to the Epetra_Comm communicator associated with this operator.

References sg_comm.

◆ HasNormInf()

bool Stokhos::ApproxJacobiPreconditioner::HasNormInf ( ) const

virtual

Returns true if the this object can provide an approximate Inf-norm, false otherwise.

References sg_op.

◆ Label()

const char * Stokhos::ApproxJacobiPreconditioner::Label ( ) const

virtual

Returns a character string describing the operator.

References label.

◆ NormInf()

double Stokhos::ApproxJacobiPreconditioner::NormInf ( ) const

virtual

Returns an approximate infinity norm of the operator matrix.

References sg_op.

◆ OperatorDomainMap()

const Epetra_Map & Stokhos::ApproxJacobiPreconditioner::OperatorDomainMap ( ) const

virtual

Returns the Epetra_Map object associated with the domain of this matrix operator.

References sg_map.

◆ OperatorRangeMap()

const Epetra_Map & Stokhos::ApproxJacobiPreconditioner::OperatorRangeMap ( ) const

virtual

Returns the Epetra_Map object associated with the range of this matrix operator.

References sg_map.

◆ setupPreconditioner()

void Stokhos::ApproxJacobiPreconditioner::setupPreconditioner	(	const Teuchos::RCP< Stokhos::SGOperator > &	sg_op,
		const Epetra_Vector &	x )

virtual

Setup preconditioner.

Implements Stokhos::SGPreconditioner.

References label, mat_free_op, mean_prec, prec_factory, sg_op, and sg_poly.

◆ SetUseTranspose()

int Stokhos::ApproxJacobiPreconditioner::SetUseTranspose ( bool UseTranspose )

virtual

Set to true if the transpose of the operator is requested.

References mat_free_op, mean_prec, sg_op, UseTranspose(), and useTranspose.

◆ UseTranspose()

bool Stokhos::ApproxJacobiPreconditioner::UseTranspose ( ) const

virtual

Returns the current UseTranspose setting.

References useTranspose.

Referenced by SetUseTranspose().

Member Data Documentation

◆ base_map

Teuchos::RCP<const Epetra_Map> Stokhos::ApproxJacobiPreconditioner::base_map

protected

Stores base map.

Referenced by ApplyInverse(), and ApproxJacobiPreconditioner().

◆ epetraCijk

Teuchos::RCP<const Stokhos::EpetraSparse3Tensor> Stokhos::ApproxJacobiPreconditioner::epetraCijk

protected

Stores Epetra Cijk tensor.

Referenced by ApplyInverse(), and ApproxJacobiPreconditioner().

◆ label

std::string Stokhos::ApproxJacobiPreconditioner::label

protected

Label for operator.

Referenced by ApproxJacobiPreconditioner(), Label(), and setupPreconditioner().

◆ mat_free_op

Teuchos::RCP<Stokhos::SGOperator> Stokhos::ApproxJacobiPreconditioner::mat_free_op

protected

SG operator to implement SG mat-vec.

Referenced by ApplyInverse(), ApproxJacobiPreconditioner(), setupPreconditioner(), and SetUseTranspose().

◆ mean_prec

Teuchos::RCP<Epetra_Operator> Stokhos::ApproxJacobiPreconditioner::mean_prec

protected

Stores mean preconditioner.

Referenced by ApplyInverse(), ApproxJacobiPreconditioner(), setupPreconditioner(), and SetUseTranspose().

◆ num_iter

int Stokhos::ApproxJacobiPreconditioner::num_iter

protected

Number of Jacobi iterations.

Referenced by ApplyInverse(), and ApproxJacobiPreconditioner().

◆ prec_factory

Teuchos::RCP<Stokhos::AbstractPreconditionerFactory> Stokhos::ApproxJacobiPreconditioner::prec_factory

protected

Stores factory for building mean preconditioner.

Referenced by ApproxJacobiPreconditioner(), and setupPreconditioner().

◆ rhs_block

Teuchos::RCP<EpetraExt::BlockMultiVector> Stokhos::ApproxJacobiPreconditioner::rhs_block

mutableprotected

Temporary vector for storing rhs in Gauss-Seidel loop.

Referenced by ApplyInverse(), and ApproxJacobiPreconditioner().

◆ sg_basis

Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> > Stokhos::ApproxJacobiPreconditioner::sg_basis

protected

Stochastic Galerking basis.

Referenced by ApproxJacobiPreconditioner().

◆ sg_comm

Teuchos::RCP<const EpetraExt::MultiComm> Stokhos::ApproxJacobiPreconditioner::sg_comm

protected

Stores SG parallel communicator.

Referenced by ApproxJacobiPreconditioner(), and Comm().

◆ sg_map

Teuchos::RCP<const Epetra_Map> Stokhos::ApproxJacobiPreconditioner::sg_map

protected

Stores SG map.

Referenced by ApplyInverse(), ApproxJacobiPreconditioner(), OperatorDomainMap(), and OperatorRangeMap().

◆ sg_op

Teuchos::RCP<Stokhos::SGOperator> Stokhos::ApproxJacobiPreconditioner::sg_op

protected

Pointer to the SG operator.

Referenced by Apply(), ApproxJacobiPreconditioner(), HasNormInf(), NormInf(), setupPreconditioner(), and SetUseTranspose().

◆ sg_poly

Teuchos::RCP<Stokhos::EpetraOperatorOrthogPoly > Stokhos::ApproxJacobiPreconditioner::sg_poly

protected

Pointer to the PCE expansion of Jacobian.

Referenced by ApproxJacobiPreconditioner(), and setupPreconditioner().

◆ useTranspose

bool Stokhos::ApproxJacobiPreconditioner::useTranspose

protected

Flag indicating whether transpose was selected.

Referenced by ApproxJacobiPreconditioner(), SetUseTranspose(), and UseTranspose().

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

Stokhos_ApproxJacobiPreconditioner.hpp
Stokhos_ApproxJacobiPreconditioner.cpp

Public Member Functions

Epetra_Operator methods

Detailed Description

Constructor & Destructor Documentation

◆ ApproxJacobiPreconditioner()

Member Function Documentation

◆ Apply()

◆ ApplyInverse()

◆ Comm()

◆ HasNormInf()

◆ Label()

◆ NormInf()

◆ OperatorDomainMap()

◆ OperatorRangeMap()

◆ setupPreconditioner()

◆ SetUseTranspose()

◆ UseTranspose()

Member Data Documentation

◆ base_map

◆ epetraCijk

◆ label

◆ mat_free_op

◆ mean_prec

◆ num_iter

◆ prec_factory

◆ rhs_block

◆ sg_basis

◆ sg_comm

◆ sg_map

◆ sg_op

◆ sg_poly

◆ useTranspose