Modified computeKSLinkQuda to compute unitarized links if required.

Previously, this function only computed the fat and long links.
Unitarization was applied by calling another routine in the milc-interface.
I have modified computeKSLinkQuda to compute the fat and long links and
then apply unitarization to the fatlink if required.