Is it possible to declare multiple MvNormal variables with different covariance matrices in a vectorized way?

It samples without error (the kron one)

edit:

I think I misunderstood @ricardoV94’s question; the model with a single MvN and batch dimension doesn’t sample and runs into the linked issue