Added a PR. Please look at the functions. There is difference between the likelihood numbers from the example and the non-vectorized version.