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

Reproducing Kernel-Based Piecewise Methods for Efficiently Solving Oscillatory Systems of Second-Order Initial Value Problems

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

A new application of the reproducing kernel method for solving linear systems of fractional order Volterra integro-differential equations

Reproducing Kernel Function-Based Filon and Levin Methods for Solving Highly Oscillatory Integral

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

A Sixth Order Kernel Functions Approach for Nonlinear Fourth Order Boundary Value Problems

A trigonometrically fitted explicit Numerov-type method for second-order initial value problems with oscillating solutions

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

A Novel Kernel Functions Algorithm for Solving Impulsive Boundary Value Problems.

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

Efficient Multistep Algorithms for First-Order IVPs with Oscillating Solutions: II Implicit and Predictor–Corrector Algorithms

Spectral-Galerkin Method for Second Kind VIEs with Highly Oscillatory Kernels of The Stationary Point

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

An efficient off-step spline algorithm for wave simulation of nonlinear Kuramoto–Sivashinsky and Korteweg-de Vries equations

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

An accelerated frequency-independent solver for oscillatory differential equations

Cost-reduction implicit exponential Runge-Kutta methods for highly oscillatory systems

Extended Phase Properties and Stability Analysis of RKN-type Integrators for Solving General Oscillatory Second-Order Initial Value Problems

An efficient numerical algorithm for solving nonlinear fractional Volterra integro-differential equation

High asymptotic order numerical methods for highly oscillatory ODEs with large initial data

Convergence Order of the Reproducing Kernel Method for Solving Boundary Value Problems