PyMC/PyTensor Error with Multivariate Normal Distribution: Incompatible Elemwise Input Shape

Try to call .eval on the cov and mu variables to see if they have the shape you expected them to have