47 #ifndef __Teko_StridedMappingStrategy_hpp__
48 #define __Teko_StridedMappingStrategy_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"
75 class StridedMappingStrategy :
public MappingStrategy {
89 StridedMappingStrategy(
const std::vector<int>& vars,
const Teuchos::RCP<const Epetra_Map>& map,
90 const Epetra_Comm& comm);
103 virtual void copyEpetraIntoThyra(
104 const Epetra_MultiVector& epetra_X,
105 const Teuchos::Ptr<Thyra::MultiVectorBase<double> >& thyra_X)
const;
115 const Teuchos::RCP<
const Thyra::MultiVectorBase<double> >& thyra_Y,
116 Epetra_MultiVector& epetra_Y)
const;
125 virtual const Teuchos::RCP<const Epetra_Map> domainMap()
const {
return domainMap_; }
134 virtual const Teuchos::RCP<const Epetra_Map> rangeMap()
const {
return rangeMap_; }
140 virtual std::string toString()
const {
return std::string(
"StridedMappingStrategy"); }
161 void buildBlockTransferData(
const std::vector<int>& vars,
162 const Teuchos::RCP<const Epetra_Map>& baseMap,
163 const Epetra_Comm& comm);
176 const std::vector<std::pair<int, Teuchos::RCP<Epetra_Map> > >& getMaps()
const {
191 const Teuchos::RCP<Thyra::BlockedLinearOpBase<double> > buildBlockedThyraOp(
192 const Teuchos::RCP<const Epetra_CrsMatrix>& mat,
const std::string& label =
"<ANYM>")
const;
203 void rebuildBlockedThyraOp(
const RCP<const Epetra_CrsMatrix>& mat,
204 const RCP<Thyra::BlockedLinearOpBase<double> >& A)
const;
213 Teuchos::RCP<const Epetra_Map> domainMap_;
214 Teuchos::RCP<const Epetra_Map> rangeMap_;
219 std::vector<std::pair<int, Teuchos::RCP<Epetra_Map> > > blockMaps_;
220 std::vector<Teuchos::RCP<Epetra_Export> > blockExport_;
221 std::vector<Teuchos::RCP<Epetra_Import> > blockImport_;
virtual void copyThyraIntoEpetra(const RCP< const Thyra::MultiVectorBase< double > > &thyraX, Epetra_MultiVector &epetraX) const =0
Copy an Thyra::MultiVectorBase into a Epetra_MultiVector.