Non-asymptotic bounds for modified tamed unadjusted Langevin algorithm in non-convex setting
We consider the problem of sampling from a target distribution $\pi_\beta$ on $\mathbb{R}^d$ with density proportional to $\theta\mapsto e^{-\beta U(\theta)}$ using explicit numerical schemes based on discretising the Langevin stochastic differential equation (SDE). In recent literature, taming has...
Saved in:
Main Author: | |
---|---|
Other Authors: | |
Format: | Final Year Project |
Language: | English |
Published: |
Nanyang Technological University
2022
|
Subjects: | |
Online Access: | https://hdl.handle.net/10356/156899 |
Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
Institution: | Nanyang Technological University |
Language: | English |
id |
sg-ntu-dr.10356-156899 |
---|---|
record_format |
dspace |
spelling |
sg-ntu-dr.10356-1568992023-02-28T23:14:40Z Non-asymptotic bounds for modified tamed unadjusted Langevin algorithm in non-convex setting Ng, Matthew Cheng En Ariel Neufeld School of Physical and Mathematical Sciences Ying Zhang ariel.neufeld@ntu.edu.sg Science::Mathematics::Probability theory Science::Mathematics::Applied mathematics We consider the problem of sampling from a target distribution $\pi_\beta$ on $\mathbb{R}^d$ with density proportional to $\theta\mapsto e^{-\beta U(\theta)}$ using explicit numerical schemes based on discretising the Langevin stochastic differential equation (SDE). In recent literature, taming has been proposed and studied as a method for ensuring stability of Langevin-based numerical schemes in the case of super-linearly growing drift coefficients for the Langevin SDE. In particular, the Tamed Unadjusted Langevin Algorithm (TULA) was proposed in recent literature to sample from target distributions with the gradient of the potential $U$ being super-linear. However, theoretical guarantees in Wasserstein distances for Langevin-based algorithms have traditionally been derived assuming strong convexity of the potential $U$. In this paper, we propose a novel taming factor and derive, under a non-convex setting, non-asymptotic theoretical bounds in Wasserstein-1 and Wasserstein-2 distances between the law of our algorithm, which we name the modified Tamed Unadjusted Langevin Algorithm (mTULA), and the target distribution $\pi_\beta$. We obtain resepctive rates of convergence $\mathcal{O}(\lambda)$ and $\mathcal{O}(\lambda^{1/2})$ in Wasserstein-1 and Wasserstein-2 distances for the discretisation error of mTULA in step size $\lambda$. Numerical simulations which support our theoretical findings are presented. Bachelor of Science in Mathematical Sciences 2022-04-27T06:17:42Z 2022-04-27T06:17:42Z 2022 Final Year Project (FYP) Ng, M. C. E. (2022). Non-asymptotic bounds for modified tamed unadjusted Langevin algorithm in non-convex setting. Final Year Project (FYP), Nanyang Technological University, Singapore. https://hdl.handle.net/10356/156899 https://hdl.handle.net/10356/156899 en application/pdf Nanyang Technological University |
institution |
Nanyang Technological University |
building |
NTU Library |
continent |
Asia |
country |
Singapore Singapore |
content_provider |
NTU Library |
collection |
DR-NTU |
language |
English |
topic |
Science::Mathematics::Probability theory Science::Mathematics::Applied mathematics |
spellingShingle |
Science::Mathematics::Probability theory Science::Mathematics::Applied mathematics Ng, Matthew Cheng En Non-asymptotic bounds for modified tamed unadjusted Langevin algorithm in non-convex setting |
description |
We consider the problem of sampling from a target distribution $\pi_\beta$ on $\mathbb{R}^d$ with density proportional to $\theta\mapsto e^{-\beta U(\theta)}$ using explicit numerical schemes based on discretising the Langevin stochastic differential equation (SDE). In recent literature, taming has been proposed and studied as a method for ensuring stability of Langevin-based numerical schemes in the case of super-linearly growing drift coefficients for the Langevin SDE. In particular, the Tamed Unadjusted Langevin Algorithm (TULA) was proposed in recent literature to sample from target distributions with the gradient of the potential $U$ being super-linear. However, theoretical guarantees in Wasserstein distances for Langevin-based algorithms have traditionally been derived assuming strong convexity of the potential $U$. In this paper, we propose a novel taming factor and derive, under a non-convex setting, non-asymptotic theoretical bounds in Wasserstein-1 and Wasserstein-2 distances between the law of our algorithm, which we name the modified Tamed Unadjusted Langevin Algorithm (mTULA), and the target distribution $\pi_\beta$. We obtain resepctive rates of convergence $\mathcal{O}(\lambda)$ and $\mathcal{O}(\lambda^{1/2})$ in Wasserstein-1 and Wasserstein-2 distances for the discretisation error of mTULA in step size $\lambda$. Numerical simulations which support our theoretical findings are presented. |
author2 |
Ariel Neufeld |
author_facet |
Ariel Neufeld Ng, Matthew Cheng En |
format |
Final Year Project |
author |
Ng, Matthew Cheng En |
author_sort |
Ng, Matthew Cheng En |
title |
Non-asymptotic bounds for modified tamed unadjusted Langevin algorithm in non-convex setting |
title_short |
Non-asymptotic bounds for modified tamed unadjusted Langevin algorithm in non-convex setting |
title_full |
Non-asymptotic bounds for modified tamed unadjusted Langevin algorithm in non-convex setting |
title_fullStr |
Non-asymptotic bounds for modified tamed unadjusted Langevin algorithm in non-convex setting |
title_full_unstemmed |
Non-asymptotic bounds for modified tamed unadjusted Langevin algorithm in non-convex setting |
title_sort |
non-asymptotic bounds for modified tamed unadjusted langevin algorithm in non-convex setting |
publisher |
Nanyang Technological University |
publishDate |
2022 |
url |
https://hdl.handle.net/10356/156899 |
_version_ |
1759855446659694592 |