PURE INTEGER FUNCTION RandExcitSymLabelProd(SymLabel1, SymLabel2)
IMPLICIT NONE
INTEGER, INTENT(IN) :: SymLabel1, SymLabel2
IF (tNoSymGenRandExcits) THEN
RandExcitSymLabelProd = 0
else if (tKPntSym) THEN
!Look up the symmetry in the product table for labels (returning labels, not syms)
RandExcitSymLabelProd = SymTableLabels(SymLabel1, SymLabel2)
ELSE
RandExcitSymLabelProd = IEOR(SymLabel1, SymLabel2)
! write(stdout,*) "***",SymLabel1,SymLabel2,RandExcitSymLabelProd
end if
END FUNCTION RandExcitSymLabelProd