LinearSolve.jl
LinearSolve.jl copied to clipboard
`sol.resid` not returning the correct thing
julia> using LinearSolve, LinearAlgebra; A=rand(4,4); b=rand(4); sol=solve(LinearProblem(A,b),KrylovJL_GMRES()); norm(A*sol.u - b)
2.546262077003306e-15
julia> sol.resid
4
Oh, did I put the iter there? https://github.com/JuliaSmoothOptimizers/Krylov.jl/issues/476