Personal tools
You are here: Home / openCMISS / Wiki / Solver Libraries
Log in

Forgot your password?

Solver Libraries


The "Trilinos Project": "is an effort to develop and implement robust parallel algorithms using modern object-oriented software design, while still leveraging the value of established numerical libraries such as PETSc, Aztec, the BLAS and LAPACK."

Includes multilevel preconditioners and nonlinear solver. Note that Jon Pearce uses this in his thesis work.

Epetra in Trilinos has experimental support for "shared memory parallelism within MPI processes":


"hypre": is "a library of high

performance preconditioners that features parallel multigrid methods for both structured and unstructured grid problems".

A subset ("most": of hypre is threaded.

Note that hypre contains BoomerAMG, which is an MPI-based Algebraic Multigrid Code. Also, Oliver and I (Travis) have some experience using hypre and have good connections to developers of hypre.

The "licence": states: "Commercialization of this product is prohibited without notifying the Department of Energy (DOE) or Lawrence Livermore National Laboratory (LLNL)."


"Aztec": "is a library that provides algorithms for the solution of large sparse linear systems arising in scientific and engineering applications".

Written in C. Callable from Fortran. Uses MPI. "No threading currently":

The "licence": is for "research, development and demonstration purposes only".


"SUNDIALS": is a SUite of Nonlinear and

DIfferential/ALgebraic equation Solvers.

Developed from Fortran codes. Now C. Note that Jon Pearce is using this for this thesis work.


"UbikSolve": is a F90/95 sparse MPI-based linear algebra library, with emphasis on preconditioned Krylov subspace methods for the solution of linear systems of equations.

I (Travis) used this at Los Alamos and found it easy to use. Also, it has an Artistic License, like perl. Furthermore, you pass a pointer to your Matrix-Vector multiply into the Krylov solver so it is easy to use with respect to different distributed matrix storage schemes.


Petsc can also be thought of as a linear solver library. In fact, it contains a variety of solver options that can probably be taken advantage of independently of the rest of Petsc.