pub(crate) fn polar_trig<const N: usize>(theta: Jet<N>) -> [Jet<N>; 2]Expand description
[cos(theta), sin(theta)] with their derivatives. The sine takes the principal
branch (nonnegative real part) and is exactly zero at the poles, where theta is a
floating-point multiple of PI.