for(mu = 0; mu < Nd; ++mu)
for(int mu = 0; mu < Nd; ++mu)
void rgauge(multi1d<LatticeColorMatrix>& u, LatticeColorMatrix& g){  START_CODE("rgauge");    // g = exp(A), where A = random traceless antihermitian matrix.    // NOTE!!!!: we are using a hack. The twelth_order exponentiation followed  //   by a reunitarization should be replace by an exact exponentiation.  //   So, the gauge field distriubution is NOT from Haar measure, but  //   if someone can think of a really good reason why it should be Haar  //   measure, then something can be done about it.    gaussian(g);  taproj(g);  expm12(g);  reunit(g);      for(mu = 0; mu < Nd; ++mu)  {    LatticeColorMatrix u_tmp = g * u[mu];    u[mu] = u_tmp * shift(adj(g), FORWARD, mu);  }      END_CODE("rgauge");}