81 typedef Xpetra::BlockedCrsMatrix<Scalar, LocalOrdinal, GlobalOrdinal, Node> BlockCrs;
83 RCP<Matrix> originalA = Get< RCP<Matrix> >(currentLevel,
"A");
84 RCP<BlockCrs> A = Teuchos::rcp_dynamic_cast<BlockCrs>(originalA);
87 TEUCHOS_TEST_FOR_EXCEPTION(A==Teuchos::null,
Exceptions::BadCast,
"MueLu::RepartitionBlockDiagonalFactory::Build: input matrix A is not of type BlockedCrsMatrix! error.");
90 RCP<BlockCrs> DiagonalMatrix = Teuchos::rcp(
new BlockCrs(A->getBlockedRangeMap(),A->getBlockedDomainMap(),0));
91 for(
size_t i=0; i< A->Rows(); i++)
92 DiagonalMatrix->setMatrix(i,i,A->getMatrix(i,i));
94 Set(currentLevel,
"A",Teuchos::rcp_dynamic_cast<Matrix>(DiagonalMatrix));