In practise, the likelihood component is optimized using negative sampling with EM for the latent slots.
we use hard EM, to select a single slot before taking gradient steps with respect to the model parameters.