if( magic_number != 20103) QDP_error_exit("readMILC: unexpected byte order of file");
bool byterev = false ; if( magic_number != 20103){  byterev = true ;  QDPUtil::byte_swap((void *)&magic_number, sizeof(int), 1 ); } if( magic_number != 20103){ QDP_error_exit("readMILC: unexpected magic number"); }
void readMILC(MILCGauge_t& header, multi1d<LatticeColorMatrix>& u, const string& cfg_file){  START_CODE("readMILC");  MILCGaugeInit(header);  // initialize the header with defaults  BinaryReader cfg_in(cfg_file); // for now, cfg_io_location not used  int magic_number;  read(cfg_in, magic_number);  if( magic_number != 20103)    QDP_error_exit("readMILC: unexpected byte order of file");  // Check lattice size  header.nrow.resize(Nd);  read(cfg_in, header.nrow, Nd);  for(int j = 0; j < Nd; ++j)    if ( header.nrow[j] != Layout::lattSize()[j] )      QDP_error_exit("readMILC: unexpected lattice size: header.nrow[%d]=%d",                     j,header.nrow[j]);  // Time stamp  char date_tmp[65];  cfg_in.readArray(date_tmp, 1, 64);  date_tmp[64] = '\0';  header.date = date_tmp;  // Site order - only support non-sitelist format  int order;  read(cfg_in, order);  if( order != 0)    QDP_error_exit("readMILC: only support non-sitelist format");  /*   * Read away...   */    // MILC format has the directions inside the sites  for(int site=0; site < Layout::vol(); ++site)  {    multi1d<int> coord = crtesn(site, Layout::lattSize()); // The coordinate          // Read in Nd SU(3) matrices.     // NOTE: the su3_matrix layout should be the same as in QDP    for(int mu=0; mu < Nd; ++mu)      read(cfg_in, u[mu], coord);    // read in a single site  }  // Go ahead an read checksums, but will not use for now  unsigned int sum29, sum31;  read(cfg_in, sum29);  read(cfg_in, sum31);  cfg_in.close();  END_CODE("readMILC");}