pub fn entanglement_entropy_density( density: &[Complex64], subsystem_dimension: usize, environment_dimension: usize, order: EntropyOrder, ) -> Result<f64>