Stable polefinding and rational least-squares fitting via eigenvalues.

A common way of finding the poles of a meromorphic function in a domain, where an explicit expression of is unknown but can be evaluated at any given , is to interpolate by a rational function such that at prescribed sample points , and then find the roots of . This is a two-step process and the type of the rational interpolant needs to be specified by the user. Many other algorithms for polefinding and rational interpolation (or least-squares fitting) have been proposed, but their numerical stability has remained largely unexplored. In this work we describe an algorithm with the following three features: (1) it automatically finds an appropriate type for a rational approximant, thereby allowing the user to input just the function , (2) it finds the poles via a generalized eigenvalue problem of matrices constructed directly from the sampled values in a one-step fashion, and (3) it computes rational approximants in a numerically stable manner, in that with small at the sample points, making it the first rational interpolation (or approximation) algorithm with guaranteed numerical stability. Our algorithm executes an implicit change of polynomial basis by the QR factorization, and allows for oversampling combined with least-squares fitting. Through experiments we illustrate the resulting accuracy and stability, which can significantly outperform existing algorithms.