Exponential node clustering at singularities for rational approximation, quadrature, and PDEs.

Rational approximations of functions with singularities can converge at a root-exponential rate if the poles are exponentially clustered. We begin by reviewing this effect in minimax, least-squares, and AAA approximations on intervals and complex domains, conformal mapping, and the numerical solution of Laplace, Helmholtz, and biharmonic equations by the "lightning" method. Extensive and wide-ranging numerical experiments are involved.

New Laplace and Helmholtz solvers.

Numerical algorithms based on rational functions are introduced that solve the Laplace and Helmholtz equations on 2D domains with corners quickly and accurately, despite the corner singularities.

Continuous analogues of matrix factorizations.

Analogues of singular value decomposition (SVD), QR, LU and Cholesky factorizations are presented for problems in which the usual discrete matrix is replaced by a 'quasimatrix', continuous in one dimension, or a 'cmatrix', continuous in both dimensions. Two challenges arise: the generalization of the notions of triangular structure and row and column pivoting to continuous variables (required in all cases except the SVD, and far from obvious), and the convergence of the infinite series that define the cmatrix factorizations. Our generalizations of triangularity and pivoting are based on a new notion of a 'triangular quasimatrix'.