[Trilinos-Users] trouble with Epetra_SerialDenseMatrix

[Nikhil Kriplani]

Hello,

I am encountering some strange behavior while trying to use
Epetra_SerialDenseMatrix. I am creating a new matrix in my code, by
doing something like
Tc = Epetra_SerialDenseMatrix(rows, cols)
However, my program crashes at this point. When I run the debugger, it
shows that control is indeed transferred to the Serial matrix
constructor (Epetra_SerialDenseMatrix.cpp line 54), the Shape function
gets called, and returns back to the matrix constructor.

At this point, for some reason the Matrix destructor function is
called (which calls CleanUpData()) and I get a segmentation fault in
libc6. It is also strange, that in my same source file, I have two
other SerialDenseMatrices just like this, and they work fine.

Sorry if this does not make any sense, but I am really quite confused.
I am using Trilinos version 9.0.2 on Linux AMD64, gcc-4.3 and get the
same behavior with version 8.0.5.

Another question: Would you recommend using Teuchos Dense
Vectors/Matrices instead of the Epetra equivalents?

Thanks,
Nik