Exponentially convergent Fourier-Chebyshev quadrature schemes on bounded and infinite intervals (Q1116643): Difference between revisions
From MaRDI portal
Revision as of 14:26, 19 June 2024
scientific article
Language | Label | Description | Also known as |
---|---|---|---|
English | Exponentially convergent Fourier-Chebyshev quadrature schemes on bounded and infinite intervals |
scientific article |
Statements
Exponentially convergent Fourier-Chebyshev quadrature schemes on bounded and infinite intervals (English)
0 references
1987
0 references
The method of \textit{C. W. Clenshaw} and \textit{A. R. Curtis} [Numer. Math. 2, 197-205 (1960; Zbl 0093.140)] for numerical integration is extended to semiinfinite ([0,\(\infty]\) and infinite [-\(\infty,\infty]\) intervals. The common framework for both these extensions and for integration on a finite interval is to (1) map the integration domain to \(t\in [0,\pi]\), (2) compute a Fourier sine or cosine approximation to the transformed integrand via interpolation, and (3) integrate the approximation. The interpolation is most easily performed via the sine or cosine cardinal functions, which are discussed in the appendix. The algorithm is mathematically equivalent to expanding the integrand in (mapped or unmapped) Chebyshev polynomials as done by Clenshaw and Curtis, but the trigonometric approach simplifies the mechanics. Like Gaussian quadratures, the error for the change-of-coordinates Fourier method decreases exponentially with N, the number of grid points, but the generalized Curtis-Clenshaw algorithm is much easier to program than Gaussian quadrature because the abscissas and weights are given by simple, explicit formulas.
0 references
Fourier approximation
0 references
infinite intervals
0 references
rational Chebyshev function
0 references
adaptive quadrature
0 references
Gaussian quadratures
0 references
Fourier method
0 references
Curtis- Clenshaw algorithm
0 references
0 references