A hamiltonian Monte Carlo method for non-smooth energy sampling

Abstract : Efficient sampling from high-dimensional distribu- tions is a challenging issue that is encountered in many large data recovery problems. In this context, sampling using Hamil- tonian dynamics is one of the recent techniques that have been proposed to exploit the target distribution geometry. Such schemes have clearly been shown to be efficient for multidimensional sam- pling but, rather, are adapted to distributions from the exponential family with smooth energy functions. In this paper, we address the problem of using Hamiltonian dynamics to sample from probabil- ity distributions having non-differentiable energy functions such as those based on the l1 norm. Such distributions are being used intensively in sparse signal and image recovery applications. The technique studied in this paper uses a modified leapfrog transform involving a proximal step. The resulting nonsmooth Hamiltonian Monte Carlo method is tested and validated on a number of exper- iments. Results show its ability to accurately sample according to various multivariate target distributions. The proposed technique is illustrated on synthetic examples and is applied to an image denoising problem.
Complete list of metadatas

Cited literature [43 references]  Display  Hide  Download

Contributor : Open Archive Toulouse Archive Ouverte (oatao) <>
Submitted on : Wednesday, October 5, 2016 - 10:31:20 AM
Last modification on : Thursday, January 23, 2020 - 6:22:13 PM
Long-term archiving on: Friday, January 6, 2017 - 1:10:54 PM


Files produced by the author(s)



Lotfi Chaari, Jean-Yves Tourneret, Caroline Chaux, Hadj Batatia. A hamiltonian Monte Carlo method for non-smooth energy sampling. IEEE Transactions on Signal Processing, Institute of Electrical and Electronics Engineers, 2016, vol. 64 (n° 21), pp. 5585-5594. ⟨10.1109/TSP.2016.2585120⟩. ⟨hal-01376544⟩



Record views


Files downloads