Barry Smith <bsm...@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.
signature.asc
Description: PGP signature