Filon-Clenshaw-Curtis rules for highly oscillatory integrals with algebraic singularities and stationary points

From MaRDI portal
Publication:2845602




Abstract: In this paper we propose and analyse composite Filon-Clenshaw-Curtis quadrature rules for integrals of the form Ik[a,b](f,g):=intabf(x)exp(mathrmikg(x))dx, where kgeq0, f may have integrable singularities and g may have stationary points. Our composite rule is defined on a mesh with M subintervals and requires MN+1 evaluations of f. It satisfies an error estimate of the form CNkrMN1+r, where r is determined by the strength of any singularity in f and the order of any stationary points in g and CN is a constant which is independent of k and M, but depends on N. The regularity requirements on f and g are explicit in the error estimates. For fixed k, the rate of convergence of the rule as Mightarrowinfty is the same as would be obtained if f was smooth. Moreover, the quadrature error decays at least as fast as kightarrowinfty as does the original integral Ik[a,b](f,g). For the case of nonlinear oscillators g, the algorithm requires the evaluation of g1 at non-stationary points. Numerical results demonstrate the sharpness of the theory. An application to the implementation of boundary integral methods for the high-frequency Helmholtz equation is given.




Cited in
(54)






This page was built for publication: Filon-Clenshaw-Curtis rules for highly oscillatory integrals with algebraic singularities and stationary points

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2845602)