[petsc-users] Cholesky matrix factorization

[Pearl Flath]

Hi all,

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

  PCFactorGetMatrix(pc,&factor);

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?

Thanks,

Pearl Flath
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120908/f859efbe/attachment.html>