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

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

Exponential Fourier Collocation Methods for Solving First-Order Differential Equations

Semi-analytical Exponential RKN Integrators for Efficiently Solving High-Dimensional Nonlinear Wave Equations Based on FFT Techniques.

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

Symplecticity-preserving continuous-stage Runge-Kutta-Nyström methods

An Efficient Numerical Method for Solving a Class of Nonlinear Fractional Differential Equations and Error Estimates

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

New efficient four-stage implicit trigonometrically fitted modified RKN for second-order ODEs

Novel exponentially fitted two-derivative Runge–Kutta methods for solving the radial Schrödinger equation

Solve Partial Differential Equations by Two or More Radial Basis Functions

An accelerated frequency-independent solver for oscillatory differential equations

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

An Exponential Time Differencing Runge–Kutta Method ETDRK32 for Phase Field Models

A novel fourth-order scheme for two-dimensional Riesz space fractional nonlinear reaction-diffusion equations and its optimal preconditioned solver

Runge–Kutta methods for linear ordinary differential equations

Nonlinear Stability and Convergence of ERKN Integrators for Solving Nonlinear Multi-Frequency Highly Oscillatory Second-Order ODEs with Applications to Semi-Linear Wave Equations

Parallel-in-time Multilevel Krylov Methods: A Prototype

An innovative fourth-order numerical scheme with error analysis for Lane-Emden-Fowler type systems

Two new compact finite-difference schemes for the solution of boundary value problems in second-order non-linear ordinary differential equations, using non-uniform grids

Two fast and unconditionally stable finite difference methods for Riesz fractional diffusion equations with variable coefficients

A Fractional Spectral Collocation for Solving Second Kind Nonlinear Volterra Integral Equations with Weakly Singular Kernels