A Riemann–Hilbert approach to Jacobi operators and Gaussian quadrature. arXiv:1311.5838

The computation of the entries of Jacobi operators associated with orthogonal poly-nomials has important applications in numerical analysis. From truncating the oper-ator to form a Jacobi matrix, one can apply the Golub–Welsh algorithm to compute the Gaussian quadrature weights and nodes. Furthermore, the entries of the Jacobi operator are the coefficients in the three-term recurrence relationship for the polyno-mials. This provides an efficient method for evaluating the orthogonal polynomials. Here, we present an Ø(N) method to compute the first N rows of Jacobi operators from the associated weight. The method exploits the Riemann–Hilbert representation of the polynomials by solving a deformed Riemann–Hilbert problem numerically. We further adapt this computational approach to certain entire weights that are beyond the reach of current asymptotic Riemann–Hilbert techniques.