42#ifndef STOKHOS_KRONECKER_PRODUCT_PRECONDITIONER_HPP
43#define STOKHOS_KRONECKER_PRODUCT_PRECONDITIONER_HPP
45#include "Teuchos_RCP.hpp"
48#include "EpetraExt_MultiComm.h"
53#include "Teuchos_ParameterList.hpp"
67 const Teuchos::RCP<const EpetraExt::MultiComm>&
sg_comm,
69 const Teuchos::RCP<const Stokhos::EpetraSparse3Tensor>&
epetraCijk,
70 const Teuchos::RCP<const Epetra_Map>&
base_map,
71 const Teuchos::RCP<const Epetra_Map>&
sg_map,
73 const Teuchos::RCP<Stokhos::AbstractPreconditionerFactory>&
G_prec_factory,
74 const Teuchos::RCP<Teuchos::ParameterList>&
params);
110 virtual double NormInf()
const;
113 virtual const char*
Label ()
const;
164 Teuchos::RCP<const EpetraExt::MultiComm>
sg_comm;
167 Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> >
sg_basis;
185 Teuchos::RCP<Teuchos::ParameterList>
params;
197 Teuchos::RCP<Stokhos::SGOperator>
sg_op;
200 Teuchos::RCP<Stokhos::EpetraOperatorOrthogPoly >
sg_poly;
206 Teuchos::RCP<const Cijk_type>
Cijk;
209 Teuchos::RCP<Epetra_CrsMatrix>
G;
An Epetra operator representing applying the mean in a block stochastic Galerkin expansion.
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< Epetra_Operator > mean_prec
Stores mean preconditioner.
std::string label
Label for operator.
virtual double NormInf() const
Returns an approximate infinity norm of the operator matrix.
Teuchos::RCP< const EpetraExt::MultiComm > sg_comm
Stores SG parallel communicator.
virtual bool UseTranspose() const
Returns the current UseTranspose setting.
bool useTranspose
Flag indicating whether transpose was selected.
virtual bool HasNormInf() const
Returns true if the this object can provide an approximate Inf-norm, false otherwise.
bool only_use_linear
Limit construction of G to linear terms.
virtual void setupPreconditioner(const Teuchos::RCP< Stokhos::SGOperator > &sg_op, const Epetra_Vector &x)
Setup preconditioner.
Teuchos::RCP< const Epetra_Map > base_map
Stores base map.
virtual ~KroneckerProductPreconditioner()
Destructor.
Teuchos::RCP< Teuchos::ParameterList > params
Preconditioner parameters.
Teuchos::RCP< const Cijk_type > Cijk
Stores triple product tensor.
virtual const Epetra_Comm & Comm() const
Returns a reference to the Epetra_Comm communicator associated with this operator.
bool scale_op
Flag indicating whether operator be scaled with <\psi_i^2>
Teuchos::RCP< Stokhos::AbstractPreconditionerFactory > mean_prec_factory
Stores factory for building mean preconditioner.
Teuchos::RCP< const Stokhos::EpetraSparse3Tensor > epetraCijk
Stores Epetra Cijk tensor.
Teuchos::RCP< Epetra_MultiVector > result_MVT
Teuchos::RCP< Epetra_Operator > G_prec
Stores G preconditioner.
double MatrixTrace(const Epetra_CrsMatrix &A, const Epetra_CrsMatrix &B) const
Compute trace of matrix A'B.
Teuchos::RCP< Epetra_CrsMatrix > G
Pointer to CrsMatrix G.
Stokhos::Sparse3Tensor< int, double > Cijk_type
Short-hand for Cijk.
virtual const Epetra_Map & OperatorDomainMap() const
Returns the Epetra_Map object associated with the domain of this matrix operator.
Teuchos::RCP< Stokhos::SGOperator > sg_op
Pointer to the SG operator.
virtual const Epetra_Map & OperatorRangeMap() const
Returns the Epetra_Map object associated with the range of this matrix operator.
Teuchos::RCP< const Stokhos::OrthogPolyBasis< int, double > > sg_basis
Stochastic Galerking basis.
Teuchos::RCP< Stokhos::AbstractPreconditionerFactory > G_prec_factory
Stores factory for building G preconditioner.
Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > sg_poly
Pointer to the PCE expansion of Jacobian.
virtual int SetUseTranspose(bool UseTranspose)
Set to true if the transpose of the operator is requested.
Teuchos::RCP< const Epetra_Map > sg_map
Stores SG map.
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 ...
KroneckerProductPreconditioner(const KroneckerProductPreconditioner &)
Private to prohibit copying.
virtual const char * Label() const
Returns a character string describing the operator.
KroneckerProductPreconditioner & operator=(const KroneckerProductPreconditioner &)
Private to prohibit copying.
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.