tweak mu to fix JuliaApproximation#30

jamiecjx · Sep 8, 2023 · 897a052 · 897a052
1 parent dbc28d2
commit 897a052
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/recurrence.jl b/src/recurrence.jl
@@ -118,7 +118,7 @@ function backwardrecurrence!(K, A, B, C, z, nN::AbstractUnitRange, j...)
             μ = min(abs(data[k,j...]/data[k-1,j...]), abs(data[k-1,j...]/data[k-2,j...]))
             # data[k] * μ^M ≤ ε
             #  M ≥ log(ε/data[k])/log(μ)
-            N = ceil(Int, max(2N, min(maxiterations, log(eps(real(T))/100)/log(μ))))
+            N = ceil(Int, max(2N, min(maxiterations, log(eps(real(T))/100)/log(μ+eps(real(T))))))
             _growdata!(K, N, j...)
             resize!(u, N)
             data = K.data