I’m not sure how easy it is to sample with constraints
I think this is what the first model in this example is doing with a univariate normal? Censored Data Models — PyMC example gallery
I guess you can do the same for multivariate normal, but then you have to sample one multivariate per observed point. No idea how that fares vs numerical integration.