Match abstol with the rest of NonlinearSolve

src/simple.jl

@@ -50,7 +50,7 @@ function ___kinsol(f,
     linear_solver = :Dense,
     jac_upper = 0,
     jac_lower = 0,
-    abstol::Float64 = 1e-6)
+    abstol::Float64 = eps(Float64) ^ (4 // 5))
     # f, Function to be optimized of the form f(y::Vector{Float64}, fy::Vector{Float64})
     #    where `y` is the input vector, and `fy` is the result of the function
     # y0, Vector of initial values