pub fn apply_fermionic_subsystem_phases_density( density: &mut [Complex64], local_dimensions: &[usize], retained_sites: &[usize], ) -> Result<()>
Apply the same fermionic mode permutation to both density-matrix axes.