
<html><body><div style="font-family: times new roman, new york, times, serif; font-size: 12pt; color: #000000"><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(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/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 ?<br></div><div><br></div><div>Franck<br></div></div></body></html>