Gibbs sampler working and tested against ensemble sampler
debugging