Tibo
Tibo

Reputation: 305

Iterative Solver Library in Fortran

I wrote a code which solves large systems of PDEs using some discretization methods which basically involves solving large, sparse systems Ax=b many times at each time steps.

I currently use the PARDISO solver (from the intel MKL library) which is a direct LU factorization of A to solve the system. I would like to compare this method with the use of iterative solvers (which, with the use of preconditioners, might perform better since I could use the same preconditioner over many time steps if my Jacobian matrix does not vary too much).

My question is then, what library do you suggest for sparse iterative solvers in fortran ? I found one (SLATEC) which was written in 1993 so I wonder if there is something more performant which was written more recently ?

Thanks :)

Upvotes: 3

Views: 2833

Answers (2)

oz123
oz123

Reputation: 28868

I would also add:
LIS
AGMG

Oh, well ... the complete list of Linear Algebra Solvers

Upvotes: 1

Tibo
Tibo

Reputation: 305

Thanks for the comments, PETSc seems to be exactly what I was looking for, just need to learn how to link C calls into fortran now :)

Upvotes: 0

Related Questions