Improved polytope volume calculations based on Hamiltonian Monte Carlo with boundary reflections and sweet arithmetics - Université Côte d'Azur Accéder directement au contenu
Pré-Publication, Document De Travail Année : 2021

Improved polytope volume calculations based on Hamiltonian Monte Carlo with boundary reflections and sweet arithmetics

Amélioration des calculs de volume de polytope basés sur le Monte Carlo Hamiltonien avec des réflexions sur les bords et des arithmétiques édulcorées

Résumé

Computing the volume of a high dimensional polytope is a fundamental problem in geometry, also connected to the calculation of densities of states in statistical physics, and a central building block of such algorithms is the method used to sample a target probability distribution. This paper studies Hamiltonian Monte Carlo (HMC) with reflections on the boundary of a domain, providing an enhanced alternative to Hit-and-run (HAR) to sample a target distribution restricted to the polytope. We make three contributions. First, we provide a convergence bound, paving the way to more precise mixing time analysis. Second, we present a robust implementation based on multi-precision arithmetic, a mandatory ingredient to guarantee exact predicates and robust constructions. We however allow controlled failures to happen, introducing the Sweeten Exact Geometric Computing (SEGC) paradigm. Third, we use our HMC random walk to perform H-polytope volume calculations, using it as an alternative to HAR within the volume algorithm by Cousins and Vempala. The systematic tests conducted up to dimension $n$ = 100 on the cube, the isotropic and the standard simplex show that HMC significantly outperforms HAR both in terms of accuracy and running time. Additional tests show that calculations may be handled up to dimension $n$ = 500. These tests also establish that multiprecision is mandatory to avoid exits from the polytope.
Fichier principal
Vignette du fichier
convexvol-nocolor.pdf (1.33 Mo) Télécharger le fichier
Origine : Fichiers produits par l'(les) auteur(s)

Dates et versions

hal-03048725 , version 1 (09-12-2020)
hal-03048725 , version 2 (14-12-2020)
hal-03048725 , version 3 (11-05-2021)
hal-03048725 , version 4 (07-10-2021)

Identifiants

  • HAL Id : hal-03048725 , version 3

Citer

Augustin Chevallier, Sylvain Pion, Frédéric Cazals. Improved polytope volume calculations based on Hamiltonian Monte Carlo with boundary reflections and sweet arithmetics. 2021. ⟨hal-03048725v3⟩

Relations

323 Consultations
332 Téléchargements

Partager

Gmail Facebook X LinkedIn More