DeepAI AI Chat
Log In Sign Up

Multilevel adaptive sparse Leja approximations for Bayesian inverse problems

04/27/2019
by   Ionut-Gabriel Farcas, et al.
Technische Universität München
0

Deterministic interpolation and quadrature methods are often unsuitable to address Bayesian inverse problems depending on computationally expensive forward mathematical models. While interpolation may give precise posterior approximations, deterministic quadrature is usually unable to efficiently investigate an informative and thus concentrated likelihood. This leads to a large number of required expensive evaluations of the mathematical model. To overcome these challenges, we formulate and test a multilevel adaptive sparse Leja algorithm. At each level, adaptive sparse grid interpolation and quadrature are used to approximate the posterior and perform all quadrature operations, respectively. Specifically, our algorithm uses coarse discretizations of the underlying mathematical model to investigate the parameter space and to identify areas of high posterior probability. Adaptive sparse grid algorithms are then used to place points in these areas, and ignore other areas of small posterior probability. The points are weighted Leja points. As the model discretization is coarse, the construction of the sparse grid is computationally efficient. On this sparse grid, the posterior measure can be approximated accurately with few expensive, fine model discretizations. The efficiency of the algorithm can be enhanced further by exploiting more than two discretization levels. We apply the proposed multilevel adaptive sparse Leja algorithm in numerical experiments involving elliptic inverse problems in 2D and 3D space, in which we compare it with Markov chain Monte Carlo sampling and a standard multilevel approximation.

READ FULL TEXT
11/22/2021

Bayesian Inversion of Log-normal Eikonal Equations

We study the Bayesian inverse problem for inferring the log-normal slown...
07/19/2019

Adaptive sampling-based quadrature rules for efficient Bayesian prediction

A novel method is proposed to infer Bayesian predictions of computationa...
05/08/2021

Spline-Based Bayesian Emulators for Large Scale Spatial Inverse Problems

A Bayesian approach to nonlinear inverse problems is considered where th...
04/15/2019

Approximate Bayesian Inference via Sparse grid Quadrature Evaluation for Hierarchical Models

We combine conditioning techniques with sparse grid quadrature rules to ...
04/28/2022

Multilevel Optimization for Inverse Problems

Inverse problems occur in a variety of parameter identification tasks in...
04/05/2021

Multilevel Stein variational gradient descent with applications to Bayesian inverse problems

This work presents a multilevel variant of Stein variational gradient de...