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

Performs EM update for iterative calculation of Bayesian multiscale MAP estimate in Nowak/Kolaczyk algorithm. Called by BMSPoisInv.m -Performs EM update for iterative calculation of Bayesian multiscale MAP estimate in Nowak/Kolaczy