<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Apr 8, 2015 at 11:06 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">David:<span><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"><span><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"><span><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><br></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></span><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></div></div></blockquote><div><br></div></span><div>When user requests Cholesky factorization, PETSc sets  "SYM=2" as default; If user sets flag of spd, then the interface uses "SYM=1" in order to get an LDL^T factorization. See MatGetFactor_xxx_mumps() in</div><div>petsc/src/mat/impls/aij/mpi/mumps/mumps.c:</div><div><br></div><div>...</div><div> B->factortype = MAT_FACTOR_CHOLESKY;<br></div><div><div>  if (A->spd_set && A->spd) mumps->sym = 1;</div><div>  else                      mumps->sym = 2;</div></div><div>...</div></div></div></div></blockquote><div><br></div><div><br></div><div>OK, got it, thanks!</div><div><br></div><div>David</div><div><br></div><div></div></div><br></div><div class="gmail_extra"><br></div></div>