The tamed unadjusted Langevin algorithm
In this article, we consider the problem of sampling from a probability measure π having a density on Rd proportional to x↦e−U(x). The Euler discretization of the Langevin stochastic differential equation (SDE) is known to be unstable, when the potential U is superlinear. Based on previous works on...
Saved in:
Published in | Stochastic processes and their applications Vol. 129; no. 10; pp. 3638 - 3663 |
---|---|
Main Authors | , , , |
Format | Journal Article |
Language | English |
Published |
Elsevier B.V
01.10.2019
Elsevier |
Subjects | |
Online Access | Get full text |
Cover
Loading…
Summary: | In this article, we consider the problem of sampling from a probability measure π having a density on Rd proportional to x↦e−U(x). The Euler discretization of the Langevin stochastic differential equation (SDE) is known to be unstable, when the potential U is superlinear. Based on previous works on the taming of superlinear drift coefficients for SDEs, we introduce the Tamed Unadjusted Langevin Algorithm (TULA) and obtain non-asymptotic bounds in V-total variation norm and Wasserstein distance of order 2 between the iterates of TULA and π, as well as weak error bounds. Numerical experiments are presented which support our findings. |
---|---|
ISSN: | 0304-4149 1879-209X |
DOI: | 10.1016/j.spa.2018.10.002 |