Constraints for multivariate prior in Bayesian update scheme

Ok, I switched to EMCEE with rejection sampling. Not the fastest solution, but it works