pure function TotSymRep()
! Our definition of the totally symmetric representation
! changes according to whether we're using Abelian/k-point
! symmetry or the standard symmetry. It's just a matter of
! convenience, rather than some deep theoretical insight!
Type(Symmetry) TotSymRep
if (TAbelian .or. tUEG) then
TotSymRep%s = 0
else
TotSymRep%s = 1
end if
end function TotSymRep