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" 51 #include "Epetra_Map.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);
85 const Epetra_Vector&
x);
99 virtual int Apply(
const Epetra_MultiVector& Input,
100 Epetra_MultiVector& Result)
const;
107 Epetra_MultiVector& Y)
const;
110 virtual double NormInf()
const;
113 virtual const char*
Label ()
const;
128 virtual const Epetra_Comm &
Comm()
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;
212 #endif // STOKHOS_APPROX_GAUSS_SEIDEL_PRECONDITIONER_HPP Teuchos::RCP< const Epetra_Map > sg_map
Stores SG map.
virtual bool HasNormInf() const
Returns true if the this object can provide an approximate Inf-norm, false otherwise.
A stochastic preconditioner based on applying one iteration of approximate Gauss-Seidel.
virtual double NormInf() const
Returns an approximate infinity norm of the operator matrix.
virtual bool UseTranspose() const
Returns the current UseTranspose setting.
Teuchos::RCP< EpetraExt::BlockMultiVector > rhs_block
Temporary vector for storing rhs in Gauss-Seidel loop.
ApproxGaussSeidelPreconditioner & operator=(const ApproxGaussSeidelPreconditioner &)
Private to prohibit copying.
Teuchos::RCP< Epetra_Operator > mean_prec
Stores mean preconditioner.
Teuchos::RCP< const EpetraExt::MultiComm > sg_comm
Stores SG parallel communicator.
ApproxGaussSeidelPreconditioner(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 > ¶ms)
Constructor.
virtual ~ApproxGaussSeidelPreconditioner()
Destructor.
bool useTranspose
Flag indicating whether transpose was selected.
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 ...
virtual const char * Label() const
Returns a character string describing the operator.
Stokhos::Sparse3Tensor< int, double > Cijk_type
Short-hand for Cijk.
const IndexType const IndexType const IndexType const IndexType const ValueType const ValueType * x
bool symmetric
Use symmetric Gauss-Seidel.
Teuchos::RCP< Stokhos::AbstractPreconditionerFactory > prec_factory
Stores factory for building mean preconditioner.
bool scale_op
Flag indicating whether operator be scaled with <^2>
Teuchos::RCP< const Cijk_type > Cijk
Pointer to triple product.
Top-level namespace for Stokhos classes and functions.
bool only_use_linear
Limit Gauss-Seidel loop to linear terms.
Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > sg_poly
Pointer to the PCE expansion of Jacobian.
std::string label
Label for operator.
Teuchos::RCP< const Stokhos::OrthogPolyBasis< int, double > > sg_basis
Stochastic Galerking basis.
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.
virtual int SetUseTranspose(bool UseTranspose)
Set to true if the transpose of the operator is requested.
virtual const Epetra_Comm & Comm() const
Returns a reference to the Epetra_Comm communicator associated with this operator.
virtual void setupPreconditioner(const Teuchos::RCP< Stokhos::SGOperator > &sg_op, const Epetra_Vector &x)
Setup preconditioner.
An abstract class to represent a generic stochastic Galerkin preconditioner as an Epetra_Operator...
Teuchos::RCP< Stokhos::SGOperator > sg_op
Pointer to the SG operator.
virtual const Epetra_Map & OperatorDomainMap() const
Returns the Epetra_Map object associated with the domain of this matrix operator. ...
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 ...
Teuchos::RCP< const Stokhos::EpetraSparse3Tensor > epetraCijk
Stores Epetra Cijk tensor.
Teuchos::RCP< Epetra_MultiVector > mat_vec_tmp
Temporary vector for storing matrix-vector products.