pub fn apply_indirection( mat: &DMatrix<u32>, indices_orig2destin: &[i32], removal_policy: &IndirRemovalPolicy, ) -> (DMatrix<u32>, Vec<bool>)