47 #ifndef __Teko_BlockedEpetraOperator_hpp__
48 #define __Teko_BlockedEpetraOperator_hpp__
51 #include "Epetra_Operator.h"
54 #include "Teuchos_RCP.hpp"
56 #include "Thyra_LinearOpBase.hpp"
59 #include "Teko_BlockedReordering.hpp"
60 #include "Teko_EpetraOperatorWrapper.hpp"
61 #include "Teko_BlockedMappingStrategy.hpp"
80 const Teuchos::RCP<const Epetra_Operator> &content,
81 const std::string &label =
"<ANYM>");
91 virtual void SetContent(
const std::vector<std::vector<int> > &vars,
92 const Teuchos::RCP<const Epetra_Operator> &content);
99 virtual const Teuchos::RCP<const Epetra_Operator> GetContent()
const {
return fullContent_; }
101 virtual const Teuchos::RCP<const Epetra_Operator> GetContent() {
return fullContent_; }
103 const Teuchos::RCP<const Epetra_Operator> GetBlock(
int i,
int j)
const;
108 void Reorder(
const BlockReorderManager &brm);
115 virtual void WriteBlocks(
const std::string &prefix)
const;
126 virtual int SetUseTranspose(
bool ) {
return -1; }
128 virtual int ApplyInverse(
const Epetra_MultiVector & , Epetra_MultiVector & )
const {
129 TEUCHOS_ASSERT(
false);
133 virtual double NormInf()
const {
134 TEUCHOS_ASSERT(
false);
139 virtual bool UseTranspose()
const {
return false; }
140 virtual bool HasNormInf()
const {
return false; }
141 virtual const Epetra_Comm &Comm()
const {
return fullContent_->Comm(); }
148 Teuchos::RCP<const Epetra_Operator> fullContent_;
149 Teuchos::RCP<BlockedMappingStrategy> blockedMapping_;
150 Teuchos::RCP<Thyra::LinearOpBase<double> > blockedOperator_;
151 Teuchos::RCP<const BlockReorderManager> reorderManager_;
155 void BuildBlockedOperator();
virtual void RebuildOps()
bool testAgainstFullOperator(int count, double tol) const
Helps perform sanity checks.
void Reorder(const BlockReorderManager &brm)
Implements the Epetra_Operator interface with a Thyra LinearOperator. This enables the use of absrtac...
virtual void SetContent(const std::vector< std::vector< int > > &vars, const Teuchos::RCP< const Epetra_Operator > &content)
void RemoveReording()
Remove any reordering on this object.
Tear about a user specified Epetra_Operator (CrsMatrix) using a vector of vectors of GIDs for each bl...
BlockedEpetraOperator(const std::vector< std::vector< int > > &vars, const Teuchos::RCP< const Epetra_Operator > &content, const std::string &label="<ANYM>")
virtual void WriteBlocks(const std::string &prefix) const