LatticeSpinMatrix  S = traceColor(q3 * quarkContract13(q1*Cgm, Cgm*q2)); #endif
LatticeSpinMatrix deltaContract132(const T1& q1,				   const T2& q2, 				   const T3& q3,				   int n, int m){#if 0  /* C gamma_5 = Gamma(5) = - (C gamma_5)^T */  /* C gamma_mu = Gamma(10) * Gamma(1 << mu) = + (C gamma_mu)^T */  LatticeSpinMatrix  S =     traceColor(q3 * quarkContract13((q1*Gamma(10))*Gamma(n), Gamma(10)*(Gamma(m)*q2)));#else  SpinMatrix g_one = 1.0;  SpinMatrix Cgm = 0.5 * (Gamma(10) * (Gamma(2) * g_one  +  timesI(Gamma(1) * g_one)));  LatticeSpinMatrix  S =     traceColor(q3 * quarkContract13(q1*Cgm, Cgm*q2));#endif  return S;}