[Trilinos-Users] CrsMatrix Transpose

[Pate Motter]

Hi,

I am currently having segmentation faults when transposing a CrsMatrix
using the RowMatrixTransposer method mentioned in a previous email. For my
demo it it works for the serial version as well as mpiruns with 2 and 3
nodes. However, when running at > 3 the program segfaults at
createTranspose().

It's completely possible that I may not have the problem set up correctly.

int main(int argc, char *argv[]) {
Teuchos::GlobalMPISession mpiSession(&argc, &argv);
Platform& platform = Tpetra::DefaultPlatform::getDefaultPlatform();
RCP<const Teuchos::Comm<int> > comm = rcp (new Teuchos::MpiComm<int>
(MPI_COMM_WORLD));
RCP<NT> node = platform.getNode();

std::string filename("../large.mtx");
RCP<MAT> A, B;
A = Tpetra::MatrixMarket::Reader<MAT>::readSparseFile(filename, comm, node,
true);
Tpetra::RowMatrixTransposer<ST, LO, GO, NT> transposer(A);
B = transposer.createTranspose();
}


Thank you,
Pate Motter
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://software.sandia.gov/pipermail/trilinos-users/attachments/20150222/a9bce4f6/attachment.html>