Summary of Sparse Linear Solvers Available from PETSc
Requests and contributions welcome
| Algorithm | Matrix types | External Package* | Parallel | Complex | ||
|---|---|---|---|---|---|---|
|
* External Package If this column is filled in you must install the external package indicated to use this feature. See instructions ** A "block Jacobi" style SOR with lits local smoothing followed by update of ghost points from other processors, see MatSOR() *** Use PCType of PCSOR |
||||||
| Preconditioners | Jacobi | aij, baij, sbaij, dense | --- | X | X | |
| point block Jacobi | aij, baij, baij, kaij, mpisell, is | --- | X | X | ||
| SOR | seqdense, seqaij, seqsbaij, mpiaij ** | --- | X | |||
| point block SOR*** | seqbaij, bs = 2,3,4,5 | --- | X | |||
| block Jacobi | aij, baij, sbaij | --- | X | X | ||
| additive Schwarz | aij, baij, sbaij | --- | X | X | ||
| deflation | --- | X | X | |||
| Incomplete factorizations | ILU(k) | seqaij, seqbaij | --- | X | ||
| ICC(k) | seqaij, seqbaij, seqsbaij | --- | X | |||
| ILU dt | seqaij | SuperLU Sequential ILU solver (LBNL). | X | |||
| aij | pilut/hypre (LLNL) Deprecated. Use euclid/hypre instead | X | ||||
| Matrix-free | infrastructure | X | X | |||
| Multigrid | infrastructure | X | X | |||
| geometric (DMDA for structured grid) | X | X | ||||
| smoothed aggregation | aij | --- | X | X | ||
| structured geometric | hyprestruct | PFMG from hypre | X | |||
| classical algebraic | aij | BoomerAMG/hypre (LLNL) | X | |||
| aij | ML/Trilinos (SNL) | X | ||||
| domain decomposition-based | aij, baij, sbaij, is | HPDDM (CNRS) | X | X | ||
| Physics-based splitting | relaxation and Schur-complement | aij, baij, nest | --- | X | X | |
| least squares commutator | schurcomplement | --- | X | X | ||
| Approximate inverses | aij | Parasails/hypre (LLNL) | X | |||
| aij | SPAI | X | ||||
| Substructuring | balancing Neumann-Neumann | is | --- | X | X | |
| BDDC | is | --- | X | X | ||
| Direct solvers | LU | LU | seqaij, seqbaij | --- | X | |
| seqaij | MATLAB | X | ||||
| aij | PaStiX (Inria) | X | X | |||
| aij | SuperLU (LBNL) SuperLU Sequential LU solver / SuperLU_DIST Parallel LU solver | X | X | |||
| aij, baij | MUMPS (CERFACS, CNRS, ENS Lyon, INP Toulouse, Inria, Mumps Technologies, University of Bordeaux) | X | X | |||
| seqaij | ESSL (IBM) | |||||
| seqaij | UMFPACK Part of SuiteSparse | X | ||||
| seqaij | KLU Part of SuiteSparse | X | ||||
| seqaij | LUSOL | |||||
| seqaij, seqbaij | MKL Pardiso (Intel) | X | ||||
| mpiaij, mpibaij | MKL CPardiso (Intel) | X | X | |||
| dense | Elemental | X | X | |||
| Cholesky | Cholesky | seqaij, seqsbaij | --- | X | ||
| sbaij | PaStiX (Inria) | X | X | |||
| sbaij | MUMPS (CERFACS, CNRS, ENS Lyon, INP Toulouse, Inria, Mumps Technologies, University of Bordeaux) | X | X | |||
| seqaij, seqsbaij | CHOLMOD Part of SuiteSparse | X | ||||
| dense | Elemental | X | X | |||
| seqsbaij | MKL Pardiso (Intel) | |||||
| mpisbaij | MKL CPardiso (Intel) | X | ||||
| QR | matlab | MATLAB | ||||
| XXt and XYt | aij | --- | X | |||
| Krylov methods | Richardson | --- | X | X | ||
| Chebyshev | --- | X | X | |||
| conjugate gradients | --- | X | X | |||
| GMRES | --- | X | X | |||
| bi-CG-stab | --- | X | X | |||
| transpose free QMR | --- | X | X | |||
| conjugate residuals | --- | X | X | |||
| conjugate gradient squared | --- | X | X | |||
| bi-conjugate gradient | --- | X | X | |||
| MINRES | --- | X | X | |||
| flexible GMRES | --- | X | X | |||
| LSQR (supports rectangular matrices) | --- | X | X | |||
| CG for least squares (supports rectangular matrices) | --- | X | X | |||
| SYMMLQ | --- | X | X | |||
| LGMRES | --- | X | X | |||
| GCR | --- | X | X | |||
| conjugate gradient on the normal equations | --- | X | X | |||
| GCR with inner orthogonalization and deflated restarting | HPDDM (CNRS) | X | X | |||
This page was inspired by Jack Dongarra's freely available software for linear algebra list