portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-25 13:15:55 +01:00
Code Releases Activity

First version of conserved current contraction for Wilson type quarks

Browse Source
This commit is contained in:
Fabian Joswig 2022-03-31 17:02:09 +01:00
parent 605cf401e1
commit 0542eaf1da
1 changed files with 31 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -603,7 +603,37 @@ void WilsonFermion<Impl>::ContractConservedCurrent(PropagatorField &q_in_1,
conformable(_grid, q_in_1.Grid());
conformable(_grid, q_in_2.Grid());
conformable(_grid, q_out.Grid());
assert(0);
PropagatorField tmp_shifted(UGrid);
PropagatorField g5Lg5(UGrid);
PropagatorField R(UGrid);
PropagatorField gmuR(UGrid);
Gamma::Algebra Gmu [] = {
Gamma::Algebra::GammaX,
Gamma::Algebra::GammaY,
Gamma::Algebra::GammaZ,
Gamma::Algebra::GammaT,
};
Gamma gmu=Gamma(Gmu[mu]);
tmp_shifted=Cshift(q_in_1,mu,1);
g5Lg5=g5*tmp_shifted*g5;
R=q_in_2;
gmuR=gmu*R;
qout=adj(g5Lg5)*R;
qout+=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;
qout-=adj(g5Lg5)*R;
qout+=adj(g5Lg5)*gmuR;
qout/=2;
}