please index `grouping_matrix` with `grouping_matrix[i, j]`

`group_map[grouping[i_sample]]` is == to `group_map[grouping_i]`