From e1301cb4de60d2a1135be285ecf126d80492ef4c 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/src/SimpleNonlinearSolve.jl b/src/SimpleNonlinearSolve.jl
index 5c500d5..85052e5 100644
--- a/src/SimpleNonlinearSolve.jl
+++ b/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