Speedup implementation of multivariate_normal and allow `method` of covariance decomposition#1203

ricardoV94 merged 2 commits intopymc-devs:mainpymc-devs/pytensor:mainfrom

ricardoV94:speedup_mvnormalricardoV94/pytensor:speedup_mvnormalCopy head branch name to clipboard