Piecewise smooth chebfuns
Ricardo Pachόn‚ Rodrigo Platte and Nick Trefethen
Abstract
Algorithms are described that make it possible to manipulate piecewise-smooth functions on real intervals numerically with close to machine precision. Breakpoints are introduced in some such calculations at points determined by numerical rootfinding, 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].
Institution
Oxford University Computing Laboratory
Month
May
Number
NA−08/07
Year
2008