Volesti: Volume Approximation and Sampling for Convex Polytopes in R
arXiv:2007.01578 · doi:10.32614/RJ-2021-077
Abstract
Sampling from high dimensional distributions and volume approximation of convex bodies are fundamental operations that appear in optimization, finance, engineering, artificial intelligence and machine learning. In this paper we present volesti, an R package that provides efficient, scalable algorithms for volume estimation, uniform and Gaussian sampling from convex polytopes. volesti scales to hundreds of dimensions, handles efficiently three different types of polyhedra and provides non existing sampling routines to R. We demonstrate the power of volesti by solving several challenging problems using the R language.
arXiv comment19 pages, 8 figures, 3 tables
References in corpus (5)
- A Conceptual Introduction to Hamiltonian Monte Carlo
- Efficient Random-Walk Methods for Approximating Polytope Volume
- Fast MCMC sampling algorithms on polytopes
- Theoretical guarantees for approximate sampling from smooth and log-concave densities
- Practical volume computation of structured convex bodies, and an application to modeling portfolio dependencies and financial crises