Bug in `probe_with_posterior_samples`
Please look at the output of demos/wiener_filter_via_hamiltonian.py in the branch work_on_demos. There is a huge discrepancy between posterior_mean.png and reconstruction.png.
The formulae in StatCalculator seem to be correct.