Multilevel and Multi-index Monte Carlo methods for the McKean-Vlasov equation
We address the approximation of functionals depending on a system of particles, described by stochastic differential equations (SDEs), in the mean-field limit when the number of particles approaches infinity. This problem is equivalent to estimating the weak solution of the limiting McKean-Vlasov SD...
| Main Authors: | , |
|---|---|
| Format: | Journal article |
| Published: |
Springer Verlag
2017
|
| _version_ | 1826262265189367808 |
|---|---|
| author | Haji-Ali, A Tempone, R |
| author_facet | Haji-Ali, A Tempone, R |
| author_sort | Haji-Ali, A |
| collection | OXFORD |
| description | We address the approximation of functionals depending on a system of particles, described by stochastic differential equations (SDEs), in the mean-field limit when the number of particles approaches infinity. This problem is equivalent to estimating the weak solution of the limiting McKean-Vlasov SDE. To that end, our approach uses systems with finite numbers of particles and a time-stepping scheme. In this case, there are two discretization parameters: the number of time steps and the number of particles. Based on these two parameters, we consider different variants of the Monte Carlo and Multilevel Monte Carlo (MLMC) methods and show that, in the best case, the optimal work complexity of MLMC, to estimate the functional in one typical setting with an error tolerance of TOL, is O (TOL-3) when using the partitioning estimator and the Milstein time-stepping scheme. We also consider a method that uses the recent Multi-index Monte Carlo method and show an improved work complexity in the same typical setting of O (TOL-2 log(TOL-1)2). Our numerical experiments are carried out on the so-called Kuramoto model, a system of coupled oscillators. |
| first_indexed | 2024-03-06T19:33:40Z |
| format | Journal article |
| id | oxford-uuid:1e529d93-9ddd-436d-8a5e-e458d0dcba15 |
| institution | University of Oxford |
| last_indexed | 2024-03-06T19:33:40Z |
| publishDate | 2017 |
| publisher | Springer Verlag |
| record_format | dspace |
| spelling | oxford-uuid:1e529d93-9ddd-436d-8a5e-e458d0dcba152022-03-26T11:15:45ZMultilevel and Multi-index Monte Carlo methods for the McKean-Vlasov equationJournal articlehttp://purl.org/coar/resource_type/c_dcae04bcuuid:1e529d93-9ddd-436d-8a5e-e458d0dcba15Symplectic Elements at OxfordSpringer Verlag2017Haji-Ali, ATempone, RWe address the approximation of functionals depending on a system of particles, described by stochastic differential equations (SDEs), in the mean-field limit when the number of particles approaches infinity. This problem is equivalent to estimating the weak solution of the limiting McKean-Vlasov SDE. To that end, our approach uses systems with finite numbers of particles and a time-stepping scheme. In this case, there are two discretization parameters: the number of time steps and the number of particles. Based on these two parameters, we consider different variants of the Monte Carlo and Multilevel Monte Carlo (MLMC) methods and show that, in the best case, the optimal work complexity of MLMC, to estimate the functional in one typical setting with an error tolerance of TOL, is O (TOL-3) when using the partitioning estimator and the Milstein time-stepping scheme. We also consider a method that uses the recent Multi-index Monte Carlo method and show an improved work complexity in the same typical setting of O (TOL-2 log(TOL-1)2). Our numerical experiments are carried out on the so-called Kuramoto model, a system of coupled oscillators. |
| spellingShingle | Haji-Ali, A Tempone, R Multilevel and Multi-index Monte Carlo methods for the McKean-Vlasov equation |
| title | Multilevel and Multi-index Monte Carlo methods for the McKean-Vlasov equation |
| title_full | Multilevel and Multi-index Monte Carlo methods for the McKean-Vlasov equation |
| title_fullStr | Multilevel and Multi-index Monte Carlo methods for the McKean-Vlasov equation |
| title_full_unstemmed | Multilevel and Multi-index Monte Carlo methods for the McKean-Vlasov equation |
| title_short | Multilevel and Multi-index Monte Carlo methods for the McKean-Vlasov equation |
| title_sort | multilevel and multi index monte carlo methods for the mckean vlasov equation |
| work_keys_str_mv | AT hajialia multilevelandmultiindexmontecarlomethodsforthemckeanvlasovequation AT temponer multilevelandmultiindexmontecarlomethodsforthemckeanvlasovequation |