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

How big are we talking for N and q? A bunch of independent MvN distributions is equivalent to a single MvN, if you stack up the means and concatenate the covariance matrices block-diagonally. It might get slow if N * q is too huge, though.