func (m *MoM) Gibbs() {
for d := 0; d < m.M; d++ {
m.unsetZ(d)
z := util.SampleLCounts(m.conditional(d))
m.setZ(d, z)
}
}
func (m *DPMoM) Gibbs() {
for d := 0; d < m.M; d++ {
m.unsetZ(d)
lcounts, mapping := m.conditional(d)
i := util.SampleLCounts(lcounts)
m.setZ(d, mapping[i])
}
}