[petsc-users] MUMPS with symmetric matrices

[Gianluca Meneghello]

Dear all,

I am trying to solve a linear system for a symmetric matrix with MUMPS.

Is there a way to tell MUMPS that the matrix is indeed symmetric?

The way I build the matrix is

Mat A,AT,ATA
MatHermitianTranspose(A,MAT_INITIAL_MATRIX,&AT);
MatMatMult(AT,A,MAT_INITIAL_MATRIX,7,&ATA);
MatSetOption(ATA,MAT_SYMMETRY_ETERNAL,PETSC_TRUE);

but MUMPS returns

L U Solver for unsymmetric matrices

Of course, any suggestion of a better/more efficient way to build ATA or
store only half of it, that is more than welcome.

Thanks for your help,

Gianluca
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20151104/479aee70/attachment.html>