OK, thank you for the clarification. I don't really need the full strength of matrix groups. If the command " orbit<Scalar>(Array<Matrix<Scalar>> G, Vector<Scalar> V) The orbit of a vector V under a group generated by G. " works, then I actually only need the permutation group that acts on...