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

Whenever someone takes Implement `einsum` equivalent 路 Issue #57 路 pymc-devs/pytensor 路 GitHub :smiley: