42 #ifndef STOKHOS_APPROX_JACOBI_PRECONDITIONER_HPP 43 #define STOKHOS_APPROX_JACOBI_PRECONDITIONER_HPP 45 #include "Teuchos_RCP.hpp" 48 #include "EpetraExt_MultiComm.h" 51 #include "Epetra_Map.h" 53 #include "Teuchos_ParameterList.hpp" 55 #include "EpetraExt_BlockMultiVector.h" 69 const Teuchos::RCP<const EpetraExt::MultiComm>&
sg_comm,
71 const Teuchos::RCP<const Stokhos::EpetraSparse3Tensor>&
epetraCijk,
72 const Teuchos::RCP<const Epetra_Map>&
base_map,
73 const Teuchos::RCP<const Epetra_Map>&
sg_map,
74 const Teuchos::RCP<Stokhos::AbstractPreconditionerFactory>&
prec_factory,
75 const Teuchos::RCP<Teuchos::ParameterList>& params);
86 const Epetra_Vector&
x);
100 virtual int Apply(
const Epetra_MultiVector& Input,
101 Epetra_MultiVector& Result)
const;
108 Epetra_MultiVector& Y)
const;
111 virtual double NormInf()
const;
114 virtual const char*
Label ()
const;
129 virtual const Epetra_Comm &
Comm()
const;
159 Teuchos::RCP<const EpetraExt::MultiComm>
sg_comm;
162 Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> >
sg_basis;
186 Teuchos::RCP<Stokhos::SGOperator>
sg_op;
189 Teuchos::RCP<Stokhos::EpetraOperatorOrthogPoly >
sg_poly;
195 mutable Teuchos::RCP<EpetraExt::BlockMultiVector>
rhs_block;
201 #endif // STOKHOS_APPROX_JACOBI_PRECONDITIONER_HPP Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > sg_poly
Pointer to the PCE expansion of Jacobian.
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.
A stochastic preconditioner based on applying two iterations of approximate Jacobi.
Teuchos::RCP< const Epetra_Map > sg_map
Stores SG map.
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 ...
virtual ~ApproxJacobiPreconditioner()
Destructor.
virtual bool UseTranspose() const
Returns the current UseTranspose setting.
Teuchos::RCP< const Stokhos::EpetraSparse3Tensor > epetraCijk
Stores Epetra Cijk tensor.
Teuchos::RCP< const Epetra_Map > base_map
Stores base map.
Teuchos::RCP< Stokhos::SGOperator > sg_op
Pointer to the SG operator.
Teuchos::RCP< Stokhos::SGOperator > mat_free_op
SG operator to implement SG mat-vec.
int num_iter
Number of Jacobi iterations.
const IndexType const IndexType const IndexType const IndexType const ValueType const ValueType * x
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 ...
Top-level namespace for Stokhos classes and functions.
virtual bool HasNormInf() const
Returns true if the this object can provide an approximate Inf-norm, false otherwise.
virtual int SetUseTranspose(bool UseTranspose)
Set to true if the transpose of the operator is requested.
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.
virtual double NormInf() const
Returns an approximate infinity norm of the operator matrix.
ApproxJacobiPreconditioner & operator=(const ApproxJacobiPreconditioner &)
Private to prohibit copying.
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 > ¶ms)
Constructor.
virtual const Epetra_Map & OperatorRangeMap() const
Returns the Epetra_Map object associated with the range of this matrix operator.
virtual const Epetra_Map & OperatorDomainMap() const
Returns the Epetra_Map object associated with the domain of this matrix operator. ...
virtual const Epetra_Comm & Comm() const
Returns a reference to the Epetra_Comm communicator associated with this operator.
std::string label
Label for operator.
Teuchos::RCP< Epetra_Operator > mean_prec
Stores mean preconditioner.
An abstract class to represent a generic stochastic Galerkin preconditioner as an Epetra_Operator...
bool useTranspose
Flag indicating whether transpose was selected.
Teuchos::RCP< Stokhos::AbstractPreconditionerFactory > prec_factory
Stores factory for building mean preconditioner.