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