Abstract:Oscillatory second order linear ordinary differential equations arise in many scientific calculations. Because the running times of standard solvers increase linearly with frequency when they are applied to such problems, a variety of specialized methods, most of them quite complicated, have been proposed. Here, we point out that one of the simplest approaches not only works, but yields a scheme for solving oscillatory second order linear ordinary differential equations which is significantly faster than current state-of-the-art techniques. Our method, which operates by constructing a slowly varying phase function representing a basis of solutions of the differential equation, runs in time independent of the frequency and can be applied to second order equations whose solutions are oscillatory in some regions and slowly varying in others. In the high-frequency regime, our algorithm discretizes the nonlinear Riccati equation satisfied by the derivative of the phase function via a Chebyshev spectral collocation method and applies the Newton-Kantorovich method to the resulting system of nonlinear algebraic equations. We prove that the iterates converge quadratically to a nonoscillatory solution of the Riccati equation. The quadratic convergence of the Newton-Kantorovich method and the simple form of the linearized equations ensure that this procedure is extremely efficient. Our algorithm then extends the slowly varying phase function calculated in the high-frequency regime throughout the solution domain by solving a certain third order linear ordinary differential equation related to the Riccati equation. We describe the results of numerical experiments showing that our algorithm is orders of magnitude faster than existing schemes, including the modified Magnus method [18], the current state-of-the-art approach [7] and the recently introduced ARDC method [1].

An Hermite–Obreshkov method for 2nd-order linear initial-value problems for ODE

An efficient ODE-solving method based on heuristic and statistical computations: αII-(2 + 3)P method

Numerical Solution of Higher-Order Linear and Nonlinear Ordinary Differential Equations with Orthogonal Rational Legendre Functions

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

Explicit numerical methods for solving singular initial value problems for systems of second-order nonlinear ODEs

An accelerated frequency-independent solver for oscillatory differential equations

A Matrix Framework for the Solution of ODEs: Initial-, Boundary-, and Inner-Value Problems

An Adaptive Spectral Method for Oscillatory Second-Order Linear ODEs with Frequency-Independent Cost

A New Methodology for the Development of Efficient Multistep Methods for First-Order IVPs with Oscillating Solutions

An analytical approximation scheme to two point boundary value problems of ordinary differential equations

Parallel-in-time high-order multiderivative IMEX solvers

A Parametric Six-Step Method for Second-Order IVPs with Oscillating Solutions

A multistep algorithm for ODEs

A solver for linear scalar ordinary differential equations whose running time is bounded independent of frequency

The explicit solution of linear, dissipative, second-order initial-boundary value problems with variable coefficients

Runge–Kutta methods for linear ordinary differential equations

New efficient numerical methods for some systems of linear ordinary differential equations

On the computation of accurate initial conditions for linear higher-index differential-algebraic equations and its application in initial value solvers

Solving nonlinear ODEs with the ultraspherical spectral method

Symmetric and Arbitrarily High-Order Birkhoff-Hermite Time Integrators and Their Long-Time Behaviour for Solving Nonlinear Klein-Gordon Equations

Numerical solution for high-order ordinary differential equations using H-ELM algorithm