Function bessel_array 

pub fn bessel_array(
    orders: Vec<f64>,
    arguments: Vec<Complex>,
    kind: Bessel,
    spherical: bool,
    derivative: u8,
) -> Result<(Vec<Complex>, BesselResidual)>

Evaluate elementwise Bessel functions, broadcasting either scalar input.

kind and spherical select the function as in bessel. derivative is 0 or 1, one below the cap of bessel, because the pullback evaluates derivative + 1.

Upstream: the functions of bessel over arrays.