Faster polytope rounding, sampling, and volume computation via a sublinear "Ball Walk"
Abstract
We study the problem of "isotropically rounding" a polytope K⊂Rn, that is, computing a linear transformation which makes the uniform distribution on the polytope have roughly identity covariance matrix. We assume K is defined by m linear inequalities, with guarantee that rB⊂ K⊂ RB, where B is the unit ball. We introduce a new variant of the ball walk Markov chain and show that, roughly, the expected number of arithmetic operations per-step of this Markov chain is O(m) that is sublinear in the input size mn--the per-step time of all prior Markov chains. Subsequently, we give a rounding algorithm that succeeds with probability 1- in O(mn4.5polylog(1,Rr)) arithmetic operations. This gives a factor of n improvement on the previous bound of O(mn5polylog(1,Rr)) for rounding, which uses the hit-and-run algorithm. Since the rounding preprocessing step is in many cases the bottleneck in improving sampling or volume computation, our results imply these tasks can also be achieved in roughly O(mn4.5polylog(1,Rr)+mn4δ-2) operations for computing the volume of K up to a factor 1+δ and O(mn4.5polylog(1,Rr))) for uniformly sampling on K with TV error . This improves on the previous bounds of O(mn5polylog(1,Rr)+mn4δ-2) for volume computation when roughly m≥ n2.5, and O(mn5polylog(1,Rr)) for sampling when roughly m≥ n1.5. We achieve this improvement by a novel method of computing polytope membership, where one avoids checking inequalities estimated to have a very low probability of being violated.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.