omit pointless extra Newton iteration for λ (#112)

JuliaMath · Jul 27, 2024 · 6ed7307 · 6ed7307
1 parent 2e6fc15
commit 6ed7307
Showing 1 changed file with 0 additions and 3 deletions.
diff --git a/src/gausskronrod.jl b/src/gausskronrod.jl
@@ -139,9 +139,6 @@ function eignewt(H::AbstractSymTri{T}, n::Integer) where {T<:Real}
             if isfinite(δλ)
                 lambda[i] -= δλ
                 if abs(δλ) ≤ 10 * eps(lambda[i])
-                    # do one final Newton iteration for luck and profit:
-                    δλ = eigpolyrat(H,lambda[i])
-                    lambda[i] -= δλ
                     break
                 end
             else