[Boost-users] cyl_bessel_j performance

Hi, cyl_bessel_j is several times slower than jn from glibc for me. What do I do wrong? Bessel function from boost.math seems to be pretty good coded and optimized but in a runtime it baffles me.