bessel_k works and is fast