46 #ifndef MUELU_UTILITIES_KOKKOS_DECL_HPP 47 #define MUELU_UTILITIES_KOKKOS_DECL_HPP 50 #if defined(HAVE_MUELU_KOKKOS_REFACTOR) 54 #include <Teuchos_ScalarTraits.hpp> 55 #include <Teuchos_ParameterList.hpp> 57 #include <Xpetra_BlockedCrsMatrix_fwd.hpp> 58 #include <Xpetra_CrsMatrix_fwd.hpp> 59 #include <Xpetra_CrsMatrixWrap_fwd.hpp> 60 #include <Xpetra_ExportFactory.hpp> 61 #include <Xpetra_ImportFactory_fwd.hpp> 62 #include <Xpetra_MapFactory_fwd.hpp> 63 #include <Xpetra_Map_fwd.hpp> 64 #include <Xpetra_MatrixFactory_fwd.hpp> 65 #include <Xpetra_Matrix_fwd.hpp> 66 #include <Xpetra_MultiVectorFactory_fwd.hpp> 67 #include <Xpetra_MultiVector_fwd.hpp> 68 #include <Xpetra_Operator_fwd.hpp> 69 #include <Xpetra_VectorFactory_fwd.hpp> 70 #include <Xpetra_Vector_fwd.hpp> 72 #include <Xpetra_IO.hpp> 74 #ifdef HAVE_MUELU_EPETRA 75 #include <Epetra_MultiVector.h> 76 #include <Epetra_CrsMatrix.h> 77 #include <Xpetra_EpetraCrsMatrix_fwd.hpp> 78 #include <Xpetra_EpetraMultiVector_fwd.hpp> 82 #include "MueLu_Utilities.hpp" 83 #include "MueLu_UtilitiesBase.hpp" 85 #ifdef HAVE_MUELU_TPETRA 86 #include <Tpetra_CrsMatrix.hpp> 87 #include <Tpetra_Map.hpp> 88 #include <Tpetra_MultiVector.hpp> 89 #include <Xpetra_TpetraCrsMatrix_fwd.hpp> 90 #include <Xpetra_TpetraMultiVector_fwd.hpp> 103 template <
class Scalar,
104 class LocalOrdinal = int,
105 class GlobalOrdinal = LocalOrdinal,
106 class Node = KokkosClassic::DefaultNode::DefaultNodeType>
108 #undef MUELU_UTILITIES_KOKKOS_SHORT 112 typedef typename Teuchos::ScalarTraits<SC>::magnitudeType Magnitude;
114 #ifdef HAVE_MUELU_EPETRA 117 static RCP<const Epetra_MultiVector> MV2EpetraMV(RCP<MultiVector>
const vec) {
return Utilities::MV2EpetraMV(vec); }
120 static const Epetra_MultiVector& MV2EpetraMV(
const MultiVector& vec) {
return Utilities::MV2EpetraMV(vec); }
133 #ifdef HAVE_MUELU_TPETRA 136 static RCP<const Tpetra::MultiVector<SC,LO,GO,NO> > MV2TpetraMV(RCP<MultiVector>
const vec) {
return Utilities::MV2TpetraMV(vec); }
140 static const Tpetra::MultiVector<SC,LO,GO,NO>& MV2TpetraMV(
const MultiVector& vec) {
return Utilities::MV2TpetraMV(vec); }
143 static RCP<const Tpetra::CrsMatrix<SC,LO,GO,NO> > Op2TpetraCrs(RCP<const Matrix> Op) {
return Utilities::Op2TpetraCrs(Op); }
146 static const Tpetra::CrsMatrix<SC,LO,GO,NO>& Op2TpetraCrs(
const Matrix& Op) {
return Utilities::Op2TpetraCrs(Op); }
149 static RCP<const Tpetra::RowMatrix<SC,LO,GO,NO> > Op2TpetraRow(RCP<const Matrix> Op) {
return Utilities::Op2TpetraRow(Op); }
152 static const RCP<const Tpetra::Map<LO, GO, NO> > Map2TpetraMap(
const Map& map) {
return Utilities::Map2TpetraMap(map); }
155 static RCP<Xpetra::Matrix<SC,LO,GO,NO> > Crs2Op(RCP<CrsMatrix> Op) {
return Utilities::Crs2Op(Op); }
163 static Teuchos::ArrayRCP<SC> GetMatrixDiagonal(
const Matrix& A);
171 static RCP<Vector> GetMatrixDiagonalInverse(
const Matrix& A, Magnitude tol = Teuchos::ScalarTraits<SC>::eps()*100);
181 static Teuchos::ArrayRCP<SC> GetLumpedMatrixDiagonal(
const Matrix& A);
190 static RCP<Vector> GetMatrixOverlappedDiagonal(
const Matrix& A);
196 static Teuchos::Array<Magnitude> ResidualNorm(
const Operator& Op,
const MultiVector& X,
const MultiVector& RHS) {
200 static RCP<MultiVector> Residual(
const Operator& Op,
const MultiVector& X,
const MultiVector& RHS) {
204 static void PauseForDebugger();
221 static SC PowerMethod(
const Matrix& A,
bool scaleByDiag =
true,
222 LO niters = 10, Magnitude tolerance = 1e-2,
bool verbose =
false,
unsigned int seed = 123) {
226 static void MyOldScaleMatrix(Matrix& Op,
const Teuchos::ArrayRCP<const SC>& scalingVector,
bool doInverse =
true,
227 bool doFillComplete =
true,
bool doOptimizeStorage =
true);
229 static void MyOldScaleMatrix_Tpetra(Matrix& Op,
const Teuchos::ArrayRCP<SC>& scalingVector,
230 bool doFillComplete,
bool doOptimizeStorage);
232 static void MyOldScaleMatrix_Epetra(Matrix& Op,
const Teuchos::ArrayRCP<SC>& scalingVector,
233 bool doFillComplete,
bool doOptimizeStorage);
235 static RCP<Teuchos::FancyOStream> MakeFancy(std::ostream& os) {
return Utilities::MakeFancy(os); }
241 static typename Teuchos::ScalarTraits<Scalar>::magnitudeType Distance2(
const MultiVector& v, LocalOrdinal i0, LocalOrdinal i1) {
return Utilities::Distance2(v, i0, i1); }
268 static RCP<Matrix> Transpose(Matrix& Op,
bool optimizeTranspose =
false,
const std::string & label = std::string()) {
284 template <
class Node>
285 class Utilities_kokkos<double,int,int,Node> :
public UtilitiesBase<double,int,int,Node> {
287 typedef double Scalar;
288 typedef int LocalOrdinal;
289 typedef int GlobalOrdinal;
290 typedef typename Teuchos::ScalarTraits<Scalar>::magnitudeType Magnitude;
292 #ifdef HAVE_MUELU_EPETRA 295 static RCP<const Epetra_MultiVector> MV2EpetraMV(RCP<Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
const vec) {
296 RCP<const Xpetra::EpetraMultiVectorT<GlobalOrdinal,Node> > tmpVec = rcp_dynamic_cast<Xpetra::EpetraMultiVectorT<GlobalOrdinal,Node> >(vec);
297 if (tmpVec == Teuchos::null)
298 throw Exceptions::BadCast(
"Cast from Xpetra::MultiVector to Xpetra::EpetraMultiVector failed");
299 return tmpVec->getEpetra_MultiVector();
301 static RCP< Epetra_MultiVector> MV2NonConstEpetraMV(RCP<Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> > vec) {
302 RCP<const Xpetra::EpetraMultiVectorT<GlobalOrdinal,Node> > tmpVec = rcp_dynamic_cast<Xpetra::EpetraMultiVectorT<GlobalOrdinal,Node> >(vec);
303 if (tmpVec == Teuchos::null)
304 throw Exceptions::BadCast(
"Cast from Xpetra::MultiVector to Xpetra::EpetraMultiVector failed");
305 return tmpVec->getEpetra_MultiVector();
308 static const Epetra_MultiVector& MV2EpetraMV(
const Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> & vec) {
309 const Xpetra::EpetraMultiVectorT<GlobalOrdinal,Node> & tmpVec =
dynamic_cast<const Xpetra::EpetraMultiVectorT<GlobalOrdinal,Node> &
>(vec);
310 return *(tmpVec.getEpetra_MultiVector());
312 static Epetra_MultiVector& MV2NonConstEpetraMV(Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> & vec) {
313 const Xpetra::EpetraMultiVectorT<GlobalOrdinal,Node> & tmpVec =
dynamic_cast<const Xpetra::EpetraMultiVectorT<GlobalOrdinal,Node> &
>(vec);
314 return *(tmpVec.getEpetra_MultiVector());
317 static RCP<const Epetra_CrsMatrix> Op2EpetraCrs(RCP<
const Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > Op) {
318 RCP<const Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node> > crsOp = rcp_dynamic_cast<
const Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(Op);
319 if (crsOp == Teuchos::null)
320 throw Exceptions::BadCast(
"Cast from Xpetra::Matrix to Xpetra::CrsMatrixWrap failed");
321 const RCP<const Xpetra::EpetraCrsMatrixT<GlobalOrdinal,Node> >& tmp_ECrsMtx = rcp_dynamic_cast<
const Xpetra::EpetraCrsMatrixT<GlobalOrdinal,Node> >(crsOp->getCrsMatrix());
322 if (tmp_ECrsMtx == Teuchos::null)
323 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::EpetraCrsMatrix failed");
324 return tmp_ECrsMtx->getEpetra_CrsMatrix();
326 static RCP< Epetra_CrsMatrix> Op2NonConstEpetraCrs(RCP<Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > Op) {
327 RCP<const Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node> > crsOp = rcp_dynamic_cast<
const Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(Op);
328 if (crsOp == Teuchos::null)
329 throw Exceptions::BadCast(
"Cast from Xpetra::Matrix to Xpetra::CrsMatrixWrap failed");
330 const RCP<const Xpetra::EpetraCrsMatrixT<GlobalOrdinal,Node> > &tmp_ECrsMtx = rcp_dynamic_cast<
const Xpetra::EpetraCrsMatrixT<GlobalOrdinal,Node> >(crsOp->getCrsMatrix());
331 if (tmp_ECrsMtx == Teuchos::null)
332 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::EpetraCrsMatrix failed");
333 return tmp_ECrsMtx->getEpetra_CrsMatrixNonConst();
336 static const Epetra_CrsMatrix& Op2EpetraCrs(
const Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> & Op) {
338 const Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node>& crsOp =
dynamic_cast<const Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
>(Op);
340 const Xpetra::EpetraCrsMatrixT<GlobalOrdinal,Node> & tmp_ECrsMtx =
dynamic_cast<const Xpetra::EpetraCrsMatrixT<GlobalOrdinal,Node> &
>(*crsOp.getCrsMatrix());
341 return *tmp_ECrsMtx.getEpetra_CrsMatrix();
342 }
catch (std::bad_cast) {
343 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::EpetraCrsMatrix failed");
345 }
catch (std::bad_cast) {
346 throw Exceptions::BadCast(
"Cast from Xpetra::Matrix to Xpetra::CrsMatrixWrap failed");
349 static Epetra_CrsMatrix& Op2NonConstEpetraCrs(Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> & Op) {
351 Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node>& crsOp =
dynamic_cast<Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
>(Op);
353 Xpetra::EpetraCrsMatrixT<GlobalOrdinal,Node> & tmp_ECrsMtx =
dynamic_cast<Xpetra::EpetraCrsMatrixT<GlobalOrdinal,Node> &
>(*crsOp.getCrsMatrix());
354 return *tmp_ECrsMtx.getEpetra_CrsMatrixNonConst();
355 }
catch (std::bad_cast) {
356 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::EpetraCrsMatrix failed");
358 }
catch (std::bad_cast) {
359 throw Exceptions::BadCast(
"Cast from Xpetra::Matrix to Xpetra::CrsMatrixWrap failed");
363 static const Epetra_Map& Map2EpetraMap(
const Xpetra::Map<LocalOrdinal,GlobalOrdinal,Node> & map) {
364 RCP<const Xpetra::EpetraMapT<GlobalOrdinal,Node> > xeMap = rcp_dynamic_cast<
const Xpetra::EpetraMapT<GlobalOrdinal,Node> >(rcpFromRef(map));
365 if (xeMap == Teuchos::null)
366 throw Exceptions::BadCast(
"Utilities::Map2EpetraMap : Cast from Xpetra::Map to Xpetra::EpetraMap failed");
367 return xeMap->getEpetra_Map();
372 #ifdef HAVE_MUELU_TPETRA 375 static RCP<const Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> > MV2TpetraMV(RCP<Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
const vec) {
376 #ifdef HAVE_MUELU_TPETRA_INST_INT_INT 377 RCP<const Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> > tmpVec = rcp_dynamic_cast<Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(vec);
378 if (tmpVec == Teuchos::null)
379 throw Exceptions::BadCast(
"Cast from Xpetra::MultiVector to Xpetra::TpetraMultiVector failed");
380 return tmpVec->getTpetra_MultiVector();
382 throw Exceptions::RuntimeError(
"MV2TpetraMV: Tpetra has not been compiled with support for LO=GO=int.");
385 static RCP< Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> > MV2NonConstTpetraMV(RCP<Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> > vec) {
386 #ifdef HAVE_MUELU_TPETRA_INST_INT_INT 387 RCP<const Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> > tmpVec = rcp_dynamic_cast<Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(vec);
388 if (tmpVec == Teuchos::null)
389 throw Exceptions::BadCast(
"Cast from Xpetra::MultiVector to Xpetra::TpetraMultiVector failed");
390 return tmpVec->getTpetra_MultiVector();
392 throw Exceptions::RuntimeError(
"MV2NonConstTpetraMV: Tpetra has not been compiled with support for LO=GO=int.");
396 static RCP< Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> > MV2NonConstTpetraMV2(Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> & vec) {
397 #ifdef HAVE_MUELU_TPETRA_INST_INT_INT 398 const Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& tmpVec =
dynamic_cast<const Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
>(vec);
399 return tmpVec.getTpetra_MultiVector();
401 throw Exceptions::RuntimeError(
"MV2NonConstTpetraMV2: Tpetra has not been compiled with support for LO=GO=int.");
405 static const Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& MV2TpetraMV(
const Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> & vec) {
406 #ifdef HAVE_MUELU_TPETRA_INST_INT_INT 407 const Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& tmpVec =
dynamic_cast<const Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
>(vec);
408 return *(tmpVec.getTpetra_MultiVector());
410 throw Exceptions::RuntimeError(
"MV2TpetraMV: Tpetra has not been compiled with support for LO=GO=int.");
413 static Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& MV2NonConstTpetraMV(Xpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> & vec) {
414 #ifdef HAVE_MUELU_TPETRA_INST_INT_INT 415 const Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& tmpVec =
dynamic_cast<const Xpetra::TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
>(vec);
416 return *(tmpVec.getTpetra_MultiVector());
418 throw Exceptions::RuntimeError(
"MV2NonConstTpetraMV: Tpetra has not been compiled with support for LO=GO=int.");
422 static RCP<const Tpetra::CrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > Op2TpetraCrs(RCP<
const Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > Op) {
423 #ifdef HAVE_MUELU_TPETRA_INST_INT_INT 425 RCP<const Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node> > crsOp = rcp_dynamic_cast<
const Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(Op);
426 if (crsOp == Teuchos::null)
427 throw Exceptions::BadCast(
"Cast from Xpetra::Matrix to Xpetra::CrsMatrixWrap failed");
428 const RCP<const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > &tmp_ECrsMtx = rcp_dynamic_cast<
const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(crsOp->getCrsMatrix());
429 if (tmp_ECrsMtx == Teuchos::null)
430 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix failed");
431 return tmp_ECrsMtx->getTpetra_CrsMatrix();
433 throw Exceptions::RuntimeError(
"Op2TpetraCrs: Tpetra has not been compiled with support for LO=GO=int.");
436 static RCP< Tpetra::CrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > Op2NonConstTpetraCrs(RCP<Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > Op){
437 #ifdef HAVE_MUELU_TPETRA_INST_INT_INT 438 RCP<const Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node> > crsOp = rcp_dynamic_cast<
const Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(Op);
439 if (crsOp == Teuchos::null)
440 throw Exceptions::BadCast(
"Cast from Xpetra::Matrix to Xpetra::CrsMatrixWrap failed");
441 const RCP<const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > &tmp_ECrsMtx = rcp_dynamic_cast<
const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(crsOp->getCrsMatrix());
442 if (tmp_ECrsMtx == Teuchos::null)
443 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix failed");
444 return tmp_ECrsMtx->getTpetra_CrsMatrixNonConst();
446 throw Exceptions::RuntimeError(
"Op2NonConstTpetraCrs: Tpetra has not been compiled with support for LO=GO=int.");
450 static const Tpetra::CrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& Op2TpetraCrs(
const Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> & Op) {
451 #ifdef HAVE_MUELU_TPETRA_INST_INT_INT 453 const Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node>& crsOp =
dynamic_cast<const Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
>(Op);
455 const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& tmp_ECrsMtx =
dynamic_cast<const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
>(*crsOp.getCrsMatrix());
456 return *tmp_ECrsMtx.getTpetra_CrsMatrix();
457 }
catch (std::bad_cast) {
458 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix failed");
460 }
catch (std::bad_cast) {
461 throw Exceptions::BadCast(
"Cast from Xpetra::Matrix to Xpetra::CrsMatrixWrap failed");
464 throw Exceptions::RuntimeError(
"Op2TpetraCrs: Tpetra has not been compiled with support for LO=GO=int.");
467 static Tpetra::CrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& Op2NonConstTpetraCrs(Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> & Op) {
468 #ifdef HAVE_MUELU_TPETRA_INST_INT_INT 470 Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node>& crsOp =
dynamic_cast<Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
>(Op);
472 Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& tmp_ECrsMtx =
dynamic_cast<Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>&
>(*crsOp.getCrsMatrix());
473 return *tmp_ECrsMtx.getTpetra_CrsMatrixNonConst();
474 }
catch (std::bad_cast) {
475 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix failed");
477 }
catch (std::bad_cast) {
478 throw Exceptions::BadCast(
"Cast from Xpetra::Matrix to Xpetra::CrsMatrixWrap failed");
481 throw Exceptions::RuntimeError(
"Op2NonConstTpetraCrs: Tpetra has not been compiled with support for LO=GO=int.");
485 static RCP<const Tpetra::RowMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > Op2TpetraRow(RCP<
const Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > Op) {
486 #ifdef HAVE_MUELU_TPETRA_INST_INT_INT 487 RCP<const Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node> > crsOp = rcp_dynamic_cast<
const Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(Op);
488 if (crsOp == Teuchos::null)
489 throw Exceptions::BadCast(
"Cast from Xpetra::Matrix to Xpetra::CrsMatrixWrap failed");
491 RCP<const Xpetra::CrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > crsMat = crsOp->getCrsMatrix();
492 const RCP<const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > tmp_Crs = rcp_dynamic_cast<
const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(crsMat);
493 RCP<const Xpetra::TpetraBlockCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > tmp_BlockCrs;
494 if(!tmp_Crs.is_null()) {
495 return tmp_Crs->getTpetra_CrsMatrixNonConst();
498 tmp_BlockCrs= rcp_dynamic_cast<
const Xpetra::TpetraBlockCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(crsMat);
499 if (tmp_BlockCrs.is_null())
500 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix and Xpetra::TpetraBlockCrsMatrix failed");
501 return tmp_BlockCrs->getTpetra_BlockCrsMatrixNonConst();
504 throw Exceptions::RuntimeError(
"Op2TpetraRow: Tpetra has not been compiled with support for LO=GO=int.");
507 static RCP< Tpetra::RowMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > Op2NonConstTpetraRow(RCP<Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > Op) {
508 #ifdef HAVE_MUELU_TPETRA_INST_INT_INT 509 RCP<const Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node> > crsOp = rcp_dynamic_cast<
const Xpetra::CrsMatrixWrap<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(Op);
510 if (crsOp == Teuchos::null)
511 throw Exceptions::BadCast(
"Cast from Xpetra::Matrix to Xpetra::CrsMatrixWrap failed");
513 RCP<const Xpetra::CrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > crsMat = crsOp->getCrsMatrix();
514 const RCP<const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > tmp_Crs = rcp_dynamic_cast<
const Xpetra::TpetraCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(crsMat);
515 RCP<const Xpetra::TpetraBlockCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> > tmp_BlockCrs;
516 if(!tmp_Crs.is_null()) {
517 return tmp_Crs->getTpetra_CrsMatrixNonConst();
520 tmp_BlockCrs= rcp_dynamic_cast<
const Xpetra::TpetraBlockCrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >(crsMat);
521 if (tmp_BlockCrs.is_null())
522 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix and Xpetra::TpetraBlockCrsMatrix failed");
523 return tmp_BlockCrs->getTpetra_BlockCrsMatrixNonConst();
526 throw Exceptions::RuntimeError(
"Op2NonConstTpetraRow: Tpetra has not been compiled with support for LO=GO=int.");
531 static const RCP<const Tpetra::Map<LocalOrdinal,GlobalOrdinal,Node> > Map2TpetraMap(
const Xpetra::Map<LocalOrdinal,GlobalOrdinal,Node> & map) {
532 #ifdef HAVE_MUELU_TPETRA_INST_INT_INT 533 const RCP<const Xpetra::TpetraMap<LocalOrdinal,GlobalOrdinal,Node>>& tmp_TMap = rcp_dynamic_cast<
const Xpetra::TpetraMap<LocalOrdinal,GlobalOrdinal,Node> >(rcpFromRef(map));
534 if (tmp_TMap == Teuchos::null)
535 throw Exceptions::BadCast(
"Utilities::Map2TpetraMap : Cast from Xpetra::Map to Xpetra::TpetraMap failed");
536 return tmp_TMap->getTpetra_Map();
538 throw Exceptions::RuntimeError(
"Map2TpetraMap: Tpetra has not been compiled with support for LO=GO=int.");
545 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); }
548 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); }
549 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); }
552 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); }
558 static Scalar PowerMethod(
const Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& A,
bool scaleByDiag =
true,
559 LocalOrdinal niters = 10, Magnitude tolerance = 1e-2,
bool verbose =
false,
unsigned int seed = 123) {
563 static void MyOldScaleMatrix(Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& Op,
const Teuchos::ArrayRCP<const Scalar>& scalingVector,
bool doInverse =
true,
564 bool doFillComplete =
true,
bool doOptimizeStorage =
true) {
565 Scalar one = Teuchos::ScalarTraits<Scalar>::one();
566 Teuchos::ArrayRCP<Scalar> sv(scalingVector.size());
568 for (
int i = 0; i < scalingVector.size(); ++i)
569 sv[i] = one / scalingVector[i];
571 for (
int i = 0; i < scalingVector.size(); ++i)
572 sv[i] = scalingVector[i];
575 switch (Op.getRowMap()->lib()) {
576 case Xpetra::UseTpetra:
577 MyOldScaleMatrix_Tpetra(Op, sv, doFillComplete, doOptimizeStorage);
580 case Xpetra::UseEpetra:
581 MyOldScaleMatrix_Epetra(Op, sv, doFillComplete, doOptimizeStorage);
585 throw Exceptions::RuntimeError(
"Only Epetra and Tpetra matrices can be scaled.");
591 static void MyOldScaleMatrix_Tpetra(Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& Op,
const Teuchos::ArrayRCP<Scalar>& scalingVector,
592 bool doFillComplete,
bool doOptimizeStorage) {
593 #ifdef HAVE_MUELU_TPETRA 594 #ifdef HAVE_MUELU_TPETRA_INST_INT_INT 596 Tpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>& tpOp = Op2NonConstTpetraCrs(Op);
598 const RCP<const Tpetra::Map<LocalOrdinal,GlobalOrdinal,Node> > rowMap = tpOp.getRowMap();
599 const RCP<const Tpetra::Map<LocalOrdinal,GlobalOrdinal,Node> > domainMap = tpOp.getDomainMap();
600 const RCP<const Tpetra::Map<LocalOrdinal,GlobalOrdinal,Node> > rangeMap = tpOp.getRangeMap();
602 size_t maxRowSize = tpOp.getNodeMaxNumRowEntries();
603 if (maxRowSize == Teuchos::as<size_t>(-1))
606 std::vector<Scalar> scaledVals(maxRowSize);
607 if (tpOp.isFillComplete())
610 if (Op.isLocallyIndexed() ==
true) {
611 Teuchos::ArrayView<const LocalOrdinal> cols;
612 Teuchos::ArrayView<const Scalar> vals;
614 for (
size_t i = 0; i < rowMap->getNodeNumElements(); ++i) {
615 tpOp.getLocalRowView(i, cols, vals);
616 size_t nnz = tpOp.getNumEntriesInLocalRow(i);
617 if (nnz > maxRowSize) {
619 scaledVals.resize(maxRowSize);
621 for (
size_t j = 0; j < nnz; ++j)
622 scaledVals[j] = vals[j]*scalingVector[i];
625 Teuchos::ArrayView<const Scalar> valview(&scaledVals[0], nnz);
626 tpOp.replaceLocalValues(i, cols, valview);
631 Teuchos::ArrayView<const GlobalOrdinal> cols;
632 Teuchos::ArrayView<const Scalar> vals;
634 for (
size_t i = 0; i < rowMap->getNodeNumElements(); ++i) {
635 GlobalOrdinal gid = rowMap->getGlobalElement(i);
636 tpOp.getGlobalRowView(gid, cols, vals);
637 size_t nnz = tpOp.getNumEntriesInGlobalRow(gid);
638 if (nnz > maxRowSize) {
640 scaledVals.resize(maxRowSize);
643 for (
size_t j = 0; j < nnz; ++j)
644 scaledVals[j] = vals[j]*scalingVector[i];
647 Teuchos::ArrayView<const Scalar> valview(&scaledVals[0], nnz);
648 tpOp.replaceGlobalValues(gid, cols, valview);
653 if (doFillComplete) {
654 if (domainMap == Teuchos::null || rangeMap == Teuchos::null)
655 throw Exceptions::RuntimeError(
"In Utilities::Scaling: cannot fillComplete because the domain and/or range map hasn't been defined");
657 RCP<Teuchos::ParameterList> params = rcp(
new Teuchos::ParameterList());
658 params->set(
"Optimize Storage", doOptimizeStorage);
659 params->set(
"No Nonlocal Changes",
true);
660 Op.fillComplete(Op.getDomainMap(), Op.getRangeMap(), params);
663 throw Exceptions::RuntimeError(
"Only Tpetra::CrsMatrix types can be scaled (Err.1)");
666 throw Exceptions::RuntimeError(
"Matrix scaling is not possible because Tpetra has not been compiled with support for LO=GO=int.");
669 throw Exceptions::RuntimeError(
"Matrix scaling is not possible because Tpetra has not been enabled.");
673 static void MyOldScaleMatrix_Epetra (Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>& Op,
const Teuchos::ArrayRCP<Scalar>& scalingVector,
bool doFillComplete,
bool doOptimizeStorage) {
674 #ifdef HAVE_MUELU_EPETRA 677 const Epetra_CrsMatrix& epOp = Op2NonConstEpetraCrs(Op);
679 Epetra_Map
const &rowMap = epOp.RowMap();
684 for (
int i = 0; i < rowMap.NumMyElements(); ++i) {
685 epOp.ExtractMyRowView(i, nnz, vals, cols);
686 for (
int j = 0; j < nnz; ++j)
687 vals[j] *= scalingVector[i];
691 throw Exceptions::RuntimeError(
"Only Epetra_CrsMatrix types can be scaled");
694 throw Exceptions::RuntimeError(
"Matrix scaling is not possible because Epetra has not been enabled.");
695 #endif // HAVE_MUELU_EPETRA 703 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()) {
704 switch (Op.getRowMap()->lib()) {
705 case Xpetra::UseTpetra:
707 #ifdef HAVE_MUELU_TPETRA 708 #ifdef HAVE_MUELU_TPETRA_INST_INT_INT 713 RCP<Tpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > A;
714 Tpetra::RowMatrixTransposer<Scalar, LocalOrdinal, GlobalOrdinal, Node> transposer(rcpFromRef(tpetraOp),label);
715 A = transposer.createTranspose();
716 RCP<Xpetra::TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > AA = rcp(
new Xpetra::TpetraCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node>(A));
717 RCP<Xpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > AAA = rcp_implicit_cast<Xpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> >(AA);
718 RCP<Xpetra::CrsMatrixWrap<Scalar, LocalOrdinal, GlobalOrdinal, Node> > AAAA = rcp(
new Xpetra::CrsMatrixWrap<Scalar, LocalOrdinal, GlobalOrdinal, Node> (AAA));
722 catch (std::exception& e) {
723 std::cout <<
"threw exception '" << e.what() <<
"'" << std::endl;
724 throw Exceptions::RuntimeError(
"Utilities::Transpose failed, perhaps because matrix is not a Crs matrix");
727 throw Exceptions::RuntimeError(
"Utilities::Transpose: Tpetra is not compiled with LO=GO=int. Add TPETRA_INST_INT_INT:BOOL=ON to your configuration!");
730 throw Exceptions::RuntimeError(
"Utilities::Transpose: Tpetra is not compiled!");
734 case Xpetra::UseEpetra:
736 #if defined(HAVE_MUELU_EPETRA) && defined(HAVE_MUELU_EPETRAEXT) 737 Teuchos::TimeMonitor tm(*Teuchos::TimeMonitor::getNewTimer(
"ZZ Entire Transpose"));
741 Epetra_CrsMatrix * A =
dynamic_cast<Epetra_CrsMatrix*
>(&transposer(epetraOp));
744 RCP<Epetra_CrsMatrix> rcpA(A);
745 RCP<Xpetra::EpetraCrsMatrixT<GlobalOrdinal,Node> > AA = rcp(
new Xpetra::EpetraCrsMatrixT<GlobalOrdinal,Node> (rcpA));
746 RCP<Xpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> > AAA = rcp_implicit_cast<Xpetra::CrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> >(AA);
747 RCP<Xpetra::CrsMatrixWrap<Scalar, LocalOrdinal, GlobalOrdinal, Node> > AAAA = rcp(
new Xpetra::CrsMatrixWrap<Scalar, LocalOrdinal, GlobalOrdinal, Node>(AAA));
748 AAAA->fillComplete(Op.getRangeMap(), Op.getDomainMap());
752 throw Exceptions::RuntimeError(
"Epetra (Err. 2)");
757 throw Exceptions::RuntimeError(
"Only Epetra and Tpetra matrices can be transposed.");
761 return Teuchos::null;
770 template <
class View,
unsigned AppendValue >
776 template <
class DataType,
class Arg1,
class Arg2,
unsigned U,
unsigned T >
777 struct AppendTrait<
Kokkos::View< DataType, Arg1, Arg2, Kokkos::MemoryTraits<U> >, T> {
782 template <
class DataType,
class Arg1,
unsigned U,
unsigned T >
783 struct AppendTrait<
Kokkos::View< DataType, Arg1, Kokkos::MemoryTraits<U>, void >, T> {
789 template <
class DataType,
unsigned U,
unsigned T >
790 struct AppendTrait<
Kokkos::View< DataType, Kokkos::MemoryTraits<U>, void, void >, T> {
795 template <
class DataType,
class Arg1,
class Arg2,
unsigned T >
796 struct AppendTrait<
Kokkos::View< DataType, Arg1, Arg2, void >, T> {
801 template <
class DataType,
class Arg1,
unsigned T >
802 struct AppendTrait<
Kokkos::View< DataType, Arg1, void, void >, T> {
807 template <
class DataType,
unsigned T >
808 struct AppendTrait<
Kokkos::View< DataType, void, void, void >, T> {
815 #define MUELU_UTILITIES_KOKKOS_SHORT 817 #endif // #if defined(HAVE_MUELU_KOKKOS_REFACTOR) 819 #endif // MUELU_UTILITIES_KOKKOS_DECL_HPP 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 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 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< Epetra_CrsMatrix > Op2NonConstEpetraCrs(RCP< Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > 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 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 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.
Namespace for MueLu classes and methods.
static RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > MV2NonConstTpetraMV(RCP< Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > vec)
static void PauseForDebugger()
Kokkos::View< const bool *, typename NO::device_type > DetectDirichletRows(const Xpetra::Matrix< SC, LO, GO, NO > &A, const typename Teuchos::ScalarTraits< SC >::magnitudeType &tol)
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< Epetra_MultiVector > MV2NonConstEpetraMV(RCP< Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > vec)
static RCP< Tpetra::CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op2NonConstTpetraCrs(RCP< Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op)
static RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > MV2NonConstTpetraMV2(Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &vec)
static RCP< const Tpetra::RowMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op2TpetraRow(RCP< const Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Op)
static Teuchos::ArrayRCP< Scalar > GetLumpedMatrixDiagonal(const Matrix &A)
Extract Matrix Diagonal of lumped matrix.
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 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 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)
static RCP< Teuchos::FancyOStream > MakeFancy(std::ostream &os)
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 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< 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.