# Library to solve dense linear system with GMRES

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 dgetrf and dgetrs from LAPACK.

I want to add also that the matrix is square, with diagonal dominance and invertible.

Thanks!