From 6bac7196007a888646542cce81e442525f76a148 Mon Sep 17 00:00:00 2001
From: Kipton Barros <kbarros@gmail.com>
Date: Wed, 17 Dec 2025 09:37:39 -0700
Subject: [PATCH] Generalized eigenvalue solver at higher precision

---
 src/barycentric.jl | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/barycentric.jl b/src/barycentric.jl
index a23afee7..52d156fe 100644
--- a/src/barycentric.jl
+++ b/src/barycentric.jl
@@ -212,8 +212,8 @@ function poles(r::Barycentric{T,S}) where {T,S}
         pol = filter( isfinite, eigvals(E, B) )
     catch
         # generalized eigen not available in extended precision, so:
-        λ = filter( z->abs(z)>1e-13, eigvals(E\B) )
-        pol = 1 ./ λ
+        (; α, β) = schur(complex(E), complex(B))
+        pol = filter( isfinite, α ./ β )
     end
     return pol
 end