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

[plantbayes]

Thanks!

q is relatively small (< 5), but N is fairly big (several thousand observations).
Do you think that is probably too big? If not, is there an efficient way to convert my (N,q,q) tensor into the block-diagonal format?