pass resid from solve!

Author: utkarsh530

src/solve.jl

@@ -2,14 +2,8 @@ function SciMLBase.solve(prob::NonlinearProblem,
                          alg::AbstractNonlinearSolveAlgorithm, args...;
                          kwargs...)
   solver = init(prob, alg, args...; kwargs...)
-  sol = solve!(solver)
+  sol, resid = solve!(solver)
   if typeof(sol) <: NewtonSolution
-    resid = zero(prob.u0)
-    if isinplace(prob)
-      prob.f(resid,sol.u,prob.p)
-    else
-      resid = prob.f(sol.u,prob.p)
-    end
     return SciMLBase.build_solution(prob, alg, sol.u, resid;retcode=:Success)
   end
 end
@@ -76,7 +70,7 @@ function SciMLBase.solve!(solver::AbstractImmutableNonlinearSolver)
     @set! solver.retcode = MAXITERS_EXCEED
   end
   sol = get_solution(solver)
-  return sol
+  return sol, solver.fu
 end
 
 """