Unadjusted Hamiltonian MCMC with Stratified Monte Carlo Time Integration

Abstract

A randomized time integrator is suggested for unadjusted Hamiltonian Monte Carlo (uHMC) which involves a very minor modification to the usual Verlet time integrator, and hence, is easy to implement. For target distributions of the form μ(dx) e-U(x) dx where U: Rd R 0 is K-strongly convex but only L-gradient Lipschitz, and initial distributions with finite second moment, coupling proofs reveal that an -accurate approximation of the target distribution in L2-Wasserstein distance W2 can be achieved by the uHMC algorithm with randomized time integration using O((d/K)1/3 (L/K)5/3 -2/3 ( W2(μ, ) / )+) gradient evaluations; whereas for such rough target densities the corresponding complexity of the uHMC algorithm with Verlet time integration is in general O((d/K)1/2 (L/K)2 -1 ( W2(μ, ) / )+ ). Metropolis-adjustable randomized time integrators are also provided.

0

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.

Discussion (0)

Sign in to join the discussion.

Loading comments…