This function is very slow. The problem could also be related to the implementation of the `von_mises_cdf` function.