46 #ifndef MUELU_UTILITIES_DECL_HPP 47 #define MUELU_UTILITIES_DECL_HPP 54 #include <Teuchos_DefaultComm.hpp> 55 #include <Teuchos_ScalarTraits.hpp> 56 #include <Teuchos_ParameterList.hpp> 58 #ifdef HAVE_MUELU_TPETRA 59 #include <Xpetra_TpetraBlockCrsMatrix.hpp> 61 #include <Xpetra_BlockedCrsMatrix_fwd.hpp> 62 #include <Xpetra_CrsMatrix_fwd.hpp> 63 #include <Xpetra_CrsMatrixWrap_fwd.hpp> 64 #include <Xpetra_Map_fwd.hpp> 65 #include <Xpetra_MapFactory_fwd.hpp> 66 #include <Xpetra_Matrix_fwd.hpp> 67 #include <Xpetra_MatrixFactory_fwd.hpp> 68 #include <Xpetra_MultiVector_fwd.hpp> 69 #include <Xpetra_MultiVectorFactory_fwd.hpp> 70 #include <Xpetra_Operator_fwd.hpp> 71 #include <Xpetra_Vector_fwd.hpp> 72 #include <Xpetra_VectorFactory_fwd.hpp> 73 #include <Xpetra_ExportFactory.hpp> 75 #include <Xpetra_Import.hpp> 76 #include <Xpetra_ImportFactory.hpp> 77 #include <Xpetra_MatrixMatrix.hpp> 79 #ifdef HAVE_MUELU_EPETRA 80 #include <Xpetra_EpetraCrsMatrix_fwd.hpp> 84 #include <Xpetra_EpetraCrsMatrix.hpp> 85 #include <Xpetra_CrsMatrixWrap.hpp> 91 #ifdef HAVE_MUELU_EPETRAEXT 92 class Epetra_CrsMatrix;
93 class Epetra_MultiVector;
98 #ifdef HAVE_MUELU_TPETRA 99 #include <Tpetra_CrsMatrix.hpp> 100 #include <Tpetra_RowMatrixTransposer.hpp> 101 #include <Tpetra_Map.hpp> 102 #include <Tpetra_MultiVector.hpp> 103 #include <Xpetra_TpetraCrsMatrix_fwd.hpp> 104 #include <Xpetra_TpetraMultiVector_fwd.hpp> 107 #include <MueLu_UtilitiesBase.hpp> 112 #ifdef HAVE_MUELU_EPETRA 114 template<
typename SC,
typename LO,
typename GO,
typename NO>
115 RCP<Xpetra::CrsMatrixWrap<SC,LO,GO,NO> >
118 template<
typename SC,
typename LO,
typename GO,
typename NO>
119 RCP<Xpetra::Matrix<SC, LO, GO, NO> >
122 template<
typename SC,
typename LO,
typename GO,
typename NO>
123 RCP<Xpetra::MultiVector<SC, LO, GO, NO> >
127 #ifdef HAVE_MUELU_TPETRA 128 template<
typename SC,
typename LO,
typename GO,
typename NO>
129 RCP<Xpetra::Matrix<SC, LO, GO, NO> >
133 template<
typename SC,
typename LO,
typename GO,
typename NO>
134 RCP<Xpetra::MultiVector<SC, LO, GO, NO> >
145 template <
class Scalar,
146 class LocalOrdinal = int,
147 class GlobalOrdinal = LocalOrdinal,
148 class Node = KokkosClassic::DefaultNode::DefaultNodeType>
149 class Utilities :
public UtilitiesBase<Scalar, LocalOrdinal, GlobalOrdinal, Node> {
150 #undef MUELU_UTILITIES_SHORT 154 typedef typename Teuchos::ScalarTraits<Scalar>::magnitudeType
Magnitude;
156 #ifdef HAVE_MUELU_EPETRA 159 static RCP<const Epetra_MultiVector>
MV2EpetraMV(RCP<Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
const vec);
160 static RCP< Epetra_MultiVector>
MV2NonConstEpetraMV(RCP<Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> > vec);
162 static const Epetra_MultiVector&
MV2EpetraMV(
const Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& vec);
163 static Epetra_MultiVector&
MV2NonConstEpetraMV(Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& vec);
165 static RCP<const Epetra_CrsMatrix>
Op2EpetraCrs(RCP<
const Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > Op);
166 static RCP< Epetra_CrsMatrix>
Op2NonConstEpetraCrs(RCP<Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > Op);
168 static const Epetra_CrsMatrix&
Op2EpetraCrs(
const Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& Op);
169 static Epetra_CrsMatrix&
Op2NonConstEpetraCrs(Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& Op);
171 static const Epetra_Map&
Map2EpetraMap(
const Xpetra::Map<LocalOrdinal,GlobalOrdinal,Node>& map);
175 #ifdef HAVE_MUELU_TPETRA 176 static RCP<const Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
MV2TpetraMV(RCP<Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
const vec);
178 static RCP< Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
MV2NonConstTpetraMV(RCP<Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> > vec);
179 static RCP< Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
MV2NonConstTpetraMV2(Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& vec);
181 static const Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
MV2TpetraMV(
const Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& vec);
182 static Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
MV2NonConstTpetraMV(Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& vec);
184 static RCP<const Tpetra::CrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
Op2TpetraCrs(RCP<
const Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > Op);
185 static RCP< Tpetra::CrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
Op2NonConstTpetraCrs(RCP<Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > Op);
187 static const Tpetra::CrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
Op2TpetraCrs(
const Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& Op);
188 static Tpetra::CrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
Op2NonConstTpetraCrs(Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& Op);
190 static RCP<const Tpetra::RowMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
Op2TpetraRow(RCP<
const Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > Op);
191 static RCP< Tpetra::RowMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
Op2NonConstTpetraRow(RCP<Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > Op);
194 static const RCP<const Tpetra::Map<LocalOrdinal, GlobalOrdinal, Node> >
Map2TpetraMap(
const Xpetra::Map<LocalOrdinal,GlobalOrdinal,Node>& map);
199 static RCP<Xpetra::Vector<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
GetMatrixDiagonalInverse(
const Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& A, Magnitude tol = Teuchos::ScalarTraits<Scalar>::eps()*100) {
return MueLu::UtilitiesBase<Scalar,LocalOrdinal,GlobalOrdinal,Node>::GetMatrixDiagonalInverse(A,tol); }
202 static Teuchos::Array<Magnitude>
ResidualNorm(
const Xpetra::Operator<Scalar,LocalOrdinal,GlobalOrdinal,Node>& Op,
const Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& X,
const Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& RHS) {
return MueLu::UtilitiesBase<Scalar,LocalOrdinal,GlobalOrdinal,Node>::ResidualNorm(Op,X,RHS); }
203 static RCP<Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
Residual(
const Xpetra::Operator<Scalar,LocalOrdinal,GlobalOrdinal,Node>& Op,
const Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& X,
const Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& RHS) {
return MueLu::UtilitiesBase<Scalar,LocalOrdinal,GlobalOrdinal,Node>::Residual(Op,X,RHS); }
206 static typename Teuchos::ScalarTraits<Scalar>::magnitudeType
Distance2(
const Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& v, LocalOrdinal i0, LocalOrdinal i1) {
return MueLu::UtilitiesBase<Scalar,LocalOrdinal,GlobalOrdinal,Node>::Distance2(v,i0,i1); }
210 static Scalar
PowerMethod(
const Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& A,
bool scaleByDiag =
true,
211 LocalOrdinal niters = 10, Magnitude tolerance = 1e-2,
bool verbose =
false,
unsigned int seed = 123) {
215 static Scalar
Frobenius(
const Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& A,
const Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& B) {
219 static void MyOldScaleMatrix(Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& Op,
const Teuchos::ArrayRCP<const Scalar>& scalingVector,
bool doInverse =
true,
220 bool doFillComplete =
true,
bool doOptimizeStorage =
true);
222 static void MyOldScaleMatrix_Epetra(Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& Op,
const Teuchos::ArrayRCP<Scalar>& scalingVector,
223 bool doFillComplete,
bool doOptimizeStorage);
224 static void MyOldScaleMatrix_Tpetra(Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& Op,
const Teuchos::ArrayRCP<Scalar>& scalingVector,
225 bool doFillComplete,
bool doOptimizeStorage);
227 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());
232 #ifdef HAVE_MUELU_EPETRA 249 typedef typename Teuchos::ScalarTraits<Scalar>::magnitudeType
Magnitude;
252 typedef Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node>
CrsMatrixWrap;
253 typedef Xpetra::CrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>
CrsMatrix;
254 typedef Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>
Matrix;
255 typedef Xpetra::Vector<Scalar,LocalOrdinal,GlobalOrdinal,Node>
Vector;
256 typedef Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>
MultiVector;
257 typedef Xpetra::Map<LocalOrdinal,GlobalOrdinal,Node>
Map;
258 #ifdef HAVE_MUELU_EPETRA 259 typedef Xpetra::EpetraMapT<GlobalOrdinal,Node>
EpetraMap;
265 #ifdef HAVE_MUELU_EPETRA 268 static RCP<const Epetra_MultiVector>
MV2EpetraMV(RCP<MultiVector>
const vec) {
269 RCP<const EpetraMultiVector > tmpVec = rcp_dynamic_cast<EpetraMultiVector>(vec);
270 if (tmpVec == Teuchos::null)
271 throw Exceptions::BadCast(
"Cast from Xpetra::MultiVector to Xpetra::EpetraMultiVector failed");
272 return tmpVec->getEpetra_MultiVector();
275 RCP<const EpetraMultiVector> tmpVec = rcp_dynamic_cast<EpetraMultiVector>(vec);
276 if (tmpVec == Teuchos::null)
277 throw Exceptions::BadCast(
"Cast from Xpetra::MultiVector to Xpetra::EpetraMultiVector failed");
278 return tmpVec->getEpetra_MultiVector();
281 static const Epetra_MultiVector&
MV2EpetraMV(
const MultiVector& vec) {
282 const EpetraMultiVector& tmpVec =
dynamic_cast<const EpetraMultiVector&
>(vec);
283 return *(tmpVec.getEpetra_MultiVector());
286 const EpetraMultiVector& tmpVec =
dynamic_cast<const EpetraMultiVector&
>(vec);
287 return *(tmpVec.getEpetra_MultiVector());
290 static RCP<const Epetra_CrsMatrix>
Op2EpetraCrs(RCP<const Matrix> Op) {
291 RCP<const CrsMatrixWrap> crsOp = rcp_dynamic_cast<
const CrsMatrixWrap>(Op);
292 if (crsOp == Teuchos::null)
294 const RCP<const EpetraCrsMatrix>& tmp_ECrsMtx = rcp_dynamic_cast<
const EpetraCrsMatrix>(crsOp->getCrsMatrix());
295 if (tmp_ECrsMtx == Teuchos::null)
296 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::EpetraCrsMatrix failed");
297 return tmp_ECrsMtx->getEpetra_CrsMatrix();
300 RCP<const CrsMatrixWrap> crsOp = rcp_dynamic_cast<
const CrsMatrixWrap>(Op);
301 if (crsOp == Teuchos::null)
303 const RCP<const EpetraCrsMatrix> &tmp_ECrsMtx = rcp_dynamic_cast<
const EpetraCrsMatrix>(crsOp->getCrsMatrix());
304 if (tmp_ECrsMtx == Teuchos::null)
305 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::EpetraCrsMatrix failed");
306 return tmp_ECrsMtx->getEpetra_CrsMatrixNonConst();
311 const CrsMatrixWrap& crsOp =
dynamic_cast<const CrsMatrixWrap&
>(Op);
313 const EpetraCrsMatrix& tmp_ECrsMtx =
dynamic_cast<const EpetraCrsMatrix&
>(*crsOp.getCrsMatrix());
314 return *tmp_ECrsMtx.getEpetra_CrsMatrix();
315 }
catch (std::bad_cast) {
316 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::EpetraCrsMatrix failed");
318 }
catch (std::bad_cast) {
324 CrsMatrixWrap& crsOp =
dynamic_cast<CrsMatrixWrap&
>(Op);
326 EpetraCrsMatrix& tmp_ECrsMtx =
dynamic_cast<EpetraCrsMatrix&
>(*crsOp.getCrsMatrix());
327 return *tmp_ECrsMtx.getEpetra_CrsMatrixNonConst();
328 }
catch (std::bad_cast) {
329 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::EpetraCrsMatrix failed");
331 }
catch (std::bad_cast) {
337 RCP<const EpetraMap> xeMap = rcp_dynamic_cast<
const EpetraMap>(rcpFromRef(map));
338 if (xeMap == Teuchos::null)
339 throw Exceptions::BadCast(
"Utilities::Map2EpetraMap : Cast from Xpetra::Map to Xpetra::EpetraMap failed");
340 return xeMap->getEpetra_Map();
345 #ifdef HAVE_MUELU_TPETRA 346 static RCP<const Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
MV2TpetraMV(RCP<MultiVector>
const vec) {
348 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \ 349 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT)))) 352 RCP<const Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> > tmpVec = rcp_dynamic_cast<Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(vec);
353 if (tmpVec == Teuchos::null)
354 throw Exceptions::BadCast(
"Cast from Xpetra::MultiVector to Xpetra::TpetraMultiVector failed");
355 return tmpVec->getTpetra_MultiVector();
358 static RCP< Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
MV2NonConstTpetraMV(RCP<MultiVector> vec) {
359 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \ 360 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT)))) 363 RCP<const Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> > tmpVec = rcp_dynamic_cast<Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(vec);
364 if (tmpVec == Teuchos::null)
365 throw Exceptions::BadCast(
"Cast from Xpetra::MultiVector to Xpetra::TpetraMultiVector failed");
366 return tmpVec->getTpetra_MultiVector();
370 static RCP< Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
MV2NonConstTpetraMV2(MultiVector& vec) {
371 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \ 372 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT)))) 375 const Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& tmpVec =
dynamic_cast<const Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
>(vec);
376 return tmpVec.getTpetra_MultiVector();
380 static const Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
MV2TpetraMV(
const MultiVector& vec) {
381 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \ 382 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT)))) 385 const Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& tmpVec =
dynamic_cast<const Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
>(vec);
386 return *(tmpVec.getTpetra_MultiVector());
389 static Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
MV2NonConstTpetraMV(MultiVector& vec) {
390 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \ 391 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT)))) 394 const Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& tmpVec =
dynamic_cast<const Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
>(vec);
395 return *(tmpVec.getTpetra_MultiVector());
399 static RCP<const Tpetra::CrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
Op2TpetraCrs(RCP<const Matrix> Op) {
400 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \ 401 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT)))) 405 RCP<const CrsMatrixWrap> crsOp = rcp_dynamic_cast<
const CrsMatrixWrap>(Op);
406 if (crsOp == Teuchos::null)
408 const RCP<const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > &tmp_ECrsMtx = rcp_dynamic_cast<
const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(crsOp->getCrsMatrix());
409 if (tmp_ECrsMtx == Teuchos::null)
410 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix failed");
411 return tmp_ECrsMtx->getTpetra_CrsMatrix();
415 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \ 416 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT)))) 419 RCP<const CrsMatrixWrap> crsOp = rcp_dynamic_cast<
const CrsMatrixWrap>(Op);
420 if (crsOp == Teuchos::null)
422 const RCP<const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > &tmp_ECrsMtx = rcp_dynamic_cast<
const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(crsOp->getCrsMatrix());
423 if (tmp_ECrsMtx == Teuchos::null)
424 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix failed");
425 return tmp_ECrsMtx->getTpetra_CrsMatrixNonConst();
429 static const Tpetra::CrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
Op2TpetraCrs(
const Matrix& Op) {
430 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \ 431 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT)))) 435 const CrsMatrixWrap& crsOp =
dynamic_cast<const CrsMatrixWrap&
>(Op);
437 const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& tmp_ECrsMtx =
dynamic_cast<const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
>(*crsOp.getCrsMatrix());
438 return *tmp_ECrsMtx.getTpetra_CrsMatrix();
439 }
catch (std::bad_cast) {
440 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix failed");
442 }
catch (std::bad_cast) {
448 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \ 449 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT)))) 453 CrsMatrixWrap& crsOp =
dynamic_cast<CrsMatrixWrap&
>(Op);
455 Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& tmp_ECrsMtx =
dynamic_cast<Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
>(*crsOp.getCrsMatrix());
456 return *tmp_ECrsMtx.getTpetra_CrsMatrixNonConst();
457 }
catch (std::bad_cast) {
458 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix failed");
460 }
catch (std::bad_cast) {
466 static RCP<const Tpetra::RowMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
Op2TpetraRow(RCP<const Matrix> Op) {
467 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \ 468 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT)))) 471 RCP<const CrsMatrixWrap> crsOp = rcp_dynamic_cast<
const CrsMatrixWrap>(Op);
472 if (crsOp == Teuchos::null)
475 RCP<const CrsMatrix> crsMat = crsOp->getCrsMatrix();
476 const RCP<const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > tmp_Crs = rcp_dynamic_cast<
const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(crsMat);
477 RCP<const Xpetra::TpetraBlockCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > tmp_BlockCrs;
478 if(!tmp_Crs.is_null()) {
479 return tmp_Crs->getTpetra_CrsMatrixNonConst();
482 tmp_BlockCrs= rcp_dynamic_cast<
const Xpetra::TpetraBlockCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(crsMat);
483 if (tmp_BlockCrs.is_null())
484 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix and Xpetra::TpetraBlockCrsMatrix failed");
485 return tmp_BlockCrs->getTpetra_BlockCrsMatrixNonConst();
490 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \ 491 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT)))) 494 RCP<const CrsMatrixWrap> crsOp = rcp_dynamic_cast<
const CrsMatrixWrap>(Op);
495 if (crsOp == Teuchos::null)
498 RCP<const CrsMatrix> crsMat = crsOp->getCrsMatrix();
499 const RCP<const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > tmp_Crs = rcp_dynamic_cast<
const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(crsMat);
500 RCP<const Xpetra::TpetraBlockCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > tmp_BlockCrs;
501 if(!tmp_Crs.is_null()) {
502 return tmp_Crs->getTpetra_CrsMatrixNonConst();
505 tmp_BlockCrs= rcp_dynamic_cast<
const Xpetra::TpetraBlockCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(crsMat);
506 if (tmp_BlockCrs.is_null())
507 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix and Xpetra::TpetraBlockCrsMatrix failed");
508 return tmp_BlockCrs->getTpetra_BlockCrsMatrixNonConst();
514 static const RCP<const Tpetra::Map<LocalOrdinal,GlobalOrdinal,Node> >
Map2TpetraMap(
const Map& map) {
515 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \ 516 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT)))) 519 const RCP<const Xpetra::TpetraMap<LocalOrdinal,GlobalOrdinal,Node>>& tmp_TMap = rcp_dynamic_cast<
const Xpetra::TpetraMap<LocalOrdinal,GlobalOrdinal,Node> >(rcpFromRef(map));
520 if (tmp_TMap == Teuchos::null)
521 throw Exceptions::BadCast(
"Utilities::Map2TpetraMap : Cast from Xpetra::Map to Xpetra::TpetraMap failed");
522 return tmp_TMap->getTpetra_Map();
532 static Teuchos::Array<Magnitude>
ResidualNorm(
const Xpetra::Operator<Scalar,LocalOrdinal,GlobalOrdinal,Node>& Op,
const Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& X,
const Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& RHS) {
return MueLu::UtilitiesBase<Scalar,LocalOrdinal,GlobalOrdinal,Node>::ResidualNorm(Op,X,RHS); }
533 static RCP<MultiVector>
Residual(
const Xpetra::Operator<Scalar,LocalOrdinal,GlobalOrdinal,Node>& Op,
const Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& X,
const Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& RHS) {
return MueLu::UtilitiesBase<Scalar,LocalOrdinal,GlobalOrdinal,Node>::Residual(Op,X,RHS); }
540 static Scalar
PowerMethod(
const Matrix& A,
bool scaleByDiag =
true,
541 LocalOrdinal niters = 10, Magnitude tolerance = 1e-2,
bool verbose =
false,
unsigned int seed = 123) {
545 static Scalar
Frobenius(
const Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& A,
const Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& B) {
549 static void MyOldScaleMatrix(Matrix& Op,
const Teuchos::ArrayRCP<const Scalar>& scalingVector,
bool doInverse =
true,
550 bool doFillComplete =
true,
bool doOptimizeStorage =
true) {
551 Scalar one = Teuchos::ScalarTraits<Scalar>::one();
552 Teuchos::ArrayRCP<Scalar> sv(scalingVector.size());
554 for (
int i = 0; i < scalingVector.size(); ++i)
555 sv[i] = one / scalingVector[i];
557 for (
int i = 0; i < scalingVector.size(); ++i)
558 sv[i] = scalingVector[i];
561 switch (Op.getRowMap()->lib()) {
562 case Xpetra::UseTpetra:
566 case Xpetra::UseEpetra:
578 bool doFillComplete,
bool doOptimizeStorage) {
579 #ifdef HAVE_MUELU_TPETRA 580 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \ 581 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT)))) 582 throw Exceptions::RuntimeError(
"Matrix scaling is not possible because Tpetra has not been compiled with support for LO=GO=int.");
587 const RCP<const Tpetra::Map<LocalOrdinal,GlobalOrdinal,Node> > rowMap = tpOp.getRowMap();
588 const RCP<const Tpetra::Map<LocalOrdinal,GlobalOrdinal,Node> > domainMap = tpOp.getDomainMap();
589 const RCP<const Tpetra::Map<LocalOrdinal,GlobalOrdinal,Node> > rangeMap = tpOp.getRangeMap();
591 size_t maxRowSize = tpOp.getNodeMaxNumRowEntries();
592 if (maxRowSize == Teuchos::as<size_t>(-1))
595 std::vector<Scalar> scaledVals(maxRowSize);
596 if (tpOp.isFillComplete())
599 if (Op.isLocallyIndexed() ==
true) {
600 Teuchos::ArrayView<const LocalOrdinal> cols;
601 Teuchos::ArrayView<const Scalar> vals;
603 for (
size_t i = 0; i < rowMap->getNodeNumElements(); ++i) {
604 tpOp.getLocalRowView(i, cols, vals);
605 size_t nnz = tpOp.getNumEntriesInLocalRow(i);
606 if (nnz > maxRowSize) {
608 scaledVals.resize(maxRowSize);
610 for (
size_t j = 0; j < nnz; ++j)
611 scaledVals[j] = vals[j]*scalingVector[i];
614 Teuchos::ArrayView<const Scalar> valview(&scaledVals[0], nnz);
615 tpOp.replaceLocalValues(i, cols, valview);
620 Teuchos::ArrayView<const GlobalOrdinal> cols;
621 Teuchos::ArrayView<const Scalar> vals;
623 for (
size_t i = 0; i < rowMap->getNodeNumElements(); ++i) {
624 GlobalOrdinal gid = rowMap->getGlobalElement(i);
625 tpOp.getGlobalRowView(gid, cols, vals);
626 size_t nnz = tpOp.getNumEntriesInGlobalRow(gid);
627 if (nnz > maxRowSize) {
629 scaledVals.resize(maxRowSize);
632 for (
size_t j = 0; j < nnz; ++j)
633 scaledVals[j] = vals[j]*scalingVector[i];
636 Teuchos::ArrayView<const Scalar> valview(&scaledVals[0], nnz);
637 tpOp.replaceGlobalValues(gid, cols, valview);
642 if (doFillComplete) {
643 if (domainMap == Teuchos::null || rangeMap == Teuchos::null)
644 throw Exceptions::RuntimeError(
"In Utilities::Scaling: cannot fillComplete because the domain and/or range map hasn't been defined");
646 RCP<Teuchos::ParameterList> params = rcp(
new Teuchos::ParameterList());
647 params->set(
"Optimize Storage", doOptimizeStorage);
648 params->set(
"No Nonlocal Changes",
true);
649 Op.fillComplete(Op.getDomainMap(), Op.getRangeMap(), params);
660 static void MyOldScaleMatrix_Epetra (Matrix& Op,
const Teuchos::ArrayRCP<Scalar>& scalingVector,
bool doFillComplete,
bool doOptimizeStorage) {
661 #ifdef HAVE_MUELU_EPETRA 666 Epetra_Map
const &rowMap = epOp.RowMap();
671 for (
int i = 0; i < rowMap.NumMyElements(); ++i) {
672 epOp.ExtractMyRowView(i, nnz, vals, cols);
673 for (
int j = 0; j < nnz; ++j)
674 vals[j] *= scalingVector[i];
682 #endif // HAVE_MUELU_EPETRA 690 static RCP<Matrix>
Transpose(Matrix& Op,
bool optimizeTranspose =
false,
const std::string & label = std::string()) {
691 switch (Op.getRowMap()->lib()) {
692 case Xpetra::UseTpetra: {
693 #ifdef HAVE_MUELU_TPETRA 694 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \ 695 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT)))) 696 throw Exceptions::RuntimeError(
"Utilities::Transpose: Tpetra is not compiled with LO=GO=int. Add TPETRA_INST_INT_INT:BOOL=ON to your configuration!");
702 RCP<Tpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > A;
703 Tpetra::RowMatrixTransposer<Scalar, LocalOrdinal, GlobalOrdinal, Node> transposer(rcpFromRef(tpetraOp),label);
704 A = transposer.createTranspose();
705 RCP<Xpetra::TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > AA = rcp(
new Xpetra::TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(A));
706 RCP<CrsMatrix> AAA = rcp_implicit_cast<CrsMatrix>(AA);
711 catch (std::exception& e) {
712 std::cout <<
"threw exception '" << e.what() <<
"'" << std::endl;
721 case Xpetra::UseEpetra:
723 #if defined(HAVE_MUELU_EPETRA) && defined(HAVE_MUELU_EPETRAEXT) 724 Teuchos::TimeMonitor tm(*Teuchos::TimeMonitor::getNewTimer(
"ZZ Entire Transpose"));
728 Epetra_CrsMatrix * A =
dynamic_cast<Epetra_CrsMatrix*
>(&transposer(epetraOp));
731 RCP<Epetra_CrsMatrix> rcpA(A);
732 RCP<EpetraCrsMatrix> AA = rcp(
new EpetraCrsMatrix(rcpA));
733 RCP<CrsMatrix> AAA = rcp_implicit_cast<CrsMatrix>(AA);
735 AAAA->fillComplete(Op.getRangeMap(), Op.getDomainMap());
748 return Teuchos::null;
753 #endif // HAVE_MUELU_EPETRA 761 long ExtractNonSerializableData(
const Teuchos::ParameterList& inList, Teuchos::ParameterList& serialList, Teuchos::ParameterList& nonSerialList);
773 #ifdef HAVE_MUELU_EPETRA 778 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
779 RCP<Xpetra::Matrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
781 typedef Xpetra::EpetraCrsMatrixT<GlobalOrdinal, Node> XECrsMatrix;
782 typedef Xpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> XCrsMatrix;
783 typedef Xpetra::CrsMatrixWrap<Scalar, LocalOrdinal, GlobalOrdinal, Node> XCrsMatrixWrap;
785 RCP<XCrsMatrix> Atmp = rcp(
new XECrsMatrix(A));
786 return rcp(
new XCrsMatrixWrap(Atmp));
793 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
794 RCP<Xpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
796 return rcp(
new Xpetra::EpetraMultiVectorT<GlobalOrdinal, Node>(V));
800 #ifdef HAVE_MUELU_TPETRA 805 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
806 RCP<Xpetra::Matrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
808 typedef Xpetra::TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> XTCrsMatrix;
809 typedef Xpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> XCrsMatrix;
810 typedef Xpetra::CrsMatrixWrap<Scalar, LocalOrdinal, GlobalOrdinal, Node> XCrsMatrixWrap;
812 RCP<XCrsMatrix> Atmp = rcp(
new XTCrsMatrix(Atpetra));
813 return rcp(
new XCrsMatrixWrap(Atmp));
820 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
821 RCP<Xpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
823 return rcp(
new Xpetra::TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>(Vtpetra));
830 std::ostringstream buf;
835 #ifdef HAVE_MUELU_EPETRA 840 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
841 RCP<Xpetra::Matrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
848 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
849 RCP<Xpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
853 #ifdef HAVE_MUELU_TPETRA 858 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
859 RCP<Xpetra::Matrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
866 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
867 RCP<Xpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
873 #define MUELU_UTILITIES_SHORT 874 #endif // MUELU_UTILITIES_DECL_HPP static Teuchos::ArrayRCP< Scalar > GetMatrixDiagonal(const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A)
static void SetRandomSeed(const Teuchos::Comm< int > &comm)
static Teuchos::ScalarTraits< Scalar >::magnitudeType Distance2(const MultiVector &v, LocalOrdinal i0, LocalOrdinal i1)
Exception indicating invalid cast attempted.
Xpetra::EpetraMapT< GlobalOrdinal, Node > EpetraMap
static void SetRandomSeed(const Teuchos::Comm< int > &comm)
static const RCP< const Tpetra::Map< LocalOrdinal, GlobalOrdinal, Node > > Map2TpetraMap(const Xpetra::Map< LocalOrdinal, GlobalOrdinal, Node > &map)
static const Epetra_MultiVector & MV2EpetraMV(const MultiVector &vec)
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.
RCP< Xpetra::Matrix< SC, LO, GO, NO > > EpetraCrs_To_XpetraMatrix(const Teuchos::RCP< Epetra_CrsMatrix > &A)
std::string toString(const T &what)
Little helper function to convert non-string types to strings.
RCP< Xpetra::MultiVector< SC, LO, GO, NO > > TpetraCrs_To_XpetraMultiVector(const Teuchos::RCP< Tpetra::MultiVector< SC, LO, GO, NO > > &Vtpetra)
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)
Xpetra::Map< LocalOrdinal, GlobalOrdinal, Node > Map
static RCP< const Tpetra::CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op2TpetraCrs(RCP< const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op)
static const Epetra_Map & Map2EpetraMap(const Xpetra::Map< LocalOrdinal, GlobalOrdinal, Node > &map)
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())
static RCP< Matrix > Crs2Op(RCP< CrsMatrix > Op)
static RCP< Epetra_CrsMatrix > Op2NonConstEpetraCrs(RCP< Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op)
static const Tpetra::CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > & Op2TpetraCrs(const Matrix &Op)
static Scalar PowerMethod(const Matrix &A, bool scaleByDiag=true, LocalOrdinal niters=10, Magnitude tolerance=1e-2, bool verbose=false, unsigned int seed=123)
Simple transpose for Tpetra::CrsMatrix types.
static Scalar PowerMethod(const Matrix &A, bool scaleByDiag=true, LocalOrdinal niters=10, Magnitude tolerance=1e-2, bool verbose=false, unsigned int seed=123)
Xpetra::CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > CrsMatrix
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.
static Teuchos::ScalarTraits< Scalar >::magnitudeType Distance2(const Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &v, LocalOrdinal i0, LocalOrdinal i1)
static Tpetra::CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > & Op2NonConstTpetraCrs(Matrix &Op)
Teuchos::ScalarTraits< Scalar >::magnitudeType Magnitude
Xpetra::CrsMatrixWrap< Scalar, LocalOrdinal, GlobalOrdinal, Node > CrsMatrixWrap
static RCP< const Epetra_CrsMatrix > Op2EpetraCrs(RCP< const Matrix > Op)
Xpetra::Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > Vector
static Teuchos::ArrayRCP< Scalar > GetLumpedMatrixDiagonal(const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A)
Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > Matrix
static RCP< const Tpetra::RowMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op2TpetraRow(RCP< const Matrix > Op)
static Teuchos::ScalarTraits< Scalar >::magnitudeType Distance2(const Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &v, LocalOrdinal i0, LocalOrdinal i1)
Squared distance between two rows in a multivector.
static RCP< const Epetra_MultiVector > MV2EpetraMV(RCP< MultiVector > const vec)
Helper utility to pull out the underlying Epetra objects from an Xpetra object.
Namespace for MueLu classes and methods.
static const RCP< const Tpetra::Map< LocalOrdinal, GlobalOrdinal, Node > > Map2TpetraMap(const Map &map)
Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > MultiVector
bool IsParamMuemexVariable(const std::string &name)
static RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > MV2NonConstTpetraMV(RCP< Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > vec)
static void PauseForDebugger()
static RCP< Xpetra::Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > GetMatrixDiagonalInverse(const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, Magnitude tol=Teuchos::ScalarTraits< Scalar >::eps()*100)
static void MyOldScaleMatrix_Tpetra(Matrix &Op, const Teuchos::ArrayRCP< Scalar > &scalingVector, bool doFillComplete, bool doOptimizeStorage)
static Scalar Frobenius(const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B)
static RCP< Vector > GetMatrixDiagonalInverse(const Matrix &A, Magnitude tol=Teuchos::ScalarTraits< Scalar >::eps()*100)
static RCP< const Tpetra::CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op2TpetraCrs(RCP< const Matrix > Op)
static Epetra_CrsMatrix & Op2NonConstEpetraCrs(Matrix &Op)
static RCP< Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Residual(const Xpetra::Operator< Scalar, LocalOrdinal, GlobalOrdinal, Node > &Op, const Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &X, const Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &RHS)
static RCP< Tpetra::RowMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op2NonConstTpetraRow(RCP< Matrix > Op)
static RCP< Epetra_MultiVector > MV2NonConstEpetraMV(RCP< Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > vec)
static RCP< Matrix > Transpose(Matrix &Op, bool optimizeTranspose=false, const std::string &label=std::string())
Transpose a Xpetra::Matrix.
static Epetra_MultiVector & MV2NonConstEpetraMV(MultiVector &vec)
static RCP< Epetra_MultiVector > MV2NonConstEpetraMV(RCP< MultiVector > vec)
static RCP< Tpetra::CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op2NonConstTpetraCrs(RCP< Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op)
static Teuchos::ArrayRCP< Scalar > GetMatrixDiagonal(const Matrix &A)
static RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > MV2NonConstTpetraMV(RCP< MultiVector > vec)
static Teuchos::ArrayRCP< Scalar > GetLumpedMatrixDiagonal(const Matrix &A)
static RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > MV2NonConstTpetraMV2(Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &vec)
static const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > & MV2TpetraMV(const MultiVector &vec)
static Teuchos::ArrayRCP< const bool > DetectDirichletRows(const Matrix &A, const Magnitude &tol=Teuchos::ScalarTraits< Scalar >::zero())
static RCP< const Tpetra::RowMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op2TpetraRow(RCP< const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op)
Xpetra::EpetraCrsMatrixT< GlobalOrdinal, Node > EpetraCrsMatrix
static RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > MV2NonConstTpetraMV2(MultiVector &vec)
static void PauseForDebugger()
static Teuchos::ArrayRCP< Scalar > GetLumpedMatrixDiagonal(const Matrix &A)
Extract Matrix Diagonal of lumped matrix.
static Scalar Frobenius(const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B)
Frobenius inner product of two matrices.
RCP< Xpetra::CrsMatrixWrap< SC, LO, GO, NO > > Convert_Epetra_CrsMatrix_ToXpetra_CrsMatrixWrap(RCP< Epetra_CrsMatrix > &epAB)
static void PauseForDebugger()
RCP< Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > TpetraMultiVector_To_XpetraMultiVector(const Teuchos::RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &Vtpetra)
static RCP< Teuchos::FancyOStream > MakeFancy(std::ostream &os)
static Teuchos::Array< Magnitude > ResidualNorm(const Xpetra::Operator< Scalar, LocalOrdinal, GlobalOrdinal, Node > &Op, const MultiVector &X, const MultiVector &RHS)
static RCP< Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Crs2Op(RCP< Xpetra::CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op)
static void MyOldScaleMatrix_Epetra(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)
static RCP< MultiVector > Residual(const Xpetra::Operator< Scalar, LocalOrdinal, GlobalOrdinal, Node > &Op, const MultiVector &X, const MultiVector &RHS)
static RCP< Xpetra::Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > GetMatrixOverlappedDiagonal(const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A)
static Teuchos::Array< Magnitude > ResidualNorm(const Xpetra::Operator< Scalar, LocalOrdinal, GlobalOrdinal, Node > &Op, const Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &X, const Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &RHS)
static RCP< Matrix > Crs2Op(RCP< CrsMatrix > Op)
Xpetra::CrsMatrixWrap< Scalar, LocalOrdinal, GlobalOrdinal, Node > CrsMatrixWrap
RCP< Xpetra::Matrix< SC, LO, GO, NO > > TpetraCrs_To_XpetraMatrix(const Teuchos::RCP< Tpetra::CrsMatrix< SC, LO, GO, NO > > &Atpetra)
static RCP< Vector > GetMatrixOverlappedDiagonal(const Matrix &A)
static RCP< Teuchos::FancyOStream > MakeFancy(std::ostream &os)
static Teuchos::ArrayRCP< const bool > DetectDirichletRows(const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Magnitude &tol=Teuchos::ScalarTraits< Scalar >::zero())
Detect Dirichlet rows.
static RCP< Tpetra::CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op2NonConstTpetraCrs(RCP< Matrix > Op)
static Teuchos::ArrayRCP< const bool > DetectDirichletRows(const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Magnitude &tol=Teuchos::ScalarTraits< Scalar >::magnitude(0.))
Exception throws to report errors in the internal logical of the program.
static Teuchos::Array< Magnitude > ResidualNorm(const Xpetra::Operator< Scalar, LocalOrdinal, GlobalOrdinal, Node > &Op, const Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &X, const Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &RHS)
Kokkos::Compat::KokkosSerialWrapperNode EpetraNode
static Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > & MV2NonConstTpetraMV(MultiVector &vec)
static const Epetra_Map & Map2EpetraMap(const Map &map)
static Scalar Frobenius(const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B)
void TokenizeStringAndStripWhiteSpace(const std::string &stream, std::vector< std::string > &tokenList, const char *delimChars)
static void MyOldScaleMatrix(Matrix &Op, const Teuchos::ArrayRCP< const Scalar > &scalingVector, bool doInverse=true, bool doFillComplete=true, bool doOptimizeStorage=true)
static RCP< Tpetra::RowMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op2NonConstTpetraRow(RCP< Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op)
static RCP< Vector > GetMatrixOverlappedDiagonal(const Matrix &A)
Extract Overlapped Matrix Diagonal.
static RCP< Teuchos::FancyOStream > MakeFancy(std::ostream &os)
static void MyOldScaleMatrix_Tpetra(Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &Op, const Teuchos::ArrayRCP< Scalar > &scalingVector, bool doFillComplete, bool doOptimizeStorage)
RCP< Xpetra::MultiVector< SC, LO, GO, NO > > EpetraMultiVector_To_XpetraMultiVector(const Teuchos::RCP< Epetra_MultiVector > &V)
static const Epetra_CrsMatrix & Op2EpetraCrs(const Matrix &Op)
static RCP< MultiVector > Residual(const Xpetra::Operator< Scalar, LocalOrdinal, GlobalOrdinal, Node > &Op, const Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &X, const Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &RHS)
static Scalar PowerMethod(const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, bool scaleByDiag=true, LocalOrdinal niters=10, Magnitude tolerance=1e-2, bool verbose=false, unsigned int seed=123)
static RCP< Epetra_CrsMatrix > Op2NonConstEpetraCrs(RCP< Matrix > Op)
Teuchos::ScalarTraits< Scalar >::magnitudeType Magnitude
static RCP< Vector > GetMatrixDiagonalInverse(const Matrix &A, Magnitude tol=Teuchos::ScalarTraits< Scalar >::eps()*100)
Extract Matrix Diagonal.
static void SetRandomSeed(const Teuchos::Comm< int > &comm)
Set seed for random number generator.
static Teuchos::ArrayRCP< Scalar > GetMatrixDiagonal(const Matrix &A)
Extract Matrix Diagonal.
static void MyOldScaleMatrix_Epetra(Matrix &Op, const Teuchos::ArrayRCP< Scalar > &scalingVector, bool doFillComplete, bool doOptimizeStorage)
long ExtractNonSerializableData(const Teuchos::ParameterList &inList, Teuchos::ParameterList &serialList, Teuchos::ParameterList &nonSerialList)
Xpetra::EpetraMultiVectorT< GlobalOrdinal, Node > EpetraMultiVector