Hessian-based adaptive sparse quadrature for infinite-dimensional Bayesian inverse problems
From MaRDI portal
Publication:2310097
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 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.
Recommendations
- Sparse, adaptive Smolyak quadratures for Bayesian inverse problems
- Bayes-Hermite quadrature
- Multilevel adaptive sparse Leja approximations for Bayesian inverse problems
- Adaptive Hessian-Based Nonstationary Gaussian Process Response Surface Method for Probability Density Approximation with Application to Bayesian Solution of Large-Scale Inverse Problems
- Scaling limits in computational Bayesian inversion
Cites work
- A computational framework for infinite-dimensional Bayesian inverse problems. I: The linearized case, with application to global seismic inversion
- A computational framework for infinite-dimensional Bayesian inverse problems. II: stochastic Newton MCMC with application to ice sheet flow inverse problems
- A new algorithm for high-dimensional uncertainty quantification based on dimension-adaptive sparse grid approximation and reduced basis methods
- A stochastic Newton MCMC method for large-scale statistical inverse problems with application to seismic inversion
- An adaptive sparse grid algorithm for elliptic PDEs with lognormal diffusion coefficient
- An implicit restarted Lanczos method for large symmetric eigenvalue problems
- Complexity analysis of accelerated MCMC methods for Bayesian inversion
- Convergence rates of best \(N\)-term Galerkin approximations for a class of elliptic SPDEs
- Dimension-adaptive tensor-product quadrature
- Dimension-independent likelihood-informed MCMC
- Finding structure with randomness: probabilistic algorithms for constructing approximate matrix decompositions
- High-dimensional adaptive sparse polynomial interpolation and applications to parametric PDEs
- scientific article; zbMATH DE number 3622441 (Why is no real title available?)
- scientific article; zbMATH DE number 2243098 (Why is no real title available?)
- Inverse problems: a Bayesian perspective
- MAP estimators and their consistency in Bayesian nonparametric inverse problems
- Optimal low-rank approximations of Bayesian linear inverse problems
- Riemann manifold Langevin and Hamiltonian Monte Carlo methods. With discussion and authors' reply
- Scalable and efficient algorithms for the propagation of uncertainty from data through inference to prediction for large-scale problems, with application to flow of the antarctic ice sheet
- Scaling limits in computational Bayesian inversion
- Sparse deterministic approximation of Bayesian inverse problems
- Sparse, adaptive Smolyak quadratures for Bayesian inverse problems
- Sparse-grid, reduced-basis Bayesian inversion
- Sparse-grid, reduced-basis Bayesian inversion: nonaffine-parametric nonlinear equations
Cited in
(30)- Derivative-informed projected neural networks for high-dimensional parametric maps governed by PDEs
- Variational inference for nonlinear inverse problems via neural net kernels: comparison to Bayesian neural networks, application to topology optimization
- Non-asymptotic error estimates for the Laplace approximation in Bayesian inverse problems
- Optimal design of acoustic metamaterial cloaks under uncertainty
- Variational Bayesian approximation of inverse problems using sparse precision matrices
- On the convergence of the Laplace approximation and noise-level-robustness of Laplace-based Monte Carlo methods for Bayesian inverse problems
- Taylor approximation and variance reduction for PDE-constrained optimal control under uncertainty
- Bayesian inference of heterogeneous epidemic models: application to COVID-19 spread accounting for long-term care facilities
- Efficient estimation of hydraulic conductivity heterogeneity with non-redundant measurement information
- Sparse polynomial approximations for affine parametric saddle point problems
- Sparse, adaptive Smolyak quadratures for Bayesian inverse problems
- Adaptive Hessian-Based Nonstationary Gaussian Process Response Surface Method for Probability Density Approximation with Application to Bayesian Solution of Large-Scale Inverse Problems
- Sparse quadrature for high-dimensional integration with Gaussian measure
- Quasi-Monte Carlo Bayesian estimation under Besov priors in elliptic inverse problems
- Stein variational reduced basis Bayesian inversion
- Multilevel Quasi-Monte Carlo Uncertainty Quantification for Advection-Diffusion-Reaction
- Tensor train construction from tensor actions, with application to compression of large high order derivative tensors
- Taylor approximation for chance constrained optimization problems governed by partial differential equations with high-dimensional random parameters
- Multilevel adaptive sparse Leja approximations for Bayesian inverse problems
- Generalized Modes in Bayesian Inverse Problems
- Projected Wasserstein Gradient Descent for High-Dimensional Bayesian Inference
- An Offline-Online Decomposition Method for Efficient Linear Bayesian Goal-Oriented Optimal Experimental Design: Application to Optimal Sensor Placement
- Learning physics-based models from data: perspectives from inverse problems and model reduction
- Finite element discretizations for variable-order fractional diffusion problems
- A Fast and Scalable Computational Framework for Large-Scale High-Dimensional Bayesian Optimal Experimental Design
- Large-scale Bayesian optimal experimental design with derivative-informed projected neural network
- Derivative-informed neural operator: an efficient framework for high-dimensional parametric derivative learning
- How sharp are error bounds? -- Lower bounds on quadrature worst-case errors for analytic functions --
- Accurate, scalable, and efficient Bayesian optimal experimental design with derivative-informed neural operators
- Hessian initialization strategies for \(\ell \)-BFGS solving non-linear inverse problems
This page was built for publication: Hessian-based adaptive sparse quadrature for infinite-dimensional Bayesian inverse problems
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2310097)
