pub fn entanglement_spectrum_subsystem( state: &[Complex64], local_dimensions: &[usize], retained_sites: &[usize], ) -> Result<Vec<f64>>