print_point_group_matrix_rep Subroutine

    subroutine print_point_group_matrix_rep(states)
        ! output the symmetry operation of a specific point group in a given
        ! basis
        integer, intent(in) :: states(:, :)
        character(*), parameter :: this_routine = "print_point_group_matrix_rep"

        if (lat%get_ndim() == 2) then
            ! for 2D it is the 8 fold point group symmetry and S^2 EV..
            ! activate all the symmetries
            t_symmetry_mirror = .true.
            t_symmetry_rotation = .true.
            t_symmetry_inversion = .true.

            call print_d4h_pg(states)
        else
            call stop_all(this_routine, "not yet implemented!")
        end if

    end subroutine print_point_group_matrix_rep