Skip to Main content Skip to Navigation

Efficient high-dimension gaussian sampling based on matrix splitting : application to bayesian Inversion

Abstract : The thesis deals with the problem of high-dimensional Gaussian sampling.Such a problem arises for example in Bayesian inverse problems in imaging where the number of variables easily reaches an order of 106_109. The complexity of the sampling problem is inherently linked to the structure of the covariance matrix. Different solutions to tackle this problem have already been proposed among which we emphasizethe Hogwild algorithm which runs local Gibbs sampling updates in parallel with periodic global synchronisation.Our algorithm makes use of the connection between a class of iterative samplers and iterative solvers for systems of linear equations. It does not target the required Gaussian distribution, instead it targets an approximate distribution. However, we are able to control how far off the approximate distribution is with respect to the required one by means of asingle tuning parameter.We first compare the proposed sampling algorithm with the Gibbs and Hogwild algorithms on moderately sized problems for different target distributions. Our algorithm manages to out perform the Gibbs and Hogwild algorithms in most of the cases. Let us note that the performances of our algorithm are dependent on the tuning parameter.We then compare the proposed algorithm with the Hogwild algorithm on a large scalereal application, namely image deconvolution-interpolation. The proposed algorithm enables us to obtain good results, whereas the Hogwild algorithm fails to converge. Let us note that for small values of the tuning parameter our algorithm fails to converge as well.Not with standing, a suitably chosen value for the tuning parameter enables our proposed sampler to converge and to deliver good results.
Complete list of metadatas

Cited literature [45 references]  Display  Hide  Download
Contributor : Abes Star :  Contact
Submitted on : Friday, June 1, 2018 - 12:03:07 PM
Last modification on : Friday, May 17, 2019 - 11:37:25 AM
Long-term archiving on: : Sunday, September 2, 2018 - 2:08:03 PM


Version validated by the jury (STAR)


  • HAL Id : tel-01804981, version 1


Andrei-Cristian Bӑrbos. Efficient high-dimension gaussian sampling based on matrix splitting : application to bayesian Inversion. Signal and Image processing. Université de Bordeaux, 2018. English. ⟨NNT : 2018BORD0002⟩. ⟨tel-01804981⟩



Record views


Files downloads