Adaptive stepsize algorithms for Langevin dynamics
Fuente:
arXiv
Salvato in:
| Autori principali: | , , , |
|---|---|
| Natura: | Preprint |
| Pubblicazione: |
2024
|
| Soggetti: | |
| Accesso online: | |
| Tags: |
Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
|
| _version_ | 1866929478416990208 |
|---|---|
| author | Leroy, Alix Leimkuhler, Benedict Latz, Jonas Higham, Desmond J. |
| author_facet | Leroy, Alix Leimkuhler, Benedict Latz, Jonas Higham, Desmond J. |
| contents | We discuss the design of an invariant measure-preserving transformed dynamics for the numerical treatment of Langevin dynamics based on rescaling of time, with the goal of sampling from an invariant measure. Given an appropriate monitor function which characterizes the numerical difficulty of the problem as a function of the state of the system, this method allows the stepsizes to be reduced only when necessary, facilitating efficient recovery of long-time behavior. We study both the overdamped and underdamped Langevin dynamics. We investigate how an appropriate correction term that ensures preservation of the invariant measure should be incorporated into a numerical splitting scheme. Finally, we demonstrate the use of the technique in several model systems, including a Bayesian sampling problem with a steep prior. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2403_11993 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | Adaptive stepsize algorithms for Langevin dynamics Leroy, Alix Leimkuhler, Benedict Latz, Jonas Higham, Desmond J. Numerical Analysis 65C30, 65C05, 65P10 We discuss the design of an invariant measure-preserving transformed dynamics for the numerical treatment of Langevin dynamics based on rescaling of time, with the goal of sampling from an invariant measure. Given an appropriate monitor function which characterizes the numerical difficulty of the problem as a function of the state of the system, this method allows the stepsizes to be reduced only when necessary, facilitating efficient recovery of long-time behavior. We study both the overdamped and underdamped Langevin dynamics. We investigate how an appropriate correction term that ensures preservation of the invariant measure should be incorporated into a numerical splitting scheme. Finally, we demonstrate the use of the technique in several model systems, including a Bayesian sampling problem with a steep prior. |
| title | Adaptive stepsize algorithms for Langevin dynamics |
| topic | Numerical Analysis 65C30, 65C05, 65P10 |
| url | https://arxiv.org/abs/2403.11993 |