<div dir="ltr">mumps stores <span style="color:rgb(0,0,0);font-family:"times new roman","new york",times,serif;font-size:16px"> LDLt using its private internal data structure, which is not exposed to user and can be used only via mumps API. </span><div><font color="#000000" face="times new roman, new york, times, serif"><span style="font-size:16px">Either ask mumps developer about it or use a debugger to step into mumps source code.</span></font></div><div><font color="#000000" face="times new roman, new york, times, serif"><span style="font-size:16px"><br></span></font></div><div><font color="#000000" face="times new roman, new york, times, serif"><span style="font-size:16px">Hong<br></span></font><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Aug 24, 2017 at 9:46 AM, Franck Houssen <span dir="ltr"><<a href="mailto:franck.houssen@inria.fr" target="_blank">franck.houssen@inria.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div style="font-family:"times new roman","new york",times,serif;font-size:12pt;color:rgb(0,0,0)"><div>How to get L (as a Mat) and D (probably as a Vec) from cholesky factorisation ?</div><div><br></div><div>I googled and looked in the doc but didn't find any obvious answer. I am not sure to get the steps to follow: is there a dedicated example ?<br></div><div><br></div><div>To make it short, I have:<br></div><div>PCSetType(pc, PCCHOLESKY); // LDLt<br></div><div>PCFactorSetMatSolverPackage(<wbr>pc, "mumps");</div><div>PCFactorGetMatrix(pc, &f);</div><div>MatGetInertia(f, ....); // OK !<br></div><div><br></div><div>Now, I would like to get/see L and D: how to get them ?<br></div><div><br></div><div>From petsc/src/ksp/ksp/examples/<wbr>tutorials/ex52.c, I would say MatGetDiagonal(f, D) returns D: correct ? If so, how to get L as well ?<br></div><div><br></div><div>A related question would be : how to get L and U from a PCLU factorisation ?<span class="gmail-HOEnZb"><font color="#888888"><br></font></span></div><span class="gmail-HOEnZb"><font color="#888888"><div><br></div><div>Franck<br></div></font></span></div></div></blockquote></div><br></div></div></div>