Update the wrapper of Krylov.jl

Author: amontoison

Project.toml

@@ -29,7 +29,7 @@ FastLapackInterface = "1"
 GPUArraysCore = "0.1"
 IterativeSolvers = "0.9.2"
 KLU = "0.3.0, 0.4"
-Krylov = "0.7.11, 0.8"
+Krylov = "0.9"
 KrylovKit = "0.5, 0.6"
 RecursiveFactorization = "0.2.8"
 Reexport = "1"

src/iterative_wrappers.jl

@@ -17,6 +17,7 @@ end
 
 KrylovJL_CG(args...; kwargs...) = KrylovJL(args...; KrylovAlg = Krylov.cg!, kwargs...)
 KrylovJL_GMRES(args...; kwargs...) = KrylovJL(args...; KrylovAlg = Krylov.gmres!, kwargs...)
+KrylovJL_FGMRES(args...; kwargs...) = KrylovJL(args...; KrylovAlg = Krylov.fgmres!, kwargs...)
 function KrylovJL_BICGSTAB(args...; kwargs...)
     KrylovJL(args...; KrylovAlg = Krylov.bicgstab!, kwargs...)
 end
@@ -89,6 +90,10 @@ function get_KrylovJL_solver(KrylovAlg)
         Krylov.QmrSolver
     elseif (KrylovAlg === Krylov.gmres!)
         Krylov.GmresSolver
+    elseif (KrylovAlg === Krylov.fgmres!)
+        Krylov.FgmresSolver
+    elseif (KrylovAlg === Krylov.gpmr!)
+        Krylov.GpmrSolver
     elseif (KrylovAlg === Krylov.fom!)
         Krylov.FomSolver
     end
@@ -105,6 +110,8 @@ function init_cacheval(alg::KrylovJL, A, b, u, Pl, Pr, maxiters::Int, abstol, re
     solver = if (alg.KrylovAlg === Krylov.dqgmres! ||
                  alg.KrylovAlg === Krylov.diom! ||
                  alg.KrylovAlg === Krylov.gmres! ||
+                 alg.KrylovAlg === Krylov.fgmres! ||
+                 alg.KrylovAlg === Krylov.gpmr! ||
                  alg.KrylovAlg === Krylov.fom!)
         KS(A, b, memory)
     elseif (alg.KrylovAlg === Krylov.minres! ||