Fixed initialisation and speed-ups (primarily vectorisations and log-likelihood algebra)