Abstract:We establish a new perturbation theory for orthogonal polynomials using a Riemann--Hilbert approach and consider applications in numerical linear algebra and random matrix theory. This new approach shows that the orthogonal polynomials with respect to two measures can be effectively compared using the difference of their Stieltjes transforms on a suitably chosen contour. Moreover, when two measures are close and satisfy some regularity conditions, we use the theta functions of a hyperelliptic Riemann surface to derive explicit and accurate expansion formulae for the perturbed orthogonal polynomials. In contrast to other approaches, a key strength of the methodology is that estimates can remain valid as the degree of the polynomial grows. The results are applied to analyze several numerical algorithms from linear algebra, including the Lanczos tridiagonalization procedure, the Cholesky factorization and the conjugate gradient algorithm. As a case study, we investigate these algorithms applied to a general spiked sample covariance matrix model by considering the eigenvector empirical spectral distribution and its limits. For the first time, we give precise estimates on the output of the algorithms, applied to this wide class of random matrices, as the number of iterations diverges. In this setting, beyond the first order expansion, we also derive a new mesoscopic central limit theorem for the associated orthogonal polynomials and other quantities relevant to numerical algorithms.

Log Orthogonal Functions: Approximation Properties and Applications

An Exponential Convergence Approximation to Singularly Perturbed Problems by Log Orthogonal Functions

A Spectrally Accurate Approximation to Subdiffusion Equations Using the Log Orthogonal Functions

Approximation by a new sequence of operators involving Laguerre polynomials

Approximation to Distribution of Product of Random Variables Using Orthogonal Polynomials for Lognormal Density

Piecewise General Orthogonal Polynomial Operator and Its Properties

Approximation on Variable Exponent Spaces by Linear Integral Operators

New error bounds for Legendre approximations of differentiable functions

Generalized log orthogonal functions spectral collocation method for two dimensional weakly singular Volterra integral equations of the second kind

Some study of the approximate orthogonality connected with integral orthogonalities

Conservative polynomial approximations and applications to Fokker-Planck equations

Connecting Exceptional Orthogonal Polynomials of Different Kind

Differential orthogonality: Laguerre and Hermite cases with applications

Optimal approximants and orthogonal polynomials in several variables

Piecewise linear approximate solution of fractional order non-stiff and stiff differential-algebraic equations by orthogonal hybrid functions

A Riemann--Hilbert approach to the perturbation theory for orthogonal polynomials: Applications to numerical linear algebra and random matrix theory

Orthogonal Functions Based Integral Inequalities and Their Applications to Time Delay Systems

Zeros of generalized hypergeometric polynomials via finite free convolution. Applications to multiple orthogonality

Discrete orthogonality relations for multi-indexed Laguerre and Jacobi polynomials

Orthogonality of a new family of $q$-Sobolev type polynomials

Bivariate second--order linear partial differential equations and orthogonal polynomial solutions