pub(crate) fn transverse<const N: usize>(k: Jet<N>, kz: Jet<N>) -> Jet<N>
Transverse wavenumber sqrt(k² - kz²) on the branch with nonnegative imaginary part.
sqrt(k² - kz²)