Log-concave Sampling from a Convex Body with a Barrier: a Robust and Unified Dikin Walk
Abstract
We consider the problem of sampling from a d-dimensional log-concave distribution π(θ) (-f(θ)) for L-Lipschitz f, constrained to a convex body with an efficiently computable self-concordant barrier function, contained in a ball of radius R with a w-warm start. We propose a robust sampling framework that computes spectral approximations to the Hessian of the barrier functions in each iteration. We prove that for polytopes that are described by n hyperplanes, sampling with the Lee-Sidford barrier function mixes within O((d2+dL2R2)(w/δ)) steps with a per step cost of O(ndω-1), where ω≈ 2.37 is the fast matrix multiplication exponent. Compared to the prior work of Mangoubi and Vishnoi, our approach gives faster mixing time as we are able to design a generalized soft-threshold Dikin walk beyond log-barrier. We further extend our result to show how to sample from a d-dimensional spectrahedron, the constrained set of a semidefinite program, specified by the set \x∈ Rd: Σi=1d xi Ai C \ where A1,…,Ad, C are n× n real symmetric matrices. We design a walk that mixes in O((nd+dL2R2)(w/δ)) steps with a per iteration cost of O(nω+n2d3ω-5). We improve the mixing time bound of prior best Dikin walk due to Narayanan and Rakhlin that mixes in O((n2d3+n2dL2R2)(w/δ)) steps.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.