Fast methods for multivariate normal distributions