def as_group_action_matrices(F, space, list_of_group_elements, basis): n = len(list_of_group_elements) d = len(space) A = [matrix(F, d, d) for i in range(n)] for i, g in enumerate(list_of_group_elements): for j, omega in enumerate(space): omega1 = omega.group_action(g) v1 = omega1.coordinates(basis = basis) A[i][:, j] = vector(v1) return A