printf("gaugerec %e max %e factor %e\n", sum1/((double)(VOLUME*g_nproc))/4., max2, g_rgi_C1*g_rgi_C1*g_beta*g_beta/9);
printf("gaugeplaq %e max %e factor %e\n", sum/((double)(VOLUME*g_nproc))/4., max, g_rgi_C0*g_rgi_C0*g_beta*g_beta/9);
void gauge_momenta(double step) {  int i,mu;  static su3 v, w;  su3 *z;  static su3adj deriv;  su3adj *xm;  static double st, st1;#ifdef _COMP_FORCE  double sum=0., sum1=0., max=0., max2=0.;  double sum2=0.;#endif  st  = -step * g_rgi_C0 * g_beta/3.0;  st1 = -step * g_rgi_C1 * g_beta/3.0;  for(i = 0; i < VOLUME; i++){    for(mu=0;mu<4;mu++){      z=&g_gauge_field[i][mu];      xm=&moment[i][mu];      v=get_staples(i,mu);       _su3_times_su3d(w,*z,v);      _trace_lambda(deriv,w);#ifdef _COMP_FORCE      sum2 = _su3adj_square_norm(deriv);      sum+= sum2;      if(sum2 > max) max = sum2;#endif      _minus_const_times_mom(*xm,st,deriv);      if(g_rgi_C1 > 0. || g_rgi_C1 < 0.) {	get_rectangle_staples(&v, i, mu);	_su3_times_su3d(w, *z, v);	_trace_lambda(deriv, w);#ifdef _COMP_FORCE	sum2 =_su3adj_square_norm(deriv);	sum1+= sum2;	if(sum2 > max2) max2 = sum2;#endif	_minus_const_times_mom(*xm, st1, deriv);      }    }  }#ifdef _COMP_FORCE#ifdef MPI  MPI_Reduce(&sum, &sum2, 1, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD);  sum = sum2;  MPI_Reduce(&max, &sum2, 1, MPI_DOUBLE, MPI_MAX, 0, MPI_COMM_WORLD);  max = sum2;  if(g_rgi_C1 > 0. || g_rgi_C1 < 0.) {    MPI_Reduce(&sum1, &sum2, 1, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD);    sum1 = sum2;    MPI_Reduce(&max2, &sum2, 1, MPI_DOUBLE, MPI_MAX, 0, MPI_COMM_WORLD);    max2 = sum2;    if(g_proc_id == 0) {      printf("gaugerec %e max %e factor %e\n", sum1/((double)(VOLUME*g_nproc))/4., max2,	     g_rgi_C1*g_rgi_C1*g_beta*g_beta/9);      fflush(stdout);    }  }#endif  if(g_proc_id == 0) {    printf("gaugeplaq %e max %e factor %e\n", sum/((double)(VOLUME*g_nproc))/4., max, 	   g_rgi_C0*g_rgi_C0*g_beta*g_beta/9);    fflush(stdout);  }#endif}