A transport-based multifidelity preconditioner for Markov chain Monte Carlo

08/28/2018
by   Benjamin Peherstorfer, et al.
0

Markov chain Monte Carlo (MCMC) sampling of posterior distributions arising in Bayesian inverse problems is challenging when evaluations of the forward model are computationally expensive. Replacing the forward model with a low-cost, low-fidelity model often significantly reduces computational cost; however, employing a low-fidelity model alone means that the stationary distribution of the MCMC chain is the posterior distribution corresponding to the low-fidelity model, rather than the original posterior distribution corresponding to the high-fidelity model. We propose a multifidelity approach that combines, rather than replaces, the high-fidelity model with a low-fidelity model. First, the low-fidelity model is used to construct a transport map that deterministically couples a reference Gaussian distribution with an approximation of the low-fidelity posterior. Then, the high-fidelity posterior distribution is explored using a non-Gaussian proposal distribution derived from the transport map. This multifidelity "preconditioned" MCMC approach seeks efficient sampling via a proposal that is explicitly tailored to the posterior at hand and that is constructed efficiently with the low-fidelity model. By relying on the low-fidelity model only to construct the proposal distribution, our approach guarantees that the stationary distribution of the MCMC chain is the high-fidelity posterior. In our numerical examples, our multifidelity approach achieves significant speedups compared to single-fidelity MCMC sampling methods.

READ FULL TEXT

page 13

page 16

research
12/06/2017

Inverse modeling of hydrologic systems with adaptive multi-fidelity Markov chain Monte Carlo simulations

Markov chain Monte Carlo (MCMC) simulation methods are widely used to as...
research
02/05/2021

Invertible Neural Networks versus MCMC for Posterior Reconstruction in Grazing Incidence X-Ray Fluorescence

Grazing incidence X-ray fluorescence is a non-destructive technique for ...
research
04/20/2018

A Bayesian Framework for Assessing the Strength Distribution of Composite Structures with Random Defects

This paper presents a novel stochastic framework to quantify the knock d...
research
12/06/2017

Inverse modeling of hydrologic systems with adaptive multi-fidelity simulations

Markov chain Monte Carlo (MCMC) simulation methods are widely used to as...
research
10/04/2022

Multi-fidelity Monte Carlo: a pseudo-marginal approach

Markov chain Monte Carlo (MCMC) is an established approach for uncertain...
research
05/11/2023

Object based Bayesian full-waveform inversion for shear elastography

We develop a computational framework to quantify uncertainty in shear el...

Please sign up or login with your details

Forgot password? Click here to reset