Function bessel_y_impl 

pub fn bessel_y_impl<P, E, V, const N1: usize, const N2: usize, const N3: usize, const NH: usize, const J1: usize, const J2: usize, const JH: usize, const ORDER_ONE: bool>(
    x: V,
    t: &BesselY<E, N1, N2, N3, NH>,
    tj: &BesselJ<E, J1, J2, JH>,
) -> V
where
    E: FloatElement,
    V: FloatVector<Element = E> + TranscendentalMathWithPolicy,
    P: Policy,

$Y_0(x)$ and $Y_1(x)$, selected by ORDER_ONE.

Both are $\frac{2}{\pi}\ln(x/x_k)J_\nu(x) + \text{(root-factored rational)}$ below 8, and the shared Hankel above it. Undefined for x <= 0: NaN there, and $-\infty$ at 0.