Projected Langevin Monte Carlo algorithms in non-convex and super-linear setting

Abstract

It is of significant interest in many applications to sample from a high-dimensional target distribution π with the density π(d x) e-U(x) (d x) , based on the temporal discretization of the Langevin stochastic differential equations (SDEs). In this paper, we propose an explicit projected Langevin Monte Carlo (PLMC) algorithm with non-convex potential U and super-linear gradient of U and investigate the non-asymptotic analysis of its sampling error in total variation distance. Equipped with time-independent regularity estimates for the associated Kolmogorov equation, we derive the non-asymptotic bounds on the total variation distance between the target distribution of the Langevin SDEs and the law induced by the PLMC scheme with order O(d\3γ/2 , 2γ-1 \ h | h|), where d is the dimension of the target distribution and γ ≥ 1 characterizes the growth of the gradient of U. In addition, if the gradient of U is globally Lipschitz continuous, an improved convergence order of O(d3/2 h) for the classical Langevin Monte Carlo (LMC) scheme is derived with a refinement of the proof based on Malliavin calculus techniques. To achieve a given precision ε, the smallest number of iterations of the PLMC algorithm is proved to be of order O(d\3γ/2 , 2γ-1 \ε \ · (dε) · (1ε) ). In particular, the classical Langevin Monte Carlo (LMC) scheme with the non-convex potential U and the globally Lipschitz gradient of U can be guaranteed by order O(d3/2ε · (1ε) ). Numerical experiments are provided to confirm the theoretical findings.

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…