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