feat: add bracketing default algorithm

Author: avik-pal

lib/BracketingNonlinearSolve/src/BracketingNonlinearSolve.jl

@@ -19,6 +19,14 @@ include("falsi.jl")
 include("itp.jl")
 include("ridder.jl")
 
+# Default Algorithm
+function CommonSolve.solve(prob::IntervalNonlinearProblem; kwargs...)
+    return CommonSolve.solve(prob, ITP(); kwargs...)
+end
+function CommonSolve.solve(prob::IntervalNonlinearProblem, nothing, args...; kwargs...)
+    return CommonSolve.solve(prob, ITP(), args...; kwargs...)
+end
+
 @setup_workload begin
     for T in (Float32, Float64)
         prob_brack = IntervalNonlinearProblem{false}(

lib/BracketingNonlinearSolve/test/rootfind_tests.jl

@@ -7,7 +7,7 @@ end
 @testitem "Interval Nonlinear Problems" setup=[RootfindingTestSnippet] tags=[:core] begin
     using ForwardDiff
 
-    @testset for alg in (Bisection(), Falsi(), Ridder(), Brent(), ITP(), Alefeld())
+    @testset for alg in (Bisection(), Falsi(), Ridder(), Brent(), ITP(), Alefeld(), nothing)
         tspan = (1.0, 20.0)
 
         function g(p)
@@ -52,7 +52,7 @@ end
     prob = IntervalNonlinearProblem(quadratic_f, (1.0, 20.0), 2.0)
     ϵ = eps(Float64) # least possible tol for all methods
 
-    @testset for alg in (Bisection(), Falsi(), ITP())
+    @testset for alg in (Bisection(), Falsi(), ITP(), nothing)
         @testset for abstol in [0.1, 0.01, 0.001, 0.0001, 1e-5, 1e-6, 1e-7]
             sol = solve(prob, alg; abstol)
             result_tol = abs(sol.u - sqrt(2))
@@ -76,7 +76,7 @@ end
 end
 
 @testitem "Flipped Signs and Reversed Tspan" setup=[RootfindingTestSnippet] tags=[:core] begin
-    @testset for alg in (Alefeld(), Bisection(), Falsi(), Brent(), ITP(), Ridder())
+    @testset for alg in (Alefeld(), Bisection(), Falsi(), Brent(), ITP(), Ridder(), nothing)
         f1(u, p) = u * u - p
         f2(u, p) = p - u * u