Dear All, I've updated my gist of wrapping BEAGLE here:
https://gist.github.com/sdwfrost/5c574857bd91648fb7ee The script runs fine, and the transition matrix qmat is being calculated correctly, but my log likelihood is still -Inf (rather than c. -84); is there anything obviously wrong with my Julia code? The C++ code I'm trying to mirror is here: https://github.com/sdwfrost/beagle-lib/blob/master/examples/standalone/hellobeagle/src/hello.cpp and the .h https://github.com/sdwfrost/beagle-lib/blob/master/libhmsbeagle/beagle.h Best Simon
