portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-10 07:55:35 +00:00
Code Releases Activity

Merge pull request #259 from grid-test-organisation/feature/5d-improvement-fix

Browse Source 
Missing conjugate in MooeeInvDag
This commit is contained in:
Peter Boyle 2019-12-16 04:20:37 -05:00 committed by GitHub
commit 5d834486c9
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
Grid/qcd/action/fermion/implementation/CayleyFermion5Dcache.h
View File

@ -234,12 +234,12 @@ CayleyFermion5D<Impl>::MooeeInvDag (const FermionField &psi_i, FermionField &chi
res = psi(ss+Ls-1) - conjugate(puee[Ls-2])*tmp - acc; res = psi(ss+Ls-1) - conjugate(puee[Ls-2])*tmp - acc;
// Apply L_m^{-\dagger} D^{-dagger} L^{-dagger} // Apply L_m^{-\dagger} D^{-dagger} L^{-dagger}
res = (1.0/pdee[Ls-1])*res; res = conjugate(1.0/pdee[Ls-1])*res;
coalescedWrite(chi[ss+Ls-1],res); coalescedWrite(chi[ss+Ls-1],res);
spProj5m(acc,res); spProj5m(acc,res);
spProj5p(tmp,res); spProj5p(tmp,res);
for (int s=Ls-2;s>=0;s--){ for (int s=Ls-2;s>=0;s--){
res = (1.0/pdee[s])*chi(ss+s) - conjugate(plee[s])*tmp - conjugate(pleem[s])*acc; res = conjugate(1.0/pdee[s])*chi(ss+s) - conjugate(plee[s])*tmp - conjugate(pleem[s])*acc;
spProj5p(tmp,res); spProj5p(tmp,res);
coalescedWrite(chi[ss+s],res); coalescedWrite(chi[ss+s],res);
} }