From 6bd6697ca9f79cdab6d21d01fa78bd9ccf65038c Mon Sep 17 00:00:00 2001
From: oscarddssmith <oscar.smith@juliacomputing.com>
Date: Wed, 19 Mar 2025 11:16:19 -0400
Subject: [PATCH] minor itp speedups

---
 lib/BracketingNonlinearSolve/src/itp.jl | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/lib/BracketingNonlinearSolve/src/itp.jl b/lib/BracketingNonlinearSolve/src/itp.jl
index 50443e2e9..153981faf 100644
--- a/lib/BracketingNonlinearSolve/src/itp.jl
+++ b/lib/BracketingNonlinearSolve/src/itp.jl
@@ -95,15 +95,15 @@ function CommonSolve.solve(
     k2 = alg.k2
     k1 = alg.scaled_k1 * abs(right - left)^(1 - k2)
     n0 = alg.n0
-    n_h = ceil(log2(abs(right - left) / (2 * ϵ)))
     mid = (left + right) / 2
     x_f = left + (right - left) * (fl / (fl - fr))
     xt = left
     xp = left
     r = zero(left) # minmax radius
     δ = zero(left) # truncation error
-    σ = 1.0
-    ϵ_s = ϵ * 2^(n_h + n0)
+    σ = one(mid)
+    n_h = exponent(abs(right - left) / (2 * ϵ))
+    ϵ_s = ϵ * exp2(n_h + n0)
 
     i = 1
     while i ≤ maxiters