Add Gibbs updates and log-likelihood computation