From 65bcea052a9fd20633fdbed6c1c43cb910fa8b4b Mon Sep 17 00:00:00 2001
From: Christopher Rackauckas <accounts@chrisrackauckas.com>
Date: Tue, 17 Oct 2023 04:12:33 -0400
Subject: [PATCH] Set ITP as the default method for interval problems

---
 .../src/SimpleNonlinearSolve.jl                     | 13 +++++++++++++
 1 file changed, 13 insertions(+)

diff --git a/lib/SimpleNonlinearSolve/src/SimpleNonlinearSolve.jl b/lib/SimpleNonlinearSolve/src/SimpleNonlinearSolve.jl
index 5c500d500..85052e5ff 100644
--- a/lib/SimpleNonlinearSolve/src/SimpleNonlinearSolve.jl
+++ b/lib/SimpleNonlinearSolve/src/SimpleNonlinearSolve.jl
@@ -46,6 +46,19 @@ include("batched/raphson.jl")
 include("batched/dfsane.jl")
 include("batched/broyden.jl")
 
+## Default algorithm
+
+# Set the default bracketing method to ITP
+
+function SciMLBase.solve(prob::IntervalNonlinearProblem; kwargs...)
+    SciMLBase.solve(prob, ITP(); kwargs...)
+end
+
+function SciMLBase.solve(prob::IntervalNonlinearProblem, alg::Nothing,
+    args...; kwargs...)
+    SciMLBase.solve(prob, ITP(), args...; kwargs...)
+end
+
 import PrecompileTools
 
 PrecompileTools.@compile_workload begin