Exemplo n.º 1
0
const Epetra_Map& EpetraOperator::OperatorDomainMap() const {
    RCP<Matrix> A = Hierarchy_->GetLevel(0)->Get<RCP<Matrix> >("A");

    RCP<Xpetra::BlockedCrsMatrix<double, int, int> > epbA = Teuchos::rcp_dynamic_cast<Xpetra::BlockedCrsMatrix<double, int, int> >(A);
    if (epbA != Teuchos::null)
        return Xpetra::toEpetra(epbA->getDomainMap());

    RCP<Epetra_CrsMatrix> epA = Utils::Op2NonConstEpetraCrs(A);
    return epA->DomainMap();
}