Numerical algorithms based on analytic function values at roots of unity

Report

Let $f(z)$ be an analytic or meromorphic function in the closed unit disk sampled at the $n$th roots of unity. Based on these data, how can we approximately evaluate $f(z)$ or $f^{(m)}(z)$ at a point $z$ in the disk? How can we calculate the zeros or poles of $f$ in the disk? These questions exhibit in the purest form certain algorithmic issues that arise across computational science in areas including integral equations, partial differential equations, and large-scale linear algebra. We analyze the possibilities and emphasize a basic distinction between algorithms based on polynomial or rational interpolation and those based on trapezoidal rule approximations of Cauchy integrals. We then show how these developments apply to the problem of computing the eigenvalues in the disk of a matrix of large dimension.

Authors: Austin, A; Kravanja, P; Trefethen, L

• Publisher: SINUM
• Publication date: 2013-07-01