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

Can the logp of that model be evaluated, or does it trigger: Increase support for batched multivariate distributions · Issue #5383 · pymc-devs/pymc · GitHub ?