about KSP based on parallel dense matrix

Barry Smith bsmith at mcs.anl.gov
Mon Sep 22 21:02:11 CDT 2008

    The KSP methods are INDEPENDENT of the matrix format, they can all  
be used just as well with dense matrices
as sparse matrices. As with all linear operators the question is what  
are suitable preconditioners (if any) for a particular
problem. If the matrix is well conditioned then -pc_type jacobi would  
be fine; preconditioners like ILU are probably silly
for dense matrices since they are expensive.

On Sep 22, 2008, at 8:35 PM, Yujie wrote:

> To my knowledge, PETsc doesn't provide parallel dense matrix-based  
> solvers, such as for CG, GMRES and so on. If it is, how to deal with  
> this problem? Thanks.
> Regards,
> Yujie

More information about the petsc-users mailing list