pub(crate) fn polarization_from_inputs<const N: usize>( vector: [Jet<N>; 3], pol: u8, ) -> Result<[Jet<N>; 3]>
Chain the local Cartesian polarization derivative into solver parameters.