If I have an array of points at which I want to evaluate a bessel
function, I should expect better performance by using
gsl_sf_bessel_sequence_Jnu_e, as opposed to gsl_sf_bessel_Jnu, right?
Also, why does the former include a precision mode flag, while the
latter does not?