Computes multivariate normal (MVN) densities, and samples from MVN distributions, when the covariance or precision matrix is sparse.