Let $ E \subset \mathbb{ R}^2 $ be a finite set, and let $ f : E \to [0,\infty) $. In this paper, we address the algorithmic aspects of nonnegative $C^2$ interpolation in the plane. Specifically, we provide an efficient algorithm to compute a nonnegative $C^2(\mathbb{ R}^2)$ extension of $ f $ with norm within a universal constant factor of the least possible. We also provide an efficient algorithm to approximate the trace norm.
Categories: math.CA math.OC