47 #include "Teko_TpetraReorderedMappingStrategy.hpp"
49 #include "Teko_BlockedReordering.hpp"
53 using Teuchos::rcp_dynamic_cast;
54 using Teuchos::rcpFromRef;
57 namespace TpetraHelpers {
59 TpetraReorderedMappingStrategy::TpetraReorderedMappingStrategy(
60 const BlockReorderManager& brm,
const Teuchos::RCP<const MappingStrategy>& map)
61 : reorderManager_(brm), mapStrategy_(map) {
62 rangeMap_ = mapStrategy_->rangeMap();
63 domainMap_ = mapStrategy_->domainMap();
66 void TpetraReorderedMappingStrategy::copyTpetraIntoThyra(
67 const Tpetra::MultiVector<ST, LO, GO, NT>& X,
68 const Teuchos::Ptr<Thyra::MultiVectorBase<ST> >& thyra_X)
const {
69 using Teuchos::ptr_const_cast;
70 using Teuchos::rcp_const_cast;
73 RCP<Thyra::ProductMultiVectorBase<ST> > prod_X =
74 rcp_dynamic_cast<Thyra::ProductMultiVectorBase<ST> >(rcpFromRef(*thyra_X));
75 RCP<Thyra::MultiVectorBase<ST> > flat_X = buildFlatMultiVector(reorderManager_, prod_X);
78 mapStrategy_->copyTpetraIntoThyra(X, flat_X.ptr());
81 void TpetraReorderedMappingStrategy::copyThyraIntoTpetra(
82 const RCP<
const Thyra::MultiVectorBase<ST> >& thyra_Y,
83 Tpetra::MultiVector<ST, LO, GO, NT>& Y)
const {
85 RCP<const Thyra::ProductMultiVectorBase<ST> > prod_Y =
86 rcp_dynamic_cast<
const Thyra::ProductMultiVectorBase<ST> >(rcpFromRef(*thyra_Y));
87 RCP<const Thyra::MultiVectorBase<ST> > flat_Y = buildFlatMultiVector(reorderManager_, prod_Y);
90 mapStrategy_->copyThyraIntoTpetra(flat_Y, Y);