https://github.com/quantopian/bayesalpha/blob/master/bayesalpha/dists.py#L52 not sure that line does what is expected `pm.MvNormal.dist(mu=mu, chol=chol).random(point, size)` maybe?