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