Moment preserving local spline projection operators

This article describes an elementary construction of a dual basis for non-uniform B-splines that is local, L ∞-stable and reproducts polynomials of any prescribed degree. This allows to define local projection operators with near-optimal approximation properties in any L q , 1 ≤ q ≤ ∞, and high order moment preserving properties. As the dual basis functions share the same piecewise polynomial structure as the underlying splines, simple quadrature formulas can be used to compute the projected spline coefficients