2 #ifndef PACKAGES_MUELU_ADAPTERS_AZTECOO_MUELU_AZTECEPETRAOPERATOR_CPP_ 3 #define PACKAGES_MUELU_ADAPTERS_AZTECOO_MUELU_AZTECEPETRAOPERATOR_CPP_ 5 #include "Xpetra_EpetraMultiVector.hpp" 7 #include "MueLu_config.hpp" 8 #include "MueLu_RefMaxwell.hpp" 12 #if defined(HAVE_MUELU_SERIAL) and defined(HAVE_MUELU_EPETRA) 19 const Xpetra::EpetraMultiVectorT<GO,NO> eX(Teuchos::rcpFromRef(const_cast<Epetra_MultiVector&>(X)));
20 Xpetra::EpetraMultiVectorT<GO,NO> eY(Teuchos::rcpFromRef(Y));
24 Teuchos::RCP<Xpetra::EpetraMultiVectorT<GO,NO> > tmpY = Teuchos::rcp(
new Xpetra::EpetraMultiVectorT<GO,NO>(eY.getMap(), eY.getNumVectors()));
26 xOp_->apply(eX, *tmpY);
28 eY.update(1.0, *tmpY, 0.0);
35 }
catch (std::exception& e) {
37 std::cerr <<
"Caught an exception in MueLu::AztecEpetraOperator::ApplyInverse():" << std::endl
38 << e.what() << std::endl;
44 const Epetra_Comm& AztecEpetraOperator::Comm()
const {
45 const Epetra_Map emap = Xpetra::toEpetra(xOp_->getDomainMap());
49 const Epetra_Map& AztecEpetraOperator::OperatorDomainMap()
const {
52 RCP<const Xpetra::CrsMatrixWrap<SC,LO,GO,NO> > crsOp = rcp_dynamic_cast<
const Xpetra::CrsMatrixWrap<SC,LO,GO,NO> >(A);
53 if (crsOp == Teuchos::null)
54 throw Exceptions::BadCast(
"Cast from Xpetra::Matrix to Xpetra::CrsMatrixWrap failed");
55 const RCP<const Xpetra::EpetraCrsMatrixT<GO,NO>> &tmp_ECrsMtx = rcp_dynamic_cast<
const Xpetra::EpetraCrsMatrixT<GO,NO>>(crsOp->getCrsMatrix());
56 if (tmp_ECrsMtx == Teuchos::null)
57 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::EpetraCrsMatrix failed");
58 return tmp_ECrsMtx->getEpetra_CrsMatrixNonConst()->DomainMap();
63 Teuchos::RCP<const Map> map = xOp_->getDomainMap();
64 return Xpetra::toEpetra(map);
67 const Epetra_Map & AztecEpetraOperator::OperatorRangeMap()
const {
71 RCP<const Xpetra::CrsMatrixWrap<SC,LO,GO,NO> > crsOp = rcp_dynamic_cast<
const Xpetra::CrsMatrixWrap<SC,LO,GO,NO> >(A);
72 if (crsOp == Teuchos::null)
73 throw Exceptions::BadCast(
"Cast from Xpetra::Matrix to Xpetra::CrsMatrixWrap failed");
74 const RCP<const Xpetra::EpetraCrsMatrixT<GO,NO>> &tmp_ECrsMtx = rcp_dynamic_cast<
const Xpetra::EpetraCrsMatrixT<GO,NO>>(crsOp->getCrsMatrix());
75 if (tmp_ECrsMtx == Teuchos::null)
76 throw Exceptions::BadCast(
"Cast from Xpetra::CrsMatrix to Xpetra::EpetraCrsMatrix failed");
77 return tmp_ECrsMtx->getEpetra_CrsMatrixNonConst()->RangeMap();
82 Teuchos::RCP<const Map> map = xOp_->getRangeMap();
83 return Xpetra::toEpetra(map);
Namespace for MueLu classes and methods.
Preconditioner (wrapped as a Xpetra::Operator) for Maxwell's equations in curl-curl form...
const Epetra_Comm & Comm() const