47 #ifndef __Teko_BlockedMappingStrategy_hpp__
48 #define __Teko_BlockedMappingStrategy_hpp__
54 #include "Teuchos_RCP.hpp"
57 #include "Thyra_MultiVectorBase.hpp"
58 #include "Thyra_LinearOpBase.hpp"
59 #include "Thyra_BlockedLinearOpBase.hpp"
62 #include "Epetra_Comm.h"
63 #include "Epetra_Map.h"
64 #include "Epetra_MultiVector.h"
65 #include "Epetra_Import.h"
66 #include "Epetra_Export.h"
67 #include "Epetra_CrsMatrix.h"
70 #include "Teko_EpetraOperatorWrapper.hpp"
71 #include "Teko_BlockingEpetra.hpp"
76 class BlockedMappingStrategy :
public MappingStrategy {
90 BlockedMappingStrategy(
const std::vector<std::vector<int> >& vars,
91 const Teuchos::RCP<const Epetra_Map>& baseMap,
const Epetra_Comm& comm);
104 virtual void copyEpetraIntoThyra(
105 const Epetra_MultiVector& epetra_X,
106 const Teuchos::Ptr<Thyra::MultiVectorBase<double> >& thyra_X)
const;
116 const Teuchos::RCP<
const Thyra::MultiVectorBase<double> >& thyra_Y,
117 Epetra_MultiVector& epetra_Y)
const;
126 virtual const Teuchos::RCP<const Epetra_Map> domainMap()
const {
return domainMap_; }
135 virtual const Teuchos::RCP<const Epetra_Map> rangeMap()
const {
return rangeMap_; }
141 virtual std::string toString()
const {
return std::string(
"BlockedMappingStrategy"); }
162 void buildBlockTransferData(
const std::vector<std::vector<int> >& vars,
163 const Teuchos::RCP<const Epetra_Map>& baseMap,
164 const Epetra_Comm& comm);
177 const std::vector<Blocking::MapPair>& getMaps()
const {
return blockMaps_; }
190 const Teuchos::RCP<Thyra::BlockedLinearOpBase<double> > buildBlockedThyraOp(
191 const Teuchos::RCP<const Epetra_CrsMatrix>& mat,
const std::string& label =
"<ANYM>")
const;
202 void rebuildBlockedThyraOp(
const RCP<const Epetra_CrsMatrix>& mat,
203 const RCP<Thyra::BlockedLinearOpBase<double> >& A)
const;
212 Teuchos::RCP<const Epetra_Map> domainMap_;
213 Teuchos::RCP<const Epetra_Map> rangeMap_;
218 std::vector<Blocking::MapPair> blockMaps_;
219 std::vector<Teuchos::RCP<Epetra_Import> > blockImport_;
220 std::vector<Teuchos::RCP<Epetra_Export> > blockExport_;
virtual void copyThyraIntoEpetra(const RCP< const Thyra::MultiVectorBase< double > > &thyraX, Epetra_MultiVector &epetraX) const =0
Copy an Thyra::MultiVectorBase into a Epetra_MultiVector.