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

[ricardoV94]

It accepts an ndarray with dimensions (..., m, m). Not sure what .tolist() does, but if you wrap it in np.asarray and make sure the core dims are the rightmost ones, it will work.