Higher Order Langevin Monte Carlo Algorithm

Abstract

A new (unadjusted) Langevin Monte Carlo (LMC) algorithm with improved rates in total variation and in Wasserstein distance is presented. All these are obtained in the context of sampling from a target distribution π that has a density π on Rd known up to a normalizing constant. Moreover, - π is assumed to have a locally Lipschitz gradient and its third derivative is locally H\"older continuous with exponent β ∈ (0,1]. Non-asymptotic bounds are obtained for the convergence to stationarity of the new sampling method with convergence rate 1+ β/2 in Wasserstein distance, while it is shown that the rate is 1 in total variation even in the absence of convexity. Finally, in the case where - π is strongly convex and its gradient is Lipschitz continuous, explicit constants are provided.

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…