Approximate Vertex Enumeration

TL;DR

This paper introduces approximate vertex enumeration algorithms for polytopes, allowing for near solutions that are computationally easier, especially in lower dimensions, with proven correctness in specific cases.

Contribution

The paper presents two novel algorithms for approximate vertex enumeration, including an approximate variant of Motzkin's method and a dimension-restricted shortcut algorithm, with proven correctness in certain conditions.

Findings

Algorithms are correct in 2D and 3D cases.

Approximate methods are often easier to compute than exact enumeration.

Numerical examples demonstrate practical advantages of approximation.

Abstract

The problem to compute the vertices of a polytope given by affine inequalities is called vertex enumeration. The inverse problem, which is equivalent by polarity, is called the convex hull problem. We introduce `approximate vertex enumeration' as the problem to compute the vertices of a polytope which is close to the original polytope given by affine inequalities. In contrast to exact vertex enumerations, both polytopes are not required to be combinatorially equivalent. Two algorithms for this problem are introduced. The first one is an approximate variant of Motzkin's double description method. Only under certain strong conditions, which are not acceptable for practical reasons, we were able to prove correctness of this method for polytopes of arbitrary dimension. The second method, called shortcut algorithm, is based on constructing a plane graph and is restricted to polytopes of dimension 2 and 3. We prove correctness of the shortcut algorithm. As a consequence, we also obtain correctness of the approximate double description method, only for dimension 2 and 3 but without any restricting conditions as still required for higher dimensions. We show that for dimension 2 and 3 both algorithm remain correct if imprecise arithmetic is used and the computational error caused by imprecision is not too high. Both algorithms were implemented. The numerical examples motivate the approximate vertex enumeration problem by showing that the approximate problem is often easier to solve than the exact vertex enumeration problem. It remains open whether or not the approximate double description method (without any restricting condition) is correct for polytopes of dimension 4 and higher.