On the Numerical Calculation of the Roots of Special Functions Satisfying Second Order Ordinary Differential Equations

James Bremer
DOI: https://doi.org/10.1137/16m1057139
IF: 2.968
2017-01-01
SIAM Journal on Scientific Computing
Abstract:We describe a method for calculating the roots of special functions satisfying second order ordinary differential equations. It exploits the recent observation that the solutions of equations of this type can be represented via nonoscillatory phase functions, even in the high-frequency regime. Our approach requires $\mathcal{O}(1)$ operations per root and achieves near machine precision accuracy. Moreover, despite its great generality, our approach is competitive with (and in many cases, faster than) specialized, state-of-the-art methods for the construction of Gaussian quadrature rules of large orders when it is used in such a capacity. The performance of the scheme is illustrated with several numerical experiments.
mathematics, applied
What problem does this paper attempt to address?