A numerical study of the Legendre-Galerkin method for the evaluation of the prolate spheroidal wave functions

• Published in: Numerical algorithms 2015-04, Vol.68 (4), p.691-710
• Main Authors: Schmutzhard, Sebastian, Hrycak, Tomasz, Feichtinger, Hans G.
• Format: Article
• Language: English
• Subjects: Algebra; Algorithms; Approximation; Boundary conditions; Computer Science; Constants; Continuity; Differential equations; Eigenvectors; Estimates; Galerkin method; Mathematical models; Numeric Computing; Numerical Analysis; Operators; Operators (mathematics); Original Paper; Prolate spheroids; Theory of Computation; Upper bounds; Wave functions
• ISSN: 1017-1398; 1572-9265
• DOI: 10.1007/s11075-014-9867-3

We conduct a numerical study of the Legendre-Galerkin method for the evaluation of the prolate spheroidal wave functions (PSWFs) viewed as the eigenfunctions of the prolate differential operator with boundary conditions of continuity. Our experiments indicate that the minimal dimension N of the Legendre-Galerkin matrix for the evaluation of the n th prolate with precision ε is O n + nc , as n , c → ∞ , where c >0 is the bandwidth parameter. The behavior of N , when either c or n is held constant, is also examined. As a consequence, we obtain an upper bound on the complexity of the evaluation of the prolates. We also study the condition number of the approximate Legendre coefficients, computed as an eigenvector of the Legendre-Galerkin matrix. We observe experimentally that for fixed precision ε , an error estimate based on this condition number is O n + c as n , c → ∞ . We conclude that the Legendre-Galerkin method is accurate for fairly large values of n and c .