42#ifndef STOKHOS_APPROX_GAUSS_SEIDEL_PRECONDITIONER_HPP
43#define STOKHOS_APPROX_GAUSS_SEIDEL_PRECONDITIONER_HPP
45#include "Teuchos_RCP.hpp"
48#include "EpetraExt_MultiComm.h"
53#include "Teuchos_ParameterList.hpp"
54#include "EpetraExt_BlockMultiVector.h"
68 const Teuchos::RCP<const EpetraExt::MultiComm>&
sg_comm,
70 const Teuchos::RCP<const Stokhos::EpetraSparse3Tensor>&
epetraCijk,
71 const Teuchos::RCP<const Epetra_Map>&
base_map,
72 const Teuchos::RCP<const Epetra_Map>&
sg_map,
73 const Teuchos::RCP<Stokhos::AbstractPreconditionerFactory>&
prec_factory,
74 const Teuchos::RCP<Teuchos::ParameterList>& params);
110 virtual double NormInf()
const;
113 virtual const char*
Label ()
const;
158 Teuchos::RCP<const EpetraExt::MultiComm>
sg_comm;
161 Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> >
sg_basis;
182 Teuchos::RCP<Stokhos::SGOperator>
sg_op;
185 Teuchos::RCP<Stokhos::EpetraOperatorOrthogPoly >
sg_poly;
191 Teuchos::RCP<const Cijk_type >
Cijk;
206 mutable Teuchos::RCP<EpetraExt::BlockMultiVector>
rhs_block;
A stochastic preconditioner based on applying one iteration of approximate Gauss-Seidel.
virtual bool UseTranspose() const
Returns the current UseTranspose setting.
Teuchos::RCP< Epetra_Operator > mean_prec
Stores mean preconditioner.
Teuchos::RCP< const Cijk_type > Cijk
Pointer to triple product.
virtual bool HasNormInf() const
Returns true if the this object can provide an approximate Inf-norm, false otherwise.
Teuchos::RCP< const Stokhos::EpetraSparse3Tensor > epetraCijk
Stores Epetra Cijk tensor.
Teuchos::RCP< Epetra_MultiVector > mat_vec_tmp
Temporary vector for storing matrix-vector products.
virtual const Epetra_Map & OperatorDomainMap() const
Returns the Epetra_Map object associated with the domain of this matrix operator.
bool only_use_linear
Limit Gauss-Seidel loop to linear terms.
virtual const Epetra_Comm & Comm() const
Returns a reference to the Epetra_Comm communicator associated with this operator.
virtual int SetUseTranspose(bool UseTranspose)
Set to true if the transpose of the operator is requested.
bool useTranspose
Flag indicating whether transpose was selected.
virtual void setupPreconditioner(const Teuchos::RCP< Stokhos::SGOperator > &sg_op, const Epetra_Vector &x)
Setup preconditioner.
Teuchos::RCP< const Stokhos::OrthogPolyBasis< int, double > > sg_basis
Stochastic Galerking basis.
Teuchos::RCP< Stokhos::SGOperator > sg_op
Pointer to the SG operator.
virtual ~ApproxGaussSeidelPreconditioner()
Destructor.
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 ...
ApproxGaussSeidelPreconditioner & operator=(const ApproxGaussSeidelPreconditioner &)
Private to prohibit copying.
virtual double NormInf() const
Returns an approximate infinity norm of the operator matrix.
ApproxGaussSeidelPreconditioner(const ApproxGaussSeidelPreconditioner &)
Private to prohibit copying.
bool symmetric
Use symmetric Gauss-Seidel.
Teuchos::RCP< Stokhos::AbstractPreconditionerFactory > prec_factory
Stores factory for building mean preconditioner.
std::string label
Label for operator.
Teuchos::RCP< const Epetra_Map > sg_map
Stores SG map.
Teuchos::RCP< const EpetraExt::MultiComm > sg_comm
Stores SG parallel communicator.
Teuchos::RCP< EpetraExt::BlockMultiVector > rhs_block
Temporary vector for storing rhs in Gauss-Seidel loop.
Stokhos::Sparse3Tensor< int, double > Cijk_type
Short-hand for Cijk.
virtual const char * Label() const
Returns a character string describing the operator.
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 ...
bool scale_op
Flag indicating whether operator be scaled with <\psi_i^2>
Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > sg_poly
Pointer to the PCE expansion of Jacobian.
Teuchos::RCP< const Epetra_Map > base_map
Stores base map.
virtual const Epetra_Map & OperatorRangeMap() const
Returns the Epetra_Map object associated with the range of this matrix operator.
Abstract base class for multivariate orthogonal polynomials.
An abstract class to represent a generic stochastic Galerkin preconditioner as an Epetra_Operator.
Data structure storing a sparse 3-tensor C(i,j,k) in a a compressed format.
Top-level namespace for Stokhos classes and functions.