On the abs-polynomial expansion of piecewise smooth functions

  • Tom Streubel has observed that for functions in abs-normal form, generalized Taylor expansions of arbitrary order $\bar d-1$ can be generated by algorithmic piecewise differentiation. Abs-normal form means that the real or vector valued function is defined by an evaluation procedure that involves the absolute value function $|...|$ apart from arithmetic operations and $\bar d$ times continuously differentiable univariate intrinsic functions. The additive terms in Streubel's expansion are abs-polynomial, i.e. involve neither divisions nor intrinsics. When and where no absolute values occur, Moore's recurrences can be used to propagate univariate Taylor polynomials through the evaluation procedure with a computational effort of $\mathcal O({\bar d}^2)$, provided all univariate intrinsics are defined as solutions of linear ODEs. This regularity assumption holds for all standard intrinsics, but for irregular elementaries one has to resort to Faa di Bruno's formula, which has exponential complexity in $\bar d$. As already conjectured we show that the Moore recurrences can be adapted for regular intrinsics to the abs-normal case. Finally, we observe that where the intrinsics are real analytic the expansions can be extended to infinite series that converge absolutely on spherical domains.

This document is embargoed until:


Export metadata

Additional Services

Share in Twitter Search Google Scholar
Author:Andreas GriewankORCiD, Tom StreubelORCiD, Caren TischendorfORCiD
Document Type:ZIB-Report
Tag:Moore recurrences; Nonsmooth Taylor polynomial/series; abs-linear form; abs-normal form; absolute convergence; forward mode propagation; quadratic complexity
Date of first Publication:2020/06/05
Series (Serial Number):ZIB-Report (20-14)