Spurious global test variable...
Added some comments about unsuccesful BLAS / Chelesky attempts
Finally a faster blas, use chol(W[k]) and trmm!
In stats(), replace a hcat(broadcast) by three loops
write out d/vg.β[k] + vg.ν[k] * sum(Δ .* Δ, 2)
Added some @inbounds, might help
Added testbayes for profiling the variational bayes methods.