Structured Logconcave Sampling with a Restricted Gaussian Oracle

Abstract

We give algorithms for sampling several structured logconcave families to high accuracy. We further develop a reduction framework, inspired by proximal point methods in convex optimization, which bootstraps samplers for regularized densities to improve dependences on problem conditioning. A key ingredient in our framework is the notion of a "restricted Gaussian oracle" (RGO) for g: Rd → R, which is a sampler for distributions whose negative log-likelihood sums a quadratic and g. By combining our reduction framework with our new samplers, we obtain the following bounds for sampling structured distributions to total variation distance ε. For composite densities (-f(x) - g(x)), where f has condition number and convex (but possibly non-smooth) g admits an RGO, we obtain a mixing time of O( d 3 dε), matching the state-of-the-art non-composite bound; no composite samplers with better mixing than general-purpose logconcave samplers were previously known. For logconcave finite sums (-F(x)), where F(x) = 1nΣi ∈ [n] fi(x) has condition number , we give a sampler querying O(n + (d, nd)) gradient oracles to \fi\i ∈ [n]; no high-accuracy samplers with nontrivial gradient query complexity were previously known. For densities with condition number , we give an algorithm obtaining mixing time O( d 2 dε), improving the prior state-of-the-art by a logarithmic factor with a significantly simpler analysis; we also show a zeroth-order algorithm attains the same query complexity.

0

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.

Discussion (0)

Sign in to join the discussion.

Loading comments…