in zfci and zcas we should prefer gaunt and breit from input, but take relref as backup.
essentially enables coulomb input orbitals for gaunt/breit or vice versa.
ZQvec now with gaunt interaction ; not verified for all cases but doesn't affect coulomb
some comments on qvec
use global constant for speed of light