Fix initialization for Dirichlet factor in Gibbs solver