44#ifndef STOKHOS_INTERLACED_OPERATOR_HPP
45#define STOKHOS_INTERLACED_OPERATOR_HPP
48#include "EpetraExt_MultiComm.h"
51#include "Teuchos_ParameterList.hpp"
52#include "EpetraExt_BlockCrsMatrix.h"
66 public EpetraExt::BlockCrsMatrix {
72 const Teuchos::RCP<const EpetraExt::MultiComm>&
sg_comm,
74 const Teuchos::RCP<const Stokhos::EpetraSparse3Tensor>&
epetraCijk,
75 const Teuchos::RCP<const Epetra_CrsGraph>& base_graph,
76 const Teuchos::RCP<Teuchos::ParameterList>& params);
86 const Teuchos::RCP<Stokhos::EpetraOperatorOrthogPoly >& poly);
89 virtual Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly >
93 virtual Teuchos::RCP<const Stokhos::EpetraOperatorOrthogPoly >
112 Teuchos::RCP<const EpetraExt::MultiComm>
sg_comm;
115 Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> >
sg_basis;
124 Teuchos::RCP<const Cijk_type>
Cijk;
127 Teuchos::RCP<Stokhos::EpetraOperatorOrthogPoly >
block_ops;
An Epetra operator representing the block stochastic Galerkin operator generated by fully assembling ...
Teuchos::RCP< const EpetraExt::MultiComm > sg_comm
Stores SG parallel communicator.
Stokhos::Sparse3Tensor< int, double > Cijk_type
Short-hand for Cijk.
Teuchos::RCP< const Cijk_type > Cijk
Stores triple product tensor.
bool scale_op
Flag indicating whether operator be scaled with <\psi_i^2>
Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > block_ops
Stores operators.
Teuchos::RCP< const Stokhos::OrthogPolyBasis< int, double > > sg_basis
Stochastic Galerking basis.
Teuchos::RCP< const Stokhos::EpetraSparse3Tensor > epetraCijk
Stores Epetra Cijk tensor.
void SumIntoGlobalBlock_Deterministic(double alpha, const Epetra_RowMatrix &determBlock, int Row, int Col)
Sum into global matrix.
InterlacedOperator(const InterlacedOperator &)
Private to prohibit copying.
virtual void setupOperator(const Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > &poly)
Setup operator.
virtual Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > getSGPolynomial()
Get SG polynomial.
InterlacedOperator(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_CrsGraph > &base_graph, const Teuchos::RCP< Teuchos::ParameterList > ¶ms)
Constructor.
InterlacedOperator & operator=(const InterlacedOperator &)
Private to prohibit copying.
virtual ~InterlacedOperator()
Destructor.
bool only_use_linear
Flag indicating whether to only use linear terms.
bool include_mean
Flag indicating whether to include mean term.
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.