fn duality_breaking( matrix: &DMatrix<Complex>, pol: &[u8], ) -> (f64, DMatrix<Complex>)
The fraction of scattering that flips helicity and its gradient, for a matrix at unit scale.