An Adaptive Spectral Method for Oscillatory Second-Order Linear ODEs with Frequency-Independent Cost
From MaRDI portal
Publication:6154533
Abstract: We introduce an efficient numerical method for second order linear ODEs whose solution may vary between highly oscillatory and slowly changing over the solution interval. In oscillatory regions the solution is generated via a nonoscillatory phase function that obeys the nonlinear Riccati equation. We propose a defect-correction iteration that gives an asymptotic series for such a phase function; this is numerically approximated on a Chebyshev grid with a small number of nodes. For analytic coefficients we prove that each iteration, up to a certain maximum number, reduces the residual by a factor of order of the local frequency. The algorithm adapts both the step size and the choice of method, switching to a conventional spectral collocation method away from oscillatory regions. In numerical experiments we find that our proposal outperforms other state-of-the-art oscillatory solvers, most significantly at low-to-intermediate frequencies and at low tolerances, where it may use up to times fewer function evaluations. Even in high frequency regimes, our implementation is on average 10 times faster than other specialized solvers.
Recommendations
- On the numerical solution of second order ordinary differential equations in the high-frequency regime
- Phase function methods for second order inhomogeneous linear ordinary differential equations
- Some numerical algorithms for solving the highly oscillatory second-order initial value problems
- Adiabatic Filon-type methods for highly oscillatory second-order ordinary differential equations
- A fast and well-conditioned spectral method
Cites work
- scientific article; zbMATH DE number 3867259 (Why is no real title available?)
- scientific article; zbMATH DE number 5594165 (Why is no real title available?)
- scientific article; zbMATH DE number 486467 (Why is no real title available?)
- scientific article; zbMATH DE number 1069615 (Why is no real title available?)
- scientific article; zbMATH DE number 1889799 (Why is no real title available?)
- scientific article; zbMATH DE number 1409619 (Why is no real title available?)
- A fast and well-conditioned spectral method
- Adiabatic dynamics of one-dimensional classical Hamiltonian dissipative systems
- Applied mathematics
- Approximation theory and approximation practice
- Chebfun and numerical quadrature
- Implementing Clenshaw-Curtis quadrature, I methodology and experience
- Improved estimates for nonoscillatory phase functions
- Introduction to quantum mechanics
- Methods of nonlinear analysis. Vol. 1
- NIST handbook of mathematical functions
- Numerical methods. Design, analysis, and computer implementation of algorithms.
- On the evaluation of layer potentials close to their sources
- On the existence of nonoscillatory phase functions for second order ordinary differential equations in the high-frequency regime
- On the numerical solution of second order ordinary differential equations in the high-frequency regime
- Quasilinearization approach to computations with singular potentials
- Quasilinearization method and WKB
- Solving Ordinary Differential Equations I
- Spectral Methods in MATLAB
- The chebop system for automatic solution of differential equations
- The devil's invention: Asymptotic, superasymptotic and hyperasymptotic series
- The ultraspherical spectral element method
- Theory of the alternating-gradient synchrotron
- Unified QLM for regular and arbitrary singular potentials
- WKB-based scheme with adaptive step size control for the Schrödinger equation in the highly oscillatory regime
- WKB-based schemes for the oscillatory 1D Schrödinger equation in the semiclassical limit
Cited in
(2)
This page was built for publication: An Adaptive Spectral Method for Oscillatory Second-Order Linear ODEs with Frequency-Independent Cost
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6154533)
