For such a quasi-interpolation operator $L$ which reproduces $\Pi_m(\Rd)$ and any $r\ge0$, we give an explicit construction of a quasi-interpolant $R_m^{r+m}L=L+A$ which reproduces $\Pi_{m+r}$, together with an integral error formula which involves only the $(m+r+1)$--st derivative of the function approximated. The operator $R_m^{m+r}L$ is defined on functions with $r$ additional orders of smoothness than those on which $L$ is defined. This very general construction holds in all dimensions $d$. A number of representative examples are considered.
Keywords: Quasi-interpolation, Lagrange interpolation, Bernstein polynomial, finite element method, multivariate polynomial approximation, error formula, multipoint Taylor formula, divided differences, Chu--Vandermonde convolution
Math Review Classification: Primary primary 41A80, 65D05; Secondary 41A05, 41A10.
Length: 12 pages
Last Updated: 29 January 2008