cerr << "After initialize vector of latticegauge\n";
multi1d<LatticeGauge> u(Nd);
int main(int argc, char **argv){  // Setup the geometry  const int foo[Nd] = {LX0,LX1};  multi1d<int> nrow(Nd);  nrow = foo;  geom.Init(nrow);  //! Example of calling a plaquette routine  /*! NOTE: the STL is *not* used to hold gauge fields */  cerr << "After initialize vector of latticegauge\n";  Double w_plaq, s_plaq, t_plaq, link;  cerr << "Start gaussian\n";  for(int m=0; m < u.size(); ++m)    gaussian(u[m]);  MesPlq(u, w_plaq, s_plaq, t_plaq, link);  cerr << "w_plaq = " << w_plaq << endl;  cerr << "link = " << link << endl;}