Replace zeros in joint distributions

Instead of in the constraints that result from the IPF step.
This seems to be more in line with what is described in the paper.