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

A rapid solution of a kind of 1D Fredholm oscillatory integral equation

A Universal Solution to One-Dimensional Oscillatory Integrals.

An Improved Levin Quadrature Method for Highly Oscillatory Integrals

On an Improved-Levin Oscillatory Quadrature Method

Numerical Integration of Highly Oscillatory Functions with and without Stationary Points

A Chebyshev Collocation Method for a Class of Fredholm Integral Equations with Highly Oscillatory Kernels.

Numerical integration of ordinary differential equations with rapidly oscillatory factors

Numerical solutions to Volterra integral equations of the second kind with oscillatory trigonometric kernels

Fast Evaluation of Scattering Oscillatory Integrals over Fan-Shaped Regions

An Overview of Numerical Methods for the First Kind Fredholm Integral Equation

Numerical analysis of a fast integration method for highly oscillatory functions

A Fast Numerical Solution Method for Two Dimensional Fredholm Integral Equations of the Second Kind

Asymptotics on the Fredholm Integral Equation with a Highly Oscillatory and Weakly Singular Kernel

Efficient integration for a class of highly oscillatory integrals.

Approximate solution of linear ordinary differential equations with variable coefficients

On Quadrature of Highly-oscillatory Functions

Efficient numerical simulation algorithm for highly oscillatory nonhomogeneous dynamic systems

An accelerated frequency-independent solver for oscillatory differential equations

Deep Neural Network Solutions for Oscillatory Fredholm Integral Equations

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

A Collocation Numerical Method for Highly Oscillatory Algebraic Singular Volterra Integral Equations