Abstract:The numerical evaluation of integrals of the form \begin{align*} \int_a^b f(x) e^{ikg(x)}\,dx \end{align*} is an important problem in scientific computing with significant applications in many branches of applied mathematics, science and engineering. The numerical approximation of such integrals using classical quadratures can be prohibitively expensive at high oscillation frequency ($k \gg 1$) as the number of quadrature points needed for achieving a reasonable accuracy must grow proportionally to $k$. To address this significant computational challenge, starting with Filon in 1930, several specialized quadratures have been developed to compute such oscillatory integrals efficiently. A crucial element in such Filon-type quadrature is the accurate evaluation of certain moments which poses a significant challenge when non-linear phase functions $g$ are involved. In this paper, we propose an equispaced-grid Filon-type quadrature for computing such highly oscillatory integrals that utilizes a Fourier extension of the slowly varying envelope $f$. This strategy is primarily aimed at significantly simplifying the moment calculations, even when the phase function has stationary points. Moreover, the proposed approach can also handle certain integrable singularities in the integrand. We analyze the scheme to theoretically establish high-order convergence rates. We also include a wide variety of numerical experiments, including oscillatory integrals with algebraic and logarithmic singularities, to demonstrate the performance of the quadrature.

Phenomenological Scale Factors for the Liquid-Vapor Critical Transition of Pure Fluids

A Fast Precise Algorithm for Duhamel Integrals

Improving the Accuracy of Exponentially Converging Quadratures

Quantitative estimates: How well does the discrete Fourier transform approximate the Fourier transform on $\mathbb{R}$

Criteria for the Application of Double Exponential Transformation

High-accuracy Formula for Discrete Calculation of Fourier Transforms

Some notes on the trapezoidal rule for Fourier type integrals

Computation of the exponential function of matrices by a formula without oscillatory integrals on infinite intervals

Approximating the Analytic Fourier Transform with the Discrete Fourier Transform

On Fourier re-expansions

An Efficient Quadrature Rule for the Oscillatory Infinite Generalized Bessel Transform with a General Oscillator and Its Error Analysis

On Fourier Approximation of Periodic Functions and Their Conjugates Belonging to the Weighted Lipschitz Class

Exponential sum approximations of finite completely monotonic functions

Towards non-linear quadrature formulae

Approximation of Integral Operators Using Product-Convolution Expansions

Exact Calculation of the Mean-Square Error in the Method of Expansion of Iterated Ito Stochastic integrals Based on Generalized Multiple Fourier Series

Evaluating a double integral using Euler's method and Richardson extrapolation

Parameter selection and numerical approximation properties of Fourier extensions from fixed data

Computation of highly oscillatory integrals using a Fourier extension approximation

Generalized convolution quadrature based on the trapezoidal rule