fn normal_component<const N: usize>(k: Jet<N>, q: [Jet<N>; 2]) -> Result<Jet<N>>
The normal wavenumber sqrt(k^2 - q.q) on the branch with a nonnegative imaginary part.
sqrt(k^2 - q.q)