need to run gibbs sampling from within compute_posterior; will need to save state and restore after sampling is over