Hessian-based adaptive sparse quadrature for infinite-dimensional Bayesian inverse problems

Abstract

In this work we propose and analyze a Hessian-based adaptive sparse quadrature to compute infinite-dimensional integrals with respect to the posterior distribution in the context of Bayesian inverse problems with Gaussian prior. Due to the concentration of the posterior distribution in the domain of the prior distribution, a prior-based parametrization and sparse quadrature may fail to capture the posterior distribution and lead to erroneous evaluation results. By using a parametrization based on the Hessian of the negative log-posterior, the adaptive sparse quadrature can effectively allocate the quadrature points according to the posterior distribution. A dimension-independent convergence rate of the proposed method is established under certain assumptions on the Gaussian prior and the integrands. Dimension-independent and faster convergence than $O(N^{-1/2})$ is demonstrated for a linear as well as a nonlinear inverse problem whose posterior distribution can be effectively approximated by a Gaussian distribution at the MAP point.

References in corpus (4)

• An introduction to sampling via measure transport
• High-Dimensional Stochastic Design Optimization by Adaptive-Sparse Polynomial Dimensional Decomposition
• Importance Sampling: Intrinsic Dimension and Computational Cost
• Sparse Quadrature for High-Dimensional Integration with Gaussian Measure

Cited by in corpus (14)

• Optimal design of acoustic metamaterial cloaks under uncertainty
• Taylor approximation and variance reduction for PDE-constrained optimal control under uncertainty
• Bayesian Poroelastic Aquifer Characterization from InSAR Surface Deformation Data. Part I: Maximum A Posteriori Estimate
• Bayesian inference of heterogeneous epidemic models: Application to COVID-19 spread accounting for long-term care facilities
• Variational Bayesian Approximation of Inverse Problems using Sparse Precision Matrices
• Derivative-Informed Projected Neural Networks for High-Dimensional Parametric Maps Governed by PDEs
• Generalized modes in Bayesian inverse problems
• A fast and scalable computational framework for large-scale and high-dimensional Bayesian optimal experimental design
• Multilevel adaptive sparse Leja approximations for Bayesian inverse problems
• Non-asymptotic error estimates for the Laplace approximation in Bayesian inverse problems
• Variational Inference for Nonlinear Inverse Problems via Neural Net Kernels: Comparison to Bayesian Neural Networks, Application to Topology Optimization
• Taylor approximation for chance constrained optimization problems governed by partial differential equations with high-dimensional random parameters
• How sharp are error bounds? --lower bounds on quadrature worst-case errors for analytic functions
• Stein variational reduced basis Bayesian inversion