<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Apr 8, 2015 at 9:28 PM, Hong <span dir="ltr"><<a href="mailto:hzhang@mcs.anl.gov" target="_blank">hzhang@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">David ,<span class=""><br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div><div><div><span style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:12px">I'm curious about how to do a symmetric LDL^T factorization (instead of LU) with MUMPS and SuperLU. Based on this example:</span></div><div><span style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:12px"><br></span></div><div><font color="#000000" face="Calibri, sans-serif"><span style="font-size:12px"><a href="http://www.mcs.anl.gov/petsc/petsc-3.4/src/ksp/ksp/examples/tutorials/ex52.c.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-3.4/src/ksp/ksp/examples/tutorials/ex52.c.html</a></span></font><br></div><div><font color="#000000" face="Calibri, sans-serif"><span style="font-size:12px"><br></span></font></div><div><font color="#000000" face="Calibri, sans-serif"><span style="font-size:12px">my understanding is as follows:</span></font></div><div><font color="#000000" face="Calibri, sans-serif"><span style="font-size:12px"><br></span></font></div><div><font color="#000000" face="Calibri, sans-serif"><span style="font-size:12px">- With MUMPS I gather that we need to specify:</span></font></div><div><font color="#000000" face="Calibri, sans-serif"><span style="font-size:12px"><div><br></div><div>MatSetOption(A,MAT_SPD,PETSC_TRUE);</div><div>PCSetType(pc,PCCHOLESKY);</div><div><br></div><div>I guess "-pc_type cholesky" on the command line is equivalent to the PCSetType call, right? Is specifying MAT_SPD required in order for MUMPS to do an LDL^T factorization?</div></span></font></div></div></div></div></blockquote></span><div>Mumps supports Cholesky factorization for symmetric, and symmetric+spd matrices.  You may consult mumps user manual.</div></div></div></div></blockquote><div><br></div><div><br></div><div>OK, thanks. But I was also wondering about the PETSc interface to MUMPS.</div><div><br></div><div>I consulted the MUMPS manual, but it just says that you need to specify "SYM=1" in order to get an LDL^T factorization. I gather that if we set MatSetOption(A,MAT_SPD,PETSC_TRUE); in PETSc then that will tell MUMPS to use SYM=1, right?</div><div><br></div><div>Thanks,</div><div>David</div><div><br></div><div><br></div></div></div></div>