Skip to content
Merged
Prev Previous commit
Next Next commit
tweak
  • Loading branch information
stevengj committed Feb 28, 2026
commit a4270193187b537d4e4bd244d6056972b90086bc
2 changes: 1 addition & 1 deletion src/roots.jl
Original file line number Diff line number Diff line change
Expand Up @@ -73,7 +73,7 @@ function roots(c::ChebPoly{1,<:Real}; tol::Real=5*epsvals(c.coefs), maxsize::Int
return filter_roots(c, λ)
else
# roughly halve the domain, constructing new Chebyshev polynomials on each half, and
# call roots recursively. as for chebfun, we split at an arbitrary point 0.004849834917525
# call roots recursively. Following chebfun, we split at an arbitrary point 0.004849834917525
# on [-1,1] rather than at 0 to avoid introducing additional spurious roots (since 0 is
# often a special point by symmetry).
split = oftype(tol, 1.004849834917525) * ((c.ub[1] - c.lb[1])/2) + c.lb[1]
Expand Down