Added function to compute KLD between two multivariate Gaussians