Abstract:In this paper, we combine the Smolyak technique for multi-dimensional interpolation with the Filon-Clenshaw-Curtis (FCC) rule for one-dimensional oscillatory integration, to obtain a new Filon-Clenshaw-Curtis-Smolyak (FCCS) rule for oscillatory integrals with linear phase over the d − d- dimensional cube [ − 1 , 1 ] d [-1,1]^d . By combining stability and convergence estimates for the FCC rule with error estimates for the Smolyak interpolation operator, we obtain an error estimate for the FCCS rule, consisting of the product of a Smolyak-type error estimate multiplied by a term that decreases with O ( k − d ~ ) \mathcal {O}(k^{-\widetilde {d}}) , where k k is the wavenumber and d ~ \widetilde {d} is the number of oscillatory dimensions. If all dimensions are oscillatory, a higher negative power of k k appears in the estimate. As an application, we consider the forward problem of uncertainty quantification (UQ) for a one-space-dimensional Helmholtz problem with wavenumber k k and a random heterogeneous refractive index, depending in an affine way on d d i.i.d. uniform random variables. After applying a classical hybrid numerical-asymptotic approximation, expectations of functionals of the solution of this problem can be formulated as a sum of oscillatory integrals over [ − 1 , 1 ] d [-1,1]^d , which we compute using the FCCS rule. We give numerical results for the FCCS rule which illustrate its theoretical properties and show that the accuracy of the UQ algorithm improves when both k k and the order of the FCCS rule increase. We also give results for both the quadrature and UQ problems when the underlying FCCS rule uses a dimension-adaptive Smolyak interpolation. These show increasing accuracy for the UQ problem as both the adaptive tolerance decreases and k k increases, requiring very modest increase in work as the stochastic dimension increases, for a case when the affine expansion in random variables decays quickly.

Is hyperinterpolation efficient in the approximation of singular and oscillatory functions?

A Novel Interpolation Scheme for Full-Wave Electromagnetic Simulations

Bypassing the quadrature exactness assumption of hyperinterpolation on the sphere

Hybrid hyperinterpolation over general regions

Quasi-interpolation for high-dimensional function approximation

Quadrature Methods for Highly Oscillatory Singular Integrals

Hard thresholding hyperinterpolation over general regions

An efficient iterative procedure in hyperbolic space and application to non-linear delay integral equation

Quasi-interpolators with application to postprocessing in finite element methods

A Filon-Clenshaw-Curtis-Smolyak rule for multi-dimensional oscillatory integrals with application to a UQ problem for the Helmholtz equation

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

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

Construction techniques for highly accurate quasi-interpolation operators

Computation of highly oscillatory integrals using a Fourier extension approximation

An Efficient and Accurate Interpolation Strategy for Multi-Dimensional Functions

Scaling Optimized Hermite Approximation Methods

Quasi-interpolation projectors for Subdivision Surfaces

Generalized Interpolation in H(Infinity) with A Complexity Constraint

Approximation to the K-Th Derivatives by Multiquadric Quasi-Interpolation Method

An Efficient Quadrature Rule for Highly Oscillatory Integrals with Airy Function