In this paper, we propose a Bayesian methodology for receiver function analysis, a key tool in determining the deep structure of the Earth’s crust.We exploit the assumption of sparsity for receiver functions to develop a Bayesian deconvolution

This the MATLAB code that was used to produce the figures and tables in Section V of F. Forbes and G. Fort, Combining Monte Carlo and mean-field like methods for inference in Hidden Markov Random Fields, Accepted for publication in IEEE Trans. on