Piecewise-smooth chebfuns

摘要

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].