Multilevel Markov Chain Monte Carlo
In this paper we address the problem of the prohibitively large computational cost of existing Markov chain Monte Carlo methods for large-scale applications with high-dimensional parameter spaces, e.g., in uncertainty quantification in porous media flow. We propose a new multilevel Metropolis--Hasti...
Saved in:
| Main Authors: | , , , |
|---|---|
| Format: | Article (Journal) |
| Language: | English |
| Published: |
August 7, 2019
|
| In: |
SIAM review
Year: 2019, Volume: 61, Issue: 3, Pages: 509-545 |
| ISSN: | 1095-7200 |
| DOI: | 10.1137/19M126966X |
| Online Access: | Verlag, Volltext: https://doi.org/10.1137/19M126966X Verlag, Volltext: https://epubs.siam.org/doi/10.1137/19M126966X 
|
| Author Notes: | T.J. Dodwell, C. Ketelsen, R. Scheichl, A.L. Teckentrup |
| Summary: | In this paper we address the problem of the prohibitively large computational cost of existing Markov chain Monte Carlo methods for large-scale applications with high-dimensional parameter spaces, e.g., in uncertainty quantification in porous media flow. We propose a new multilevel Metropolis--Hastings algorithm and give an abstract, problem-dependent theorem on the cost of the new multilevel estimator based on a set of simple, verifiable assumptions. For a typical model problem in subsurface flow, we then provide a detailed analysis of these assumptions and show significant gains over the standard Metropolis--Hastings estimator. Numerical experiments confirm the analysis and demonstrate the effectiveness of the method with consistent reductions of more than an order of magnitude in the cost of the multilevel estimator over the standard Metropolis--Hastings algorithm for tolerances $\varepsilon < 10^{-2}$. |
|---|---|
| Item Description: | Published electronically August 7, 2019 Gesehen am 05.02.2020 |
| Physical Description: | Online Resource |
| ISSN: | 1095-7200 |
| DOI: | 10.1137/19M126966X |