42#ifndef STOKHOS_APPROX_SCHUR_COMPLEMENT_PRECONDITIONER_HPP
43#define STOKHOS_APPROX_SCHUR_COMPLEMENT_PRECONDITIONER_HPP
45#include "Teuchos_RCP.hpp"
48#include "EpetraExt_MultiComm.h"
53#include "Teuchos_ParameterList.hpp"
54#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);
111 virtual double NormInf()
const;
114 virtual const char*
Label ()
const;
150 const EpetraExt::BlockMultiVector& Input,
151 EpetraExt::BlockMultiVector& Result)
const;
154 int row_begin,
int row_end,
155 const EpetraExt::BlockMultiVector& Input,
156 EpetraExt::BlockMultiVector& Result)
const;
172 Teuchos::RCP<const EpetraExt::MultiComm>
sg_comm;
175 Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> >
sg_basis;
196 Teuchos::RCP<Stokhos::SGOperator>
sg_op;
199 Teuchos::RCP<Stokhos::EpetraOperatorOrthogPoly >
sg_poly;
205 Teuchos::RCP<const Cijk_type >
Cijk;
230 mutable Teuchos::RCP<Epetra_MultiVector>
tmp;
233 mutable Teuchos::RCP<EpetraExt::BlockMultiVector>
rhs_block;
235 mutable Teuchos::Array<double*>
j_ptr;
A stochastic preconditioner based on applying the approximate Schur complement preconditioner as defi...
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 ...
ApproxSchurComplementPreconditioner & operator=(const ApproxSchurComplementPreconditioner &)
Private to prohibit copying.
bool only_use_linear
Limit Gauss-Seidel loop to linear terms.
Teuchos::Array< Teuchos::RCP< Cijk_type > > lower_block_Cijk
virtual bool UseTranspose() const
Returns the current UseTranspose setting.
void multiply_block(const Teuchos::RCP< const Stokhos::Sparse3Tensor< int, double > > &cijk, double alpha, const EpetraExt::BlockMultiVector &Input, EpetraExt::BlockMultiVector &Result) const
Teuchos::Array< int > mj_indices
int P
Total polynomial order.
virtual const char * Label() const
Returns a character string describing the operator.
Teuchos::Array< double * > j_ptr
Teuchos::RCP< const Stokhos::EpetraSparse3Tensor > epetraCijk
Stores Epetra Cijk tensor.
Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > sg_poly
Pointer to the PCE expansion of Jacobian.
virtual ~ApproxSchurComplementPreconditioner()
Destructor.
Teuchos::Array< int > block_indices
Starting block indices.
bool scale_op
Flag indicating whether operator be scaled with <\psi_i^2>
Teuchos::RCP< const Stokhos::OrthogPolyBasis< int, double > > sg_basis
Stochastic Galerking basis.
std::string label
Label for 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.
virtual void setupPreconditioner(const Teuchos::RCP< Stokhos::SGOperator > &sg_op, const Epetra_Vector &x)
Setup preconditioner.
Teuchos::RCP< const Cijk_type > Cijk
Pointer to triple product.
virtual double NormInf() const
Returns an approximate infinity norm of the operator matrix.
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< Stokhos::AbstractPreconditionerFactory > prec_factory
Stores factory for building mean preconditioner.
bool symmetric
Use symmetric Gauss-Seidel.
Teuchos::RCP< Epetra_MultiVector > tmp
Temporary vector for storing matrix-vector products.
Teuchos::RCP< EpetraExt::BlockMultiVector > rhs_block
Temporary vector for storing rhs in Gauss-Seidel loop.
Teuchos::RCP< Stokhos::SGOperator > sg_op
Pointer to the SG operator.
int max_num_mat_vec
Maximum number of matvecs in Apply.
Teuchos::RCP< const Epetra_Map > sg_map
Stores SG map.
void divide_diagonal_block(int row_begin, int row_end, const EpetraExt::BlockMultiVector &Input, EpetraExt::BlockMultiVector &Result) const
ApproxSchurComplementPreconditioner(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.
ApproxSchurComplementPreconditioner(const ApproxSchurComplementPreconditioner &)
Private to prohibit copying.
virtual bool HasNormInf() const
Returns true if the this object can provide an approximate Inf-norm, false otherwise.
Teuchos::Array< Teuchos::RCP< Cijk_type > > upper_block_Cijk
Triple product tensor for each sub-block.
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.
bool useTranspose
Flag indicating whether transpose was selected.
Teuchos::RCP< Epetra_Operator > mean_prec
Stores mean preconditioner.
Teuchos::RCP< const EpetraExt::MultiComm > sg_comm
Stores SG parallel communicator.
virtual int SetUseTranspose(bool UseTranspose)
Set to true if the transpose of the operator is requested.
Stokhos::Sparse3Tensor< int, double > Cijk_type
Short-hand for Cijk.
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.