<div dir="ltr">Dear all,<div><br></div><div>I am trying to solve a linear system for a symmetric matrix with MUMPS.</div><div><br></div><div>Is there a way to tell MUMPS that the matrix is indeed symmetric?</div><div><br></div><div>The way I build the matrix is</div><div><br></div><div><span style="font-family:monospace,monospace">Mat A,AT,ATA</span><br></div><div><font face="monospace, monospace">MatHermitianTranspose(A,MAT_INITIAL_MATRIX,&AT);<br></font></div><div><font face="monospace, monospace">MatMatMult(AT,A,MAT_INITIAL_MATRIX,7,&ATA);<br></font></div><div><font face="monospace, monospace">MatSetOption(ATA,MAT_SYMMETRY_ETERNAL,PETSC_TRUE);<br></font></div><div><br></div><div>but MUMPS returns </div><div><br></div><div><font face="monospace, monospace">L U Solver for unsymmetric matrices</font><br></div><div><font face="monospace, monospace"><br></font></div><div><font face="arial, helvetica, sans-serif">Of course, any suggestion of a better/more efficient way to build ATA or store only half of it, that is more than welcome.</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">Thanks for your help,</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">Gianluca</font></div></div>