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.

Gauss-Type Quadrature for Highly Oscillatory Integrals with Algebraic Singularities and Applications

Gauss quadrature for integrals and sums

Gaussian quadrature rules for composite highly oscillatory integrals

On Lebesgue Integral Quadrature

Spherical configurations and quadrature methods for integral equations of the second kind

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

Quadrature methods for singular integral equations of Mellin type based on the zeros of classical Jacobi polynomials, II

Computation of highly oscillatory integrals using a Fourier extension approximation

An Efficient Quadrature Rule for Highly Oscillatory Integrals with Airy Function

On quadrature for singular integral operators with complex symmetric quadratic forms

Singularity swapping method for nearly singular integrals based on trapezoidal rule

On the Convolution Quadrature Rule for Integral Transforms with Oscillatory Bessel Kernels

Quadrature-free polytopic discontinuous Galerkin methods for transport problems

On the Numerical Evaluation of a Class of Oscillatory Integrals in Worldline Variational Calculations

Implicitly Linear Jacobi Spectral-Collocation Methods for Two-Dimensional Weakly Singular Volterra-Hammerstein Integral Equations

Rigorous numerical integration of algebraic functions

Numerical evaluation of oscillatory integrals via automated steepest descent contour deformation

On quadrature of highly oscillatory Bessel function via asymptotic analysis of simplex integrals