pub fn entanglement_entropy_density_subsystem( density: &[Complex64], local_dimensions: &[usize], retained_sites: &[usize], order: EntropyOrder, ) -> Result<f64>