Rényi-infinity constrained sampling with membership queries

Uniform sampling over a convex body is a fundamental algorithmic problem, yet the convergence in KL or R\ényi divergence of most samplers remains poorly understood. In this work, we propose a constrained proximal sampler, a principled and simple algorithm that possesses elegant convergence guarantees. Leveraging the uniform ergodicity of this sampler, we show that it converges in the R\ényi-infinity divergence () with no query complexity overhead when starting from a warm start. This is the strongest of commonly considered performance metrics, implying rates in convergence as special cases. By applying this sampler within an annealing scheme, we propose an algorithm which can approximately sample -close to the uniform distribution on convex bodies in -divergence with query complexity. This improves on all prior results in -divergences, without resorting to any algorithmic modifications or post-processing of the sample. It also matches the prior best known complexity in total variation distance.
View on arXiv