int d(1); for(int c(0);c<mu;c++) d *= 2 ;
int d(1<<mu);
void dwf_conserved_axial_ps_corr(LatticeComplex& corr,				 const multi1d<LatticeColorMatrix>& u,				 const multi1d<LatticePropagator>& p5d, 				 const int mu){  // gamma mapping G_0 --> Gamma(1)  //               G_1 --> Gamma(2)  //               G_2 --> Gamma(4)  //               G_3 --> Gamma(8)  int d(1);  for(int c(0);c<mu;c++)    d *= 2 ;  int g5(Ns*Ns - 1) ;  LatticeComplex  C(0.0) ;  corr = 0.0 ;  int N5(p5d.size());  multi1d<LatticePropagator> us_p5d(N5) ;  for(int s(0); s<N5;s++)    us_p5d[s] = u[mu]*shift(p5d[s],FORWARD,mu) ;    for(int s(0); s<N5;s++){    // first the 1-gamma_mu term     C=0.5*(  trace(adj(   p5d[N5-1-s])*Gamma(g5)*Gamma(d)*us_p5d[s]) -	     trace(adj(   p5d[N5-1-s])*Gamma(g5)         *us_p5d[s]) +	     //now the 1+gamma_mu term	     trace(adj(us_p5d[N5-1-s])*Gamma(g5)*Gamma(d)*   p5d[s]) +	     trace(adj(us_p5d[N5-1-s])*Gamma(g5)         *   p5d[s])   );    if(s<N5/2)       corr -= C ;    else      corr += C ;  }    }