[petsc-users] Cholesky matrix factorization

[Jed Brown]

On Sat, Sep 8, 2012 at 3:19 PM, Pearl Flath <pflath at ices.utexas.edu> wrote:

> 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?
>

Looks like you want MatForwardSolve().

Note that there is also PCApplySymmetric{Left,Right}.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120908/f28d9075/attachment.html>