Convexity preserving scattered data interpolation using Powell–Sabin elements

The problem of strict convexity preserving interpolation by a piecewise quadratic smooth function is solved. A constructive procedure is provided.