English

An accelerated Levin-Clenshaw-Curtis method for the evaluation of highly oscillatory integrals

Numerical Analysis 2025-03-13 v2 Numerical Analysis

Abstract

The efficient approximation of highly oscillatory integrals plays an important role in a wide range of applications. Whilst traditional quadrature becomes prohibitively expensive in the high-frequency regime, Levin methods provide a way to approximate these integrals in many settings at uniform cost. In this work, we present an accelerated version of Levin methods that can be applied to a wide range of physically important oscillatory integrals, by exploiting the banded action of certain differential operators on a Chebyshev polynomial basis. Our proposed version of the Levin method can be computed essentially in just O(νlogν)\mathcal{O}(\nu\log\nu) operations, where ν\nu is the number of quadrature points and the dependence of the cost on a number of additional parameters is made explicit in the manuscript. This presents a significant speed-up over the direct computation of the Levin method in current state-of-the-art. We outline the construction of this accelerated method for a fairly broad class of integrals and support our theoretical description with a number of illustrative numerical examples.

Keywords

Cite

@article{arxiv.2404.11448,
  title  = {An accelerated Levin-Clenshaw-Curtis method for the evaluation of highly oscillatory integrals},
  author = {Arieh Iserles and Georg Maierhofer},
  journal= {arXiv preprint arXiv:2404.11448},
  year   = {2025}
}