Convexity preserving interpolatory subdivision with conic precision

The paper is concerned with the problem of shape preserving interpolatory subdivision. For arbitrarily spaced, planar input data an efficient non-linear subdivision algorithm reproducing conic sections and respecting the convexity properties of the initial data, is here presented. Significant numerical examples are included to illustrate the effectiveness of the proposed method and the smoothness of the limit curves.