Function zernike_basis_impl 

pub fn zernike_basis_impl<P, E, V, const L: usize, const NORM: u8, const N: usize>(
    x: V,
    y: V,
    out: &mut [V; N],
)
where
    P: Policy,
    E: FloatElement,
    V: FloatVector<Element = E> + CoreMath,

Every Zernike mode through degree L at the Cartesian point (x, y).

out[(n(n+2) + m)/2] receives $Z_n^m$ in the normalization named by NORM, and N must equal (L+1)(L+2)/2 (compile-time checked). See the module docs for the algorithm and the domain note.

The policy parameter is unused: evaluation is pure polynomial arithmetic with one fixed FMA-preferring lowering. It is accepted so the signature matches its siblings.