I've just noticed the existance of the library numba-scipy, which adds more an API for numba to use scipy functions, such as "special" functions. This would be great to speed up some functions, such as "phase_covariance", which can then be much faster in making covariance matrices, etc...