Hi all,<br><br>I have a matrix M, for which I set up a Cholesky decomposition M = L L^T  using MUMPS and PETSc. I then obtained the factored matrix with<br><br>  PCFactorGetMatrix(pc,&factor);  <br>
<br>However, I do not want to solve M*x = b, which seems to be what MatSolve does. How do I tell it to solve Lx = b instead? <br><br>Thanks,<br><br>Pearl Flath<br><br>