42 #ifndef EPETRAEXT_PRODUCT_OPERATOR_H
43 #define EPETRAEXT_PRODUCT_OPERATOR_H
316 const char *
Label()
const;
335 typedef std::vector<Teuchos::RCP<const Epetra_Operator> >
Op_t;
338 typedef std::vector<Teuchos::RCP<Epetra_Vector> >
EV_t;
402 Op_.size()==0, std::logic_error
403 ,
"Epetra::ProductOperator: Error, Client has not called initialize(...) yet!"
411 k < 0 || static_cast<int>(
Op_.size())-1 < k, std::logic_error
412 ,
"Epetra::ProductOperator: Error, k = "<<k<<
" is not in the range [0,"<<
Op_.size()-1<<
"]!"
418 #endif // EPETRAEXT_PRODUCT_OPERATOR_H
const Epetra_Map & OperatorDomainMap() const
std::vector< Teuchos::ETransp > Op_trans_t
bool UseTranspose() const
int SetUseTranspose(bool UseTranspose)
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
void initialize(const int num_Op, const Teuchos::RCP< const Epetra_Operator > Op[], const Teuchos::ETransp Op_trans[], const EApplyMode Op_inverse[])
Setup with constituent operators.
Implements Epetra_Operator as a product of one or more Epetra_Operator objects.
const char * Label() const
int ApplyInverse(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
void initializeTempVecs(bool applyInverse) const
ProductOperator()
Construct to uninitialized.
const Epetra_Comm & Comm() const
void validateIndex(int k) const
int num_Op() const
Return the number of aggregate opeators.
const Epetra_Map & OperatorRangeMap() const
Teuchos::RCP< const Epetra_Operator > Op(int k) const
Access the kth operator (zero-based).
Teuchos::ETransp Op_trans(int k) const
Access the transpose mode of the kth operator (zero-based).
int Apply(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
std::vector< EApplyMode > Op_inverse_t
void assertInitialized() const
EApplyMode Op_inverse(int k) const
Access the inverse mode of the kth operator (zero-based).
void uninitialize(int *num_Op, Teuchos::RCP< const Epetra_Operator > Op[], Teuchos::ETransp Op_trans[], EApplyMode p_inverse[])
Set to an uninitialized state and wipe out memory.
void applyConstituent(const int k, Teuchos::ETransp Op_trans, EApplyMode Op_inverse, const Epetra_MultiVector &X_k, Epetra_MultiVector *Y_k) const
Apply the kth aggregate operator M[k] correctly.
std::vector< Teuchos::RCP< const Epetra_Operator > > Op_t
std::vector< Teuchos::RCP< Epetra_Vector > > EV_t