Refine ITP return solution on hitting f(·) zero

DaniGlez · DaniGlez · commit cd79c8fd6374 · 2023-07-23T14:45:08.000+02:00
diff --git a/src/itp.jl b/src/itp.jl
@@ -125,8 +125,8 @@ function SciMLBase.solve(prob::IntervalNonlinearProblem, alg::ITP,
             fl = yp
         else
             return SciMLBase.build_solution(prob, alg, xp, yps;
-                                            retcode = ReturnCode.Success, left = left,
-                                            right = right)
+                                            retcode = ReturnCode.Success, left = xp,
+                                            right = xp)
         end
         i += 1
         mid = (left + right) / 2
