reunit(u_tmp[mu], lbad, numbad, REUNITARIZE_ERROR);
reunit(u_tmp, lbad, numbad, REUNITARIZE_ERROR);
void unitarityCheck(const multi1d<LatticeColorMatrix>& u){  LatticeBoolean lbad;  int numbad;  for (int mu=0; mu < Nd; ++mu)  {    LatticeColorMatrix u_tmp = u[mu];    reunit(u_tmp[mu], lbad, numbad, REUNITARIZE_ERROR);  }}