const int mu, const int N5)
const int mu)
void dwf_conserved_axial_ps_corr(LatticeComplex& corr,				 const multi1d<LatticeColorMatrix>& u,				 const multi1d<LatticePropagator>& p5d, 				 const int mu,				 const int N5){  // 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) ;   for(int s(0); s<N5;s++){    // first the 1-gamma_mu term     C=.5*(trace(shift(p5d[s],FORWARD,mu)*adj(p5d[s])*Gamma(g5)*Gamma(d)*u[mu])-	  trace(shift(p5d[s],FORWARD,mu)*adj(p5d[s])*Gamma(g5)         *u[mu])+	  //now th 1+gamma_mu term	  trace(p5d[s]*shift(adj(p5d[s]),FORWARD,mu)*Gamma(g5)*Gamma(d)*adj(u[mu]))+	  trace(p5d[s]*shift(adj(p5d[s]),FORWARD,mu)*Gamma(g5)         *adj(u[mu]))	  );	            if(s<N5/2)       corr -= C ;    else      corr += C ;  }  }