A new procedure for constructing double group symmetry functions is pr
esented. Using this method integrals over Hermitian operators can beco
me real quantities, even though the integrand and the functions themse
lves are complex. This is especially of interest to 4-component relati
vistic methods that use the Dirac-Coulomb Hamiltonian directly in elec
tronic structure calculations.