Can confirm - the random for multidimensional mixture is wrong 
I opened an issue: https://github.com/pymc-devs/pymc3/issues/3270 but I am not sure if it would really be fixed - as I have modified this code multiple times and each time I thought that I fixed everything…