An algorithm for convexity-preserving surface interpolation

A simple and efficient method of convexity-preserving interpolation for grid data is introduced. An earlier algorithm due to Roulier (1987), based on a shape-preserving curve interpolation scheme due to McAllister and Roulier (1981), is modified to use a simpler but equivalent curve interpolation scheme described in Iqbal (1992). Numerical examples are provided to test the performance of the method with the slopes that are further improved using iterative technique to produce more visually pleasing surfaces.