46 #ifndef MUELU_UTILITIES_DEF_HPP
47 #define MUELU_UTILITIES_DEF_HPP
49 #include <Teuchos_DefaultComm.hpp>
56 #ifdef HAVE_MUELU_EPETRA
58 #include "Epetra_MpiComm.h"
62 #if defined(HAVE_MUELU_EPETRA) && defined(HAVE_MUELU_EPETRAEXT)
75 #include <Tpetra_RowMatrixTransposer.hpp>
76 #include <TpetraExt_MatrixMatrix.hpp>
77 #include <Xpetra_TpetraMultiVector.hpp>
79 #include <Xpetra_TpetraCrsMatrix.hpp>
80 #include <Xpetra_TpetraBlockCrsMatrix.hpp>
82 #ifdef HAVE_MUELU_EPETRA
89 #include <Xpetra_Map.hpp>
90 #include <Xpetra_MapFactory.hpp>
94 #include <Xpetra_MultiVectorFactory.hpp>
102 #if defined(HAVE_MUELU_EPETRA) && defined(HAVE_MUELU_ML)
103 #include <ml_operator.h>
104 #include <ml_epetra_utils.h>
109 #ifdef HAVE_MUELU_EPETRA
114 #ifdef HAVE_MUELU_EPETRA
115 template <
typename SC,
typename LO,
typename GO,
typename NO>
117 return Xpetra::Convert_Epetra_CrsMatrix_ToXpetra_CrsMatrixWrap<SC, LO, GO, NO>(epAB);
121 #ifdef HAVE_MUELU_EPETRA
122 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
125 if (tmpVec == Teuchos::null)
126 throw Exceptions::BadCast(
"Cast from Xpetra::MultiVector to Xpetra::EpetraMultiVector failed");
127 return tmpVec->getEpetra_MultiVector();
130 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
133 if (tmpVec == Teuchos::null)
134 throw Exceptions::BadCast(
"Cast from Xpetra::MultiVector to Xpetra::EpetraMultiVector failed");
135 return tmpVec->getEpetra_MultiVector();
138 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
144 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
150 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
153 if (crsOp == Teuchos::null)
156 if (tmp_ECrsMtx == Teuchos::null)
157 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::EpetraCrsMatrix failed");
158 return tmp_ECrsMtx->getEpetra_CrsMatrix();
161 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
164 if (crsOp == Teuchos::null)
167 if (tmp_ECrsMtx == Teuchos::null)
168 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::EpetraCrsMatrix failed");
169 return tmp_ECrsMtx->getEpetra_CrsMatrixNonConst();
172 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
179 }
catch (std::bad_cast&) {
180 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::EpetraCrsMatrix failed");
182 }
catch (std::bad_cast&) {
187 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
194 }
catch (std::bad_cast&) {
195 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::EpetraCrsMatrix failed");
197 }
catch (std::bad_cast&) {
202 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
205 if (xeMap == Teuchos::null)
206 throw Exceptions::BadCast(
"Utilities::Map2EpetraMap : Cast from Xpetra::Map to Xpetra::EpetraMap failed");
207 return xeMap->getEpetra_Map();
211 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
215 if (tmpVec == Teuchos::null)
216 throw Exceptions::BadCast(
"Cast from Xpetra::MultiVector to Xpetra::TpetraMultiVector failed");
217 return tmpVec->getTpetra_MultiVector();
220 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
223 if (tmpVec == Teuchos::null)
224 throw Exceptions::BadCast(
"Cast from Xpetra::MultiVector to Xpetra::TpetraMultiVector failed");
225 return tmpVec->getTpetra_MultiVector();
228 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
234 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
240 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
247 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
251 if (crsOp == Teuchos::null)
254 if (tmp_ECrsMtx == Teuchos::null)
255 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix failed");
256 return tmp_ECrsMtx->getTpetra_CrsMatrix();
259 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
262 if (crsOp == Teuchos::null)
265 if (tmp_ECrsMtx == Teuchos::null)
266 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix failed");
267 return tmp_ECrsMtx->getTpetra_CrsMatrixNonConst();
270 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
277 }
catch (std::bad_cast&) {
278 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix failed");
280 }
catch (std::bad_cast&) {
285 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
292 }
catch (std::bad_cast&) {
293 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix failed");
295 }
catch (std::bad_cast&) {
300 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
305 if (crsOp == Teuchos::null)
308 if (tmp_ECrsMtx == Teuchos::null)
309 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraBlockCrsMatrix failed");
310 return tmp_ECrsMtx->getTpetra_BlockCrsMatrix();
313 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
317 if (crsOp == Teuchos::null)
320 if (tmp_ECrsMtx == Teuchos::null)
321 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraBlockCrsMatrix failed");
322 return tmp_ECrsMtx->getTpetra_BlockCrsMatrixNonConst();
325 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
329 const XCrsMatrixWrap& crsOp =
dynamic_cast<const XCrsMatrixWrap&
>(Op);
333 }
catch (std::bad_cast&) {
334 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraBlockCrsMatrix failed");
336 }
catch (std::bad_cast&) {
341 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
345 XCrsMatrixWrap& crsOp =
dynamic_cast<XCrsMatrixWrap&
>(Op);
349 }
catch (std::bad_cast&) {
350 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraBlockCrsMatrix failed");
352 }
catch (std::bad_cast&) {
357 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
363 if (crsOp == Teuchos::null)
370 return tmp_Crs->getTpetra_CrsMatrixNonConst();
373 if (tmp_BlockCrs.is_null())
374 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix and Xpetra::TpetraBlockCrsMatrix failed");
375 return tmp_BlockCrs->getTpetra_BlockCrsMatrixNonConst();
377 }
else if (!rmat.is_null()) {
378 return rmat->getTpetra_RowMatrix();
387 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
393 if (crsOp == Teuchos::null)
400 return tmp_Crs->getTpetra_CrsMatrixNonConst();
403 if (tmp_BlockCrs.is_null())
404 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix and Xpetra::TpetraBlockCrsMatrix failed");
405 return tmp_BlockCrs->getTpetra_BlockCrsMatrixNonConst();
407 }
else if (!rmat.is_null()) {
408 return rmat->getTpetra_RowMatrixNonConst();
417 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
420 if (tmp_TMap == Teuchos::null)
421 throw Exceptions::BadCast(
"Utilities::Map2TpetraMap : Cast from Xpetra::Map to Xpetra::TpetraMap failed");
422 return tmp_TMap->getTpetra_Map();
425 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
428 bool doOptimizeStorage) {
432 for (
int i = 0; i < scalingVector.
size(); ++i)
433 sv[i] = one / scalingVector[i];
435 for (
int i = 0; i < scalingVector.
size(); ++i)
436 sv[i] = scalingVector[i];
439 switch (Op.getRowMap()->lib()) {
441 MyOldScaleMatrix_Tpetra(Op, sv, doFillComplete, doOptimizeStorage);
445 MyOldScaleMatrix_Epetra(Op, sv, doFillComplete, doOptimizeStorage);
453 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
458 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
461 bool doOptimizeStorage) {
470 if (maxRowSize == Teuchos::as<size_t>(-1))
480 for (
size_t i = 0; i < rowMap->getLocalNumElements(); ++i) {
486 for (
size_t j = 0; j < nnz; ++j) {
487 scaledVals[j] = scalingVector[i] * vals[j];
499 for (
size_t i = 0; i < rowMap->getLocalNumElements(); ++i) {
506 for (
size_t j = 0; j < nnz; ++j)
507 scaledVals[j] = scalingVector[i] * vals[j];
515 if (doFillComplete) {
516 if (domainMap == Teuchos::null || rangeMap == Teuchos::null)
517 throw Exceptions::RuntimeError(
"In Utilities::Scaling: cannot fillComplete because the domain and/or range map hasn't been defined");
520 params->
set(
"Optimize Storage", doOptimizeStorage);
521 params->
set(
"No Nonlocal Changes",
true);
529 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
533 #if defined(HAVE_MUELU_EPETRA) && defined(HAVE_MUELU_EPETRAEXT)
534 std::string TorE =
"epetra";
536 std::string TorE =
"tpetra";
539 #if defined(HAVE_MUELU_EPETRA) && defined(HAVE_MUELU_EPETRAEXT)
548 if (TorE ==
"tpetra") {
551 if (Helpers::isTpetraCrs(Op)) {
561 transposeParams->
set(
"sort",
false);
568 if (!AAAA->isFillComplete())
571 if (Op.IsView(
"stridedMaps"))
572 AAAA->CreateView(
"stridedMaps", Teuchos::rcpFromRef(Op),
true );
575 }
else if (Helpers::isTpetraBlockCrs(Op)) {
590 transposeParams->
set(
"sort",
false);
598 if (Op.IsView(
"stridedMaps"))
599 AAAA->CreateView(
"stridedMaps", Teuchos::rcpFromRef(Op),
true );
608 std::cout <<
"Utilities::Transpose() not implemented for Epetra" << std::endl;
609 return Teuchos::null;
613 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
618 #if defined(HAVE_XPETRA_TPETRA) && (defined(HAVE_TPETRA_INST_COMPLEX_DOUBLE) || defined(HAVE_TPETRA_INST_COMPLEX_FLOAT))
619 using range_type = Kokkos::RangePolicy<LocalOrdinal, typename Node::execution_space>;
622 if ((
typeid(
Scalar).name() ==
typeid(std::complex<double>).name()) ||
623 (
typeid(
Scalar).name() ==
typeid(std::complex<float>).name())) {
624 size_t numVecs = X->getNumVectors();
626 auto XVec = X->getDeviceLocalView(Xpetra::Access::ReadOnly);
627 auto XVecScalar = Xscalar->getDeviceLocalView(Xpetra::Access::ReadWrite);
629 Kokkos::parallel_for(
630 "MueLu:Utils::RealValuedToScalarMultiVector", range_type(0, X->getLocalLength()),
631 KOKKOS_LAMBDA(
const size_t i) {
632 for (
size_t j = 0; j < numVecs; j++)
633 XVecScalar(i, j) = XVec(i, j);
641 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
654 #if !defined(HAVE_TPETRA_EXPLICIT_INSTANTIATION) || defined(HAVE_TPETRA_INST_FLOAT)
662 #if !defined(HAVE_TPETRA_EXPLICIT_INSTANTIATION) || defined(HAVE_TPETRA_INST_DOUBLE)
668 paramList.
remove(
"Coordinates");
670 #if !defined(HAVE_TPETRA_EXPLICIT_INSTANTIATION) || defined(HAVE_TPETRA_INST_FLOAT)
674 paramList.
remove(
"Coordinates");
675 doubleCoords =
rcp(
new tdMV(floatCoords->getMap(), floatCoords->getNumVectors()));
680 if (doubleCoords != Teuchos::null) {
689 throw Exceptions::RuntimeError(
"ExtractCoordinatesFromParameterList: The coordinates vector in parameter list is expected to be a Tpetra multivector with SC=double or float.");
693 if (paramList.
isType<decltype(coordinates)>(
"Coordinates")) {
694 coordinates = paramList.
get<decltype(coordinates)>(
"Coordinates");
702 #define MUELU_UTILITIES_SHORT
703 #endif // MUELU_UTILITIES_DEF_HPP
RCP< CrsMatrix > getCrsMatrix() const
Exception indicating invalid cast attempted.
static const RCP< const Tpetra::Map< LocalOrdinal, GlobalOrdinal, Node > > Map2TpetraMap(const Xpetra::Map< LocalOrdinal, GlobalOrdinal, Node > &map)
MueLu::DefaultLocalOrdinal LocalOrdinal
static RCP< const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > MV2TpetraMV(RCP< Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node >> const vec)
Helper utility to pull out the underlying Tpetra objects from an Xpetra object.
static void MyOldScaleMatrix(Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &Op, const Teuchos::ArrayRCP< const Scalar > &scalingVector, bool doInverse=true, bool doFillComplete=true, bool doOptimizeStorage=true)
static Teuchos::RCP< MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Build(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node >> &map, size_t NumVectors, bool zeroOut=true)
static RCP< Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > RealValuedToScalarMultiVector(RCP< Xpetra::MultiVector< typename Teuchos::ScalarTraits< Scalar >::coordinateType, LocalOrdinal, GlobalOrdinal, Node >> X)
static const Epetra_Map & Map2EpetraMap(const Xpetra::Map< LocalOrdinal, GlobalOrdinal, Node > &map)
Teuchos::RCP< const map_type > getRangeMap() const override
RCP< Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getTpetra_BlockCrsMatrixNonConst() const
T & get(const std::string &name, T def_value)
RCP< Epetra_CrsMatrix > getEpetra_CrsMatrixNonConst() const
static RCP< Xpetra::MultiVector< typename Teuchos::ScalarTraits< Scalar >::magnitudeType, LocalOrdinal, GlobalOrdinal, Node > > ExtractCoordinatesFromParameterList(ParameterList ¶mList)
ParameterList & set(std::string const &name, T const &value, std::string const &docString="", RCP< const ParameterEntryValidator > const &validator=null)
static RCP< Tpetra::RowMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op2NonConstTpetraRow(RCP< Xpetra::Operator< Scalar, LocalOrdinal, GlobalOrdinal, Node >> Op)
bool is_null(const std::shared_ptr< T > &p)
static RCP< Epetra_CrsMatrix > Op2NonConstEpetraCrs(RCP< Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >> Op)
void resumeFill(const Teuchos::RCP< Teuchos::ParameterList > ¶ms=Teuchos::null)
static RCP< const Tpetra::CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op2TpetraCrs(RCP< const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >> Op)
static RCP< Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op2NonConstTpetraBlockCrs(RCP< Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >> Op)
bool isParameter(const std::string &name) const
void deep_copy(MultiVector< DS, DL, DG, DN > &dst, const MultiVector< SS, SL, SG, SN > &src)
bool remove(std::string const &name, bool throwIfNotExists=true)
virtual void fillComplete(const RCP< const Map > &domainMap, const RCP< const Map > &rangeMap, const RCP< ParameterList > ¶ms=null)=0
size_t getLocalMaxNumRowEntries() const override
bool isFillComplete() const override
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
static RCP< Epetra_MultiVector > MV2NonConstEpetraMV(RCP< Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node >> vec)
MueLu::DefaultScalar Scalar
MueLu::DefaultGlobalOrdinal GlobalOrdinal
RCP< const Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getTpetra_BlockCrsMatrix() const
static RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > MV2NonConstTpetraMV2(Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &vec)
static RCP< Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Transpose(Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &Op, bool optimizeTranspose=false, const std::string &label=std::string(), const Teuchos::RCP< Teuchos::ParameterList > ¶ms=Teuchos::null)
RCP< const Epetra_CrsMatrix > getEpetra_CrsMatrix() const
local_ordinal_type replaceGlobalValues(const global_ordinal_type globalRow, const Kokkos::View< const global_ordinal_type *, Kokkos::AnonymousSpace > &inputInds, const Kokkos::View< const impl_scalar_type *, Kokkos::AnonymousSpace > &inputVals)
RCP< Xpetra::CrsMatrixWrap< SC, LO, GO, NO > > Convert_Epetra_CrsMatrix_ToXpetra_CrsMatrixWrap(RCP< Epetra_CrsMatrix > &epAB)
size_t getNumEntriesInLocalRow(local_ordinal_type localRow) const override
Teuchos::RCP< bcrs_matrix_type > createTranspose(const Teuchos::RCP< Teuchos::ParameterList > ¶ms=Teuchos::null)
Teuchos::RCP< crs_matrix_type > createTranspose(const Teuchos::RCP< Teuchos::ParameterList > ¶ms=Teuchos::null)
static void MyOldScaleMatrix_Epetra(Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &Op, const Teuchos::ArrayRCP< Scalar > &scalingVector, bool doFillComplete, bool doOptimizeStorage)
Teuchos::RCP< const map_type > getDomainMap() const override
virtual bool isLocallyIndexed() const =0
RCP< const Tpetra::CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getTpetra_CrsMatrix() const
RCP< Tpetra::CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getTpetra_CrsMatrixNonConst() const
bool isType(const std::string &name) const
static RCP< const Tpetra::RowMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op2TpetraRow(RCP< const Xpetra::Operator< Scalar, LocalOrdinal, GlobalOrdinal, Node >> Op)
virtual void getGlobalRowView(GlobalOrdinal GlobalRow, global_inds_host_view_type &indices, values_host_view_type &values) const =0
static RCP< const Epetra_MultiVector > MV2EpetraMV(RCP< Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node >> const vec)
Helper utility to pull out the underlying Epetra objects from an Xpetra object.
RCP< Epetra_MultiVector > getEpetra_MultiVector() const
virtual Teuchos::RCP< const Map > getRangeMap() const =0
RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getTpetra_MultiVector() const
size_t getNumEntriesInGlobalRow(GlobalOrdinal globalRow) const override
Exception throws to report errors in the internal logical of the program.
virtual void getLocalRowView(LocalOrdinal LocalRow, local_inds_host_view_type &indices, values_host_view_type &values) const =0
local_ordinal_type replaceLocalValues(const local_ordinal_type localRow, const Kokkos::View< const local_ordinal_type *, Kokkos::AnonymousSpace > &inputInds, const Kokkos::View< const impl_scalar_type *, Kokkos::AnonymousSpace > &inputVals)
static RCP< const Tpetra::BlockCrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op2TpetraBlockCrs(RCP< const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >> Op)
static void MyOldScaleMatrix_Tpetra(Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &Op, const Teuchos::ArrayRCP< Scalar > &scalingVector, bool doFillComplete, bool doOptimizeStorage)
static RCP< const Epetra_CrsMatrix > Op2EpetraCrs(RCP< const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >> Op)
#define TEUCHOS_TEST_FOR_EXCEPT(throw_exception_test)
static RCP< Tpetra::CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op2NonConstTpetraCrs(RCP< Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >> Op)
virtual Teuchos::RCP< const Map > getDomainMap() const =0
static RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > MV2NonConstTpetraMV(RCP< Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node >> vec)
Teuchos::RCP< const map_type > getRowMap() const override