Add density evaluation functions

Now cover MCMC, skewed Gaussian and mixture of Gaussian densities; some
functions made more efficient (vectorize loops); added helpers