Function partial_trace_density_subsystem 

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

Mixed-state partial trace for an arbitrary retained site set.