Optimal Area-Sensitive Bounds for Polytope Approximation
Sunil Arya, Guilherme D. da Fonseca, David M. Mount
TL;DR
This work introduces an area-sensitive bound for polytope approximation of convex bodies in fixed dimension: for a body K of minimum width at least ε, there exists an outer ε-approximation polytope with O((arad(K)/ε)^{(d-1)/2}) facets. The authors convert the problem to convex-function approximation on lower-dimensional domains, employing a functional view with a Legendre dual, and leverage Macbeath regions and Mahler-volume duality to construct small hitting sets. The approach yields tight area-based bounds and a robust framework linking caps, dual caps, and polar structures, with extensions to nonuniform and cap-cover results. The work also clarifies the role of area radius and related geometric measures in controlling approximation complexity, and connects to broader themes in convex geometry and approximation theory.
Abstract
Approximating convex bodies is a fundamental problem in geometry. Given a convex body $K$ in $\mathbb{R}^d$ for a fixed dimension $d$, the objective is to minimize the number of facets of an approximating polytope for a given Hausdorff error $\varepsilon$. The best known uniform bound, due to Dudley (1974), shows that $O((\text{diam}(K)/\varepsilon)^{(d-1)/2})$ facets suffice. Although this bound is optimal for fat objects, such as Euclidean balls, it is far from optimal for ``skinny'' convex bodies. Skinniness can be characterized relative to the Euclidean ball. Given a convex body $K$, define its area radius, $\text{arad}(K)$, to be the radius of the Euclidean ball having the same surface area as $K$. It follows from generalizations of the isoperimetric inequality that $\text{diam}(K) \geq 2 \cdot \text{arad}(K)$. We show that, given a convex body whose minimum width is at least $\varepsilon$, it is possible to approximate the body by a polytope having $O((\text{arad}(K)/\varepsilon)^{(d-1)/2})$ facets. Our approach works by first reducing the problem of approximating convex bodies to that of approximating convex functions. We employ a classical concept from convexity, called Macbeath regions. We demonstrate that there is a polar relationship between the Macbeath regions of a function and the Macbeath regions of its Legendre dual. This is combined with known bounds on the Mahler volume to bound the total size of the approximation.
