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

Arbitrary-Order Trigonometric Fourier Collocation Methods for Multi-Frequency Oscillatory Systems

Trigonometric Collocation Methods Based on Lagrange Basis Polynomials for Multi-Frequency Oscillatory Second-Order Differential Equations.

Exponential Fourier Collocation Methods for Solving First-Order Differential Equations

An accelerated frequency-independent solver for oscillatory differential equations

Arbitrarily High-Order Time-Stepping Schemes Based on the Operator Spectrum Theory for High-Dimensional Nonlinear Klein-Gordon Equations

Functionally-fitted energy-preserving methods for solving oscillatory nonlinear Hamiltonian systems

A Jacobi spectral collocation method for solving multi-dimensional nonlinear fractional sub-diffusion equations

High-Order Symplectic and Symmetric Composition Methods for Multi-Frequency and Multi-Dimensional Oscillatory Hamiltonian Systems

A discrete least squares collocation method for two-dimensional nonlinear time-dependent partial differential equations

A semi-analytical collocation method for solving multi-term variable-order time fractional partial differential equations

Exponential Collocation Methods for Conservative or Dissipative Systems.

Explicit Symplectic Multidimensional Exponential Fitting Modified Runge-Kutta-Nyström Methods

An Extended Discrete Gradient Formula for Oscillatory Hamiltonian Systems

Efficient implementation of RKN-type Fourier collocation methods for second-order differential equations

Numerical methods for highly oscillatory Volterra integral equations with general oscillators

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

Improved Filon-type asymptotic methods for highly oscillatory differential equations with multiple time scales

Effective Integrators for Nonlinear Second-Order Oscillatory Systems with a Time-Dependent Frequency Matrix

Exponential Collocation Methods Based on Continuous Finite Element Approximations for Efficiently Solving the Cubic Schrödinger Equation

Direct/split invariant-preserving Fourier pseudo-spectral methods for the rotation-two-component Camassa–Holm system with peakon solitons

Piecewise orthogonal collocation for computing periodic solutions of coupled delay equations