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

Merge branch 'develop' into feature/distil

Browse Source 
* develop:
  Missing conjugate in MooeeInvDag
  Allow subspace setup to no converge
  fp16 mandatory. Use SFW is not available as hdw
This commit is contained in:
Michael Marshall 2019-12-20 13:24:59 +00:00
commit df2b0c4e79
2 changed files with 3 additions and 4 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);
} }

3
tests/Test_simd.cc
View File

@ -784,7 +784,6 @@ int main (int argc, char ** argv)
std::cout <<" OK ! "<<std::endl; std::cout <<" OK ! "<<std::endl;
// Double to Half // Double to Half
#ifdef USE_FP16
std::cout << GridLogMessage<< "Double to half" ; std::cout << GridLogMessage<< "Double to half" ;
precisionChange(&H[0],&D[0],Ndp); precisionChange(&H[0],&D[0],Ndp);
precisionChange(&DD[0],&H[0],Ndp); precisionChange(&DD[0],&H[0],Ndp);
@ -813,7 +812,7 @@ int main (int argc, char ** argv)
assert( tmp < 1.0e-3 ); assert( tmp < 1.0e-3 );
} }
std::cout <<" OK ! "<<std::endl; std::cout <<" OK ! "<<std::endl;
#endif
} }
Grid_finalize(); Grid_finalize();
} }