krylov resid

Author: vpuri3

src/iterative_wrappers.jl

@@ -137,7 +137,7 @@ function SciMLBase.solve(cache::LinearCache, alg::KrylovJL; kwargs...)
 
     args = (cache.cacheval, cache.A, cache.b)
     kwargs = (atol = atol, rtol = rtol, itmax = itmax, verbose = verbose,
-              alg.kwargs...)
+              history = true, alg.kwargs...)
 
     if cache.cacheval isa Krylov.CgSolver
         N !== I &&
@@ -159,7 +159,7 @@ function SciMLBase.solve(cache::LinearCache, alg::KrylovJL; kwargs...)
         Krylov.solve!(args...; kwargs...)
     end
 
-    resid = nothing
+    resid = cache.cacheval.stats.residuals |> last
 
     return SciMLBase.build_linear_solution(alg, cache.u, resid, cache;
                                            iters = cache.cacheval.stats.niter)