Support marginals with 0 variance in multivariate Gaussian

Numpy supports it: should we? @shakirm, please, what do you think: useful for GP, or rare enough that we don't ned to make it a priority?