[Trilinos-Users] MueLu with Tpetra::BlockCrsMatrix

[Drumm, Clif]

I have been unable to create a Tpetra Preconditioner for MueLu using a Tpetra::BlockCrsMatrix. I've followed the example

BlockCrs.cpp in the muelu examples folder


RCP<Tpetra_BlockCrsMatrix> Ablock = rcp(new Tpetra_BlockCrsMatrix(crsGraph, 2)

  fill Ablock with sumIntoLocalValues(...)

RCP<Xpetra_CrsMatrix> Axt = rcp(new Xpetra_TpetraBlockCrsMatrix(Ablock));

RCP<Xpetra_Matrix> A = rcp(new Xpetra_CrsMatrixWrap(Axt));

RCP<Tpetra_Operator>    At = MueLu::Utilities<Scalar,LocalOrdinal,GlobalOrdinal,Node>::Op2NonConstTpetraRow(A);
RCP<Tpetra_Operator>    Mt = MueLu::CreateTpetraPreconditioner(At,MueLuList);


The error that I am getting is:


Setup Smoother (MueLu::Amesos2Smoother{type = Klu})
terminate called after throwing an instance of 'MueLu::Exceptions::BadCast'
  what():  Cast from Xpetra::CrsMatrix to Xpetra::TpetraCrsMatrix failed


I'm using version 12.18.1


Thank you for any help you can provide.

Clif Drumm

Radiation Effects Theory Dept 01341

Sandia National Labs
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://trilinos.org/pipermail/trilinos-users_trilinos.org/attachments/20191204/af86d7b4/attachment.html>