Efficient algorithms for Bayesian Inverse Problems with Whittle--Mat\'ern Priors
From MaRDI portal
PDEs with randomness, stochastic partial differential equations (35R60) Bayesian problems; characterization of Bayes procedures (62C10) Bayesian inference (62F15) Iterative numerical methods for linear systems (65F10) Orthogonalization in numerical linear algebra (65F25) Computational methods for sparse matrices (65F50) Numerical methods for inverse problems for boundary value problems involving PDEs (65N21) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30)
Abstract: This paper tackles efficient methods for Bayesian inverse problems with priors based on Whittle--Mat'ern Gaussian random fields. The Whittle--Mat'ern prior is characterized by a mean function and a covariance operator that is taken as a negative power of an elliptic differential operator. This approach is flexible in that it can incorporate a wide range of prior information including non-stationary effects, but it is currently computationally advantageous only for integer values of the exponent. In this paper, we derive an efficient method for handling all admissible noninteger values of the exponent. The method first discretizes the covariance operator using finite elements and quadrature, and uses preconditioned Krylov subspace solvers for shifted linear systems to efficiently apply the resulting covariance matrix to a vector. This approach can be used for generating samples from the distribution in two different ways: by solving a stochastic partial differential equation, and by using a truncated Karhunen-Lo`eve expansion. We show how to incorporate this prior representation into the infinite-dimensional Bayesian formulation, and show how to efficiently compute the maximum a posteriori estimate, and approximate the posterior variance. Although the focus of this paper is on Bayesian inverse problems, the techniques developed here are applicable to solving systems with fractional Laplacians and Gaussian random fields. Numerical experiments demonstrate the performance and scalability of the solvers and their applicability to model and real-data inverse problems in tomography and a time-dependent heat equation.
Recommendations
- Whittle-Matérn priors for Bayesian statistical inversion with applications in electrical impedance tomography
- Semivariogram methods for modeling Whittle-Matérn priors in Bayesian inverse problems
- Bayesian inverse problems with l₁ priors: a randomize-then-optimize approach
- Hyperpriors for Matérn fields with applications in Bayesian inversion
- Stochastic collocation algorithms using \(l_1\)-minimization for Bayesian solution of inverse problems
Cited in
(4)- Uncertainty quantification in scientific computing: a short introduction
- Simultaneous estimation of piecewise constant coefficients in elliptic pdes via Bayesian level-set methods
- Bayesian inference for PDE-based inverse problems using the optimization of a discrete loss
- Subspace projection regularization for large-scale Bayesian linear inverse problems
This page was built for publication: Efficient algorithms for Bayesian Inverse Problems with Whittle--Mat\'ern Priors
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6398666)
