[petsc-users] Speed up KSPSolve of dense matrix

[Jed Brown]

Barry Smith <bsmith at mcs.anl.gov> writes:
>> I've tried to use a direct solver like suggested on pp 72, but:
>> 
>> ./petBench 600 1 -ksp_type preonly -pc_type lu
>
>   You cannot use LU with SBAIJ format. Only Cholesky. So use -pc_type cholesky

And using a sparse matrix format for a dense matrix is a recipe for poor performance.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20141104/0b4771ba/attachment.pgp>