This home web page for the project Clenshaw-Curtis rules for oscillatory with (possible) logarithmic singularities.

First, the papers where these rules were dealt:

- Domínguez, V.
*Filon-Clenshaw-Curtis rules for a class of highly-oscillatory integrals with logarithmic singularities*. Journal of Computational and Applied Mathematics , 261 (2014), pp 299–319. Preprint available in http://arxiv.org/abs/1305.1365
- Domínguez, V, Graham, I.G., and Kim, T.
*Filon-Clenshaw-Curtis rules for highly-oscillatory integrals with algebraic singularities and stationary points*. SIAM Journal of Numerical Analysis 51 (2013) no. 3, 1542–1566. Preprint available in http://arxiv.org/abs/1207.2283
- Domínguez, V., Graham, I., and Smyshlyaev, V.
*Stability and error estimates for Filon-Clenshaw-Curtis rules for highly-oscillatory integrals*. IMA Journal of Numerical Analysis 31 (2011), no 4, 1253-1280. Preprint available http://opus.bath.ac.uk/26655

The algorithms for the implementation are detailed and fully explained there.

In short, we deal with the approximation of

(Notice that in the second integral the interval is fixed). Here is real and . Actually and we have the Classical Clenshaw-Curtis rule in the first case and a product Clenshaw-Curtis rule in the second case.

These rules are based on replacing by the interpolating polynomial at the Chebyshev nodes. Hence, in the first case, and with

where

The second integral is computed *exactly. *This can be done, in theory. In practise the analytical expression are both complicate and unstable as increases. The implementation is based on the following lines:

- Compute $latex p_n$ in the Chebyshev basis , which can be done efficiently using FFT techniques
- Compute
via a fast recursive relations which are shown to be stable.

The other integral is dealt following similar ideas.

You can download an updated version, as of July 2017, of the Matlab implementation of such algorithms. Fejer-type quadrature rules are supported as well (not analysed in the cited papers, though)

I appreciate any comment, suggestion or (even much more) correction you can share with me.

### Like this:

Like Loading...