a *= 0.5;
a *= (Real(1)/Real(2));
void taproj(LatticeColorMatrix& a){  START_CODE();  // a = a - a_dagger  --- a -> antihermitian matrix  LatticeColorMatrix aux_1 = a;  a -= adj(aux_1);   if (Nc > 1)  {    // tmp = Tr[ a ] * (1/Nc)    LatticeReal tmp = imag(trace(a))*Real(1/Real(Nc));    // a = a - (1/Nc) * Tr[ a] --- a -> traceless matrix    a -= cmplx(0,tmp);  }  // Normalisation to make taproj idempotent  a *= 0.5;    END_CODE();}