Add tests for p_x
Add tests for p_x

Values of optimized p_x are compared to reference p_x (taken from the
original Gibbs sampler code). Another test should be added to test
whether the result is nan (this happens when dividing by zero).

NLeSC/cptm,NLeSC/cptm