I have a fortran 90 code and I want to solve a dense linear system with GMRES. I would prefer the restarted GMRES with preconditioning. Is there some library that you know of that I could use? Now I am solving the system with LU using
dgetrs from LAPACK.
I want to add also that the matrix is square, with diagonal dominance and invertible.