The complexity of the solutions of a differential equation of the form $y''(t) + \omega^2 q(t) y(t) =0$ depends not only on that of the coefficient $q$, but also on the magnitude of the parameter $\omega$. In the most widely-studied case, when $q$ is positive, the solutions of equations of this type oscillate at a frequency that increases linearly with $\omega$ and standard ODE solvers require $\mathcal{O}\left(\omega\right)$ time to calculate them. It is well known, though, that phase function methods can be used to solve such equations numerically in time independent of $\omega$. Unfortunately, the running time of these methods increases with $\omega$ when they are applied in the commonly-occurring case in when $q$ has zeros in the solution domain (i.e., when the differential equation has turning points). Here, we introduce a generalized phase function method adapted to equations with simple turning points. More explicitly, we show the existence of slowly-varying ``Airy phase functions'' which represent the solutions of such equations at a cost which is independent of $\omega$ and describe a numerical method for calculating these Airy phase functions in time independent of $\omega$. Using our method, initial or boundary value problems for a large class of second order linear ordinary differential equations with turning points whose coefficients depend on a parameter $\omega$ can be solved in time independent of $\omega$. We also give the results of numerical experiments conducted to demonstrate the properties of our method, including one in which we used our algorithm to rapidly calculate associated Legendre functions of a wide range of orders and degrees.
翻译:暂无翻译