44#ifndef STOKHOS_KL_REDUCED_MATRIX_FREE_OPERATOR_HPP
45#define STOKHOS_KL_REDUCED_MATRIX_FREE_OPERATOR_HPP
48#include "EpetraExt_MultiComm.h"
52#include "Teuchos_ParameterList.hpp"
53#include "Teuchos_Array.hpp"
73 const Teuchos::RCP<const EpetraExt::MultiComm>&
sg_comm,
75 const Teuchos::RCP<const Stokhos::EpetraSparse3Tensor>&
epetraCijk,
80 const Teuchos::RCP<Teuchos::ParameterList>&
params);
90 const Teuchos::RCP<Stokhos::EpetraOperatorOrthogPoly >& poly);
93 virtual Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly >
97 virtual Teuchos::RCP<const Stokhos::EpetraOperatorOrthogPoly >
123 virtual double NormInf()
const;
126 virtual const char*
Label ()
const;
176 Teuchos::RCP<const EpetraExt::MultiComm>
sg_comm;
179 Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> >
sg_basis;
200 Teuchos::RCP<const Cijk_type>
Cijk;
203 Teuchos::RCP<Stokhos::EpetraOperatorOrthogPoly >
block_ops;
206 Teuchos::RCP<Teuchos::ParameterList>
params;
224 Teuchos::RCP<Epetra_CrsMatrix>
mean;
239 Teuchos::Array< Teuchos::RCP<Epetra_CrsMatrix> >
kl_blocks;
242 Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly >
kl_ops;
An Epetra operator representing the block stochastic Galerkin operator.
virtual bool UseTranspose() const
Returns the current UseTranspose setting.
double drop_tolerance
Tolerance for dropping entries in sparse 3 tensor.
Teuchos::RCP< Epetra_Map > block_vec_map
Block map for vectorized-matrices.
Teuchos::RCP< const Stokhos::OrthogPolyBasis< int, double > > sg_basis
Stochastic Galerking basis.
Teuchos::RCP< Stokhos::Sparse3Tensor< int, double > > sparse_kl_coeffs
Sparse KL coefficients.
virtual void setupOperator(const Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > &poly)
Setup operator.
virtual bool HasNormInf() const
Returns true if the this object can provide an approximate Inf-norm, false otherwise.
virtual ~KLReducedMatrixFreeOperator()
Destructor.
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 ...
int num_KL
Number of KL terms.
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< Teuchos::ParameterList > params
Algorithmic parameters.
Teuchos::RCP< const Cijk_type > Cijk
Stores triple product tensor.
Teuchos::RCP< Stokhos::EpetraVectorOrthogPoly > block_vec_poly
Polynomial sorting vectorized matrix coefficients.
Stokhos::Sparse3Tensor< int, double > Cijk_type
Short-hand for Cijk.
KLReducedMatrixFreeOperator(const KLReducedMatrixFreeOperator &)
Private to prohibit copying.
void setup()
Setup KL blocks.
Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > block_ops
Stores operators.
std::string label
Label for operator.
int expansion_size
Number of terms in expansion.
virtual const Epetra_Map & OperatorDomainMap() const
Returns the Epetra_Map object associated with the domain of this matrix operator.
Teuchos::RCP< const Epetra_Map > range_base_map
Stores range base map.
Teuchos::Array< Teuchos::Array< double > > dot_products
Dot products of KL eigenvectors and Jacobian blocks.
Teuchos::RCP< const EpetraExt::MultiComm > sg_comm
Stores SG parallel communicator.
virtual const Epetra_Comm & Comm() const
Returns a reference to the Epetra_Comm communicator associated with this operator.
Teuchos::RCP< const Stokhos::EpetraSparse3Tensor > epetraCijk
Stores Epetra Cijk tensor.
bool do_error_tests
Whether to do KL error tests (can be expensive)
virtual Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > getSGPolynomial()
Get SG polynomial.
int num_blocks
Number of blocks.
virtual double NormInf() const
Returns an approximate infinity norm of the operator matrix.
Teuchos::RCP< Stokhos::MatrixFreeOperator > kl_mat_free_op
Matrix-Free operator using KL operators.
KLReducedMatrixFreeOperator(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 > &domain_base_map, const Teuchos::RCP< const Epetra_Map > &range_base_map, const Teuchos::RCP< const Epetra_Map > &domain_sg_map, const Teuchos::RCP< const Epetra_Map > &range_sg_map, 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.
KLReducedMatrixFreeOperator & operator=(const KLReducedMatrixFreeOperator &)
Private to prohibit copying.
bool useTranspose
Flag indicating whether transpose was selected.
virtual const char * Label() const
Returns a character string describing the operator.
int num_KL_computed
Number of computed KL terms.
Teuchos::RCP< const Epetra_Map > range_sg_map
Stores range SG map.
Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > kl_ops
KL blocks as operators.
Teuchos::RCP< const Epetra_Map > domain_base_map
Stores domain base map.
Teuchos::RCP< Epetra_CrsMatrix > mean
Mean block.
Teuchos::Array< Teuchos::RCP< Epetra_CrsMatrix > > kl_blocks
KL blocks.
virtual int SetUseTranspose(bool UseTranspose)
Set to true if the transpose of the operator is requested.
Teuchos::RCP< const Epetra_Map > domain_sg_map
Stores domain SG map.
Abstract base class for multivariate orthogonal polynomials.
An abstract class to represent a generic stochastic Galerkin operator 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.