Implementation of the approach for Uncertainty Quantification presented by Pereyra in [1]. The author used the Shepp-Logan phantom. For the reconstruction, the ADMM SALSA method was applied, this implementation uses the C-SALSA2 from [2].
[1] Pereyra, Marcelo. "Maximum-a-posteriori estimation with bayesian confidence regions." SIAM Journal on Imaging Sciences 10.1 (2017): 285-302.
[2] Afonso, Manya V., José M. Bioucas-Dias, and Mário AT Figueiredo. "An augmented Lagrangian approach to the constrained optimization formulation of imaging inverse problems." IEEE Transactions on Image Processing 20.3 (2011): 681-695.