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].

Analytic solutions to nonlinear ODEs via spectral power series

A Functional Analytic Approach to the Power Series Solutions of an Nonlinear Differential Equations

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

A frequency-independent solver for systems of first order linear ordinary differential equations

Regularity, synthesis, rigidity and analytic classification for linear ordinary differential equations of second order

Solving Nonlinear Ordinary Differential Equations Using the Invariant Manifolds and Koopman Eigenfunctions

Solving nonlinear ODEs with the ultraspherical spectral method

A search for a spectral technique to solve nonlinear fractional differential equations

On Solutions of Systems of Differential Equations on Half-Line with Summable Coefficients

A power series method for solving ordinary and partial differentials equations motivated by domain growth

Spectral Methods for Solution of Differential and Functional Equations

An accelerated frequency-independent solver for oscillatory differential equations

On the use of asymptotically motivated gauge functions to obtain convergent series solutions to nonlinear ODEs

On the Obtaining Solutions of Nonlinear Differential Equations by Means of the Solutions of Simpler Linear or Nonlinear Differential Equations

Operational Solution of Non-Integer Ordinary and Evolution-Type Partial Differential Equations

A General Approach to Obtain Series Solutions of Nonlinear Differential Equations

Numerical solution of highly oscillatory ordinary differential equations

High-accuracy power series solutions with arbitrarily large radius of convergence for fractional nonlinear differential equations

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

The Method of Polynomial Particular Solutions for Solving Nonlinear Poisson-Type Equations

Convergent Power Series for Boundary Value Problems and Eigenproblems with Application to Atmospheric and Oceanic Tides