[petsc-dev] Performance of Petsc Cholesky vs LU

[Jed Brown]

On Wed, Nov 30, 2011 at 13:57, Barry Smith <bsmith at mcs.anl.gov> wrote:

> Right, so why not just use the AIJ format if you plan to use a direct
> Cholesky solver; only the upper triangular part of the factor is actually
> stored so the memory usage of AIJ for Cholesky is just the same as as of
> SBAIJ for Cholesky.,
>
>
>   In other words, if you want to use full Cholesky just use AIJ.
>

Sure, I don't know how many people only have preallocation for SBAIJ.

What I care about is that -pc_type cholesky should do something halfway
reasonable, certainly when called with an AIJ matrix for which ND makes
sense.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20111130/aabd47e0/attachment.html>