From 427c8695feaba4de96e865a3754e07f669e6e121 Mon Sep 17 00:00:00 2001 From: Fabian Joswig Date: Fri, 1 Apr 2022 16:20:21 +0100 Subject: [PATCH] Change signs and prefactors for conserved current to mimic the 5d version. --- .../WilsonFermionImplementation.h | 20 +++++++++---------- 1 file changed, 9 insertions(+), 11 deletions(-) diff --git a/Grid/qcd/action/fermion/implementation/WilsonFermionImplementation.h b/Grid/qcd/action/fermion/implementation/WilsonFermionImplementation.h index c82dd172..40e774b3 100644 --- a/Grid/qcd/action/fermion/implementation/WilsonFermionImplementation.h +++ b/Grid/qcd/action/fermion/implementation/WilsonFermionImplementation.h @@ -618,24 +618,22 @@ void WilsonFermion::ContractConservedCurrent(PropagatorField &q_in_1, }; Gamma gmu=Gamma(Gmu[mu]); + g5Lg5=g5*q_in_1*g5; + tmp_shifted=Cshift(q_in_2,mu,1); + Impl::multLinkField(R,this->Umu,tmp_shifted,mu); + gmuR=gmu*R; + + q_out=adj(g5Lg5)*R; + q_out-=adj(g5Lg5)*gmuR; + tmp_shifted=Cshift(q_in_1,mu,1); Impl::multLinkField(g5Lg5,this->Umu,tmp_shifted,mu); g5Lg5=g5*g5Lg5*g5; R=q_in_2; gmuR=gmu*R; - q_out=adj(g5Lg5)*R; - q_out+=adj(g5Lg5)*gmuR; - - g5Lg5=g5*q_in_1*g5; - tmp_shifted=Cshift(q_in_2,mu,1); - Impl::multLinkField(R,this->Umu,tmp_shifted,mu); - gmuR=gmu*R; - q_out-=adj(g5Lg5)*R; - q_out+=adj(g5Lg5)*gmuR; - - q_out*=0.5; + q_out-=adj(g5Lg5)*gmuR; }