Variational Bayesian approximation of inverse problems using sparse precision matrices

File Description SizeFormat 
2110.11840v2.pdfFile embargoed until 08 March 20232.68 MBAdobe PDF    Request a copy
Title: Variational Bayesian approximation of inverse problems using sparse precision matrices
Authors: Povala, J
Kazlauskaite, I
Febrianto, E
Cirak, F
Girolami, M
Item Type: Journal Article
Abstract: Inverse problems involving partial differential equations (PDEs) are widely used in science and engineering. Although such problems are generally ill-posed, different regularisation approaches have been developed to ameliorate this problem. Among them is the Bayesian formulation, where a prior probability measure is placed on the quantity of interest. The resulting posterior probability measure is usually analytically intractable. The Markov Chain Monte Carlo (MCMC) method has been the go-to method for sampling from those posterior measures. MCMC is computationally infeasible for large-scale problems that arise in engineering practice. Lately, Variational Bayes (VB) has been recognised as a more computationally tractable method for Bayesian inference, approximating a Bayesian posterior distribution with a simpler trial distribution by solving an optimisation problem. In this work, we argue, through an empirical assessment, that VB methods are a flexible and efficient alternative to MCMC for this class of problems. We propose a natural choice of a family of Gaussian trial distributions parametrised by precision matrices, thus taking advantage of the inherent sparsity of the inverse problem encoded in its finite element discretisation. We utilise stochastic optimisation to efficiently estimate the variational objective and assess not only the error in the solution mean but also the ability to quantify the uncertainty of the estimate. We test this on PDEs based on the Poisson equation in 1D and 2D. A Tensorflow implementation is made publicly available on GitHub.
Issue Date: Apr-2022
Date of Acceptance: 30-Jan-2022
DOI: 10.1016/j.cma.2022.114712
ISSN: 0045-7825
Publisher: Elsevier BV
Start Page: 1
End Page: 31
Journal / Book Title: Computer Methods in Applied Mechanics and Engineering
Volume: 393
Copyright Statement: © 2022 Elsevier Ltd. All rights reserved. This manuscript is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International Licence
Keywords: stat.AP
74S05, 74S60, 62-08, 62P30
Applied Mathematics
01 Mathematical Sciences
09 Engineering
Publication Status: Published
Embargo Date: 2023-03-08
Article Number: 114712
Online Publication Date: 2022-03-09
Appears in Collections:Mathematics

This item is licensed under a Creative Commons License Creative Commons