History permutation/gsl_permute_matrix_long.h