High-Order Langevin Monte Carlo Algorithms
By: Thanh Dang , Mert Gurbuzbalaban , Mohammad Rafiqul Islam and more
Potential Business Impact:
Makes computer learning faster and better.
Langevin algorithms are popular Markov chain Monte Carlo (MCMC) methods for large-scale sampling problems that often arise in data science. We propose Monte Carlo algorithms based on the discretizations of $P$-th order Langevin dynamics for any $P\geq 3$. Our design of $P$-th order Langevin Monte Carlo (LMC) algorithms is by combining splitting and accurate integration methods. We obtain Wasserstein convergence guarantees for sampling from distributions with log-concave and smooth densities. Specifically, the mixing time of the $P$-th order LMC algorithm scales as $O\left(d^{\frac{1}{R}}/\epsilon^{\frac{1}{2R}}\right)$ for $R=4\cdot 1_{\{ P=3\}}+ (2P-1)\cdot 1_{\{ P\geq 4\}}$, which has a better dependence on the dimension $d$ and the accuracy level $\epsilon$ as $P$ grows. Numerical experiments illustrate the efficiency of our proposed algorithms.
Similar Papers
The Picard-Lagrange Framework for Higher-Order Langevin Monte Carlo
Statistics Theory
Makes computer learning faster and more accurate.
Regime-Switching Langevin Monte Carlo Algorithms
Computation
Makes computer learning faster and better.
Underdamped Langevin MCMC with third order convergence
Machine Learning (Stat)
Makes computer learning faster and more accurate.