powered by
sum of multivariate normal log densities with mean 0 and covariance Sigma, unnormalized
sldmvnorm(z_t, Sigma)