Function helixiser::bessel_utils::J1b

pub fn J1b(x: f64) -> f64

Computes J1(x), for large values of x

for 'large' values of x we use approximations described in: Harrison, John. "Fast and accurate Bessel function computation." 2009 19th IEEE Symposium on Computer Arithmetic. IEEE, 2009.