portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-09 21:50:45 +01:00
Code Releases Activity

Change signs and prefactors for conserved current to mimic the 5d

Browse Source 
version.
This commit is contained in:
Fabian Joswig 2022-04-01 16:20:21 +01:00
parent 603fd96747
commit 427c8695fe
1 changed files with 9 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
Grid/qcd/action/fermion/implementation/WilsonFermionImplementation.h
View File

@ -618,24 +618,22 @@ void WilsonFermion<Impl>::ContractConservedCurrent(PropagatorField &q_in_1,
}; };
Gamma gmu=Gamma(Gmu[mu]); 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); tmp_shifted=Cshift(q_in_1,mu,1);
Impl::multLinkField(g5Lg5,this->Umu,tmp_shifted,mu); Impl::multLinkField(g5Lg5,this->Umu,tmp_shifted,mu);
g5Lg5=g5*g5Lg5*g5; g5Lg5=g5*g5Lg5*g5;
R=q_in_2; R=q_in_2;
gmuR=gmu*R; 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)*R;
q_out+=adj(g5Lg5)*gmuR; q_out-=adj(g5Lg5)*gmuR;
q_out*=0.5;
} }