Now the sampler works and one can gibbs sample (at least using T+E)