Multivariate LogNormal distribution

Yes it was a bug, it was already fixed in Allow default_output to be any valid Python index by ricardoV94 路 Pull Request #274 路 pymc-devs/pytensor 路 GitHub

The next release of PyMC should incorporate that :slight_smile: