Faster I_0 Bessel?

[cscherrer]

For the von Mises distribution, we need to compute $$\log I_0(x)$$.

Eq 14 of this paper looks simple to implement:

We could try this and test it against the current implementation