Abstract
Algorithms are described that make it possible to manipulate piecewise-smooth functions on real intervals numerically with close to machine precision. Break points are introduced in some such calculations at points determined by numerical root finding and in others by recursive subdivision or automatic edge detection. Functions are represented on each smooth subinterval by Chebyshev series or interpolants. The algorithms are implemented in object-oriented Matlab in an extension of the chebfun system, which was previously limited to smooth functions on [- 1, 1]. 2009. Published by Oxford University Press on behalf of the Institute of Mathematics and its Applications. All rights reserved.
Original language | English (US) |
---|---|
Pages (from-to) | 898-916 |
Number of pages | 19 |
Journal | IMA Journal of Numerical Analysis |
Volume | 30 |
Issue number | 4 |
DOIs | |
State | Published - Oct 2010 |
Externally published | Yes |
Keywords
- Chebyshev series
- barycentric interpolation
- chebfun system
- edge detection
ASJC Scopus subject areas
- Mathematics(all)
- Computational Mathematics
- Applied Mathematics