portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-13 01:05:36 +00:00
Code Releases Activity

New loops. Revisit as accelerator loop in future audit

Browse Source
This commit is contained in:
Peter Boyle 2019-06-15 12:44:29 +01:00
parent 16e5d7945e
commit 5254ede2d8
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
Grid/qcd/action/fermion/GparityWilsonImpl.h
View File

@ -225,7 +225,7 @@ public:
auto Uds_v = Uds.View(); auto Uds_v = Uds.View();
auto Uconj_v = Uconj.View(); auto Uconj_v = Uconj.View();
auto Utmp_v= Utmp.View(); auto Utmp_v= Utmp.View();
thread_loop( (auto ss=U_v.begin();ss<U_v.end();ss++),{ thread_foreach(ss,U_v,{
Uds_v[ss](0)(mu) = U_v[ss](); Uds_v[ss](0)(mu) = U_v[ss]();
Uds_v[ss](1)(mu) = Uconj_v[ss](); Uds_v[ss](1)(mu) = Uconj_v[ss]();
}); });
@ -238,7 +238,7 @@ public:
Utmp = where(coor==0,Uconj,Utmp); Utmp = where(coor==0,Uconj,Utmp);
} }
thread_loop((auto ss=Utmp_v.begin();ss<Utmp_v.end();ss++),{ thread_foreach(ss,Utmp_v,{
Uds_v[ss](0)(mu+4) = Utmp_v[ss](); Uds_v[ss](0)(mu+4) = Utmp_v[ss]();
}); });
@ -247,7 +247,7 @@ public:
Utmp = where(coor==0,U,Utmp); Utmp = where(coor==0,U,Utmp);
} }
thread_loop((auto ss=Utmp_v.begin();ss<Utmp_v.end();ss++),{ thread_foreach(ss,Utmp_v,{
Uds_v[ss](1)(mu+4) = Utmp_v[ss](); Uds_v[ss](1)(mu+4) = Utmp_v[ss]();
}); });
@ -262,7 +262,7 @@ public:
auto tmp = TraceIndex<SpinIndex>(outerProduct(Btilde, A)); auto tmp = TraceIndex<SpinIndex>(outerProduct(Btilde, A));
auto link_v = link.View(); auto link_v = link.View();
auto tmp_v = tmp.View(); auto tmp_v = tmp.View();
thread_loop((auto ss = tmp_v.begin(); ss < tmp_v.end(); ss++), { thread_foreach(ss,tmp_v,{
link_v[ss]() = tmp_v[ss](0, 0) + conjugate(tmp_v[ss](1, 1)); link_v[ss]() = tmp_v[ss](0, 0) + conjugate(tmp_v[ss](1, 1));
}); });
PokeIndex<LorentzIndex>(mat, link, mu); PokeIndex<LorentzIndex>(mat, link, mu);
@ -297,7 +297,7 @@ public:
auto tmp_v = tmp.View(); auto tmp_v = tmp.View();
auto Atilde_v = Atilde.View(); auto Atilde_v = Atilde.View();
auto Btilde_v = Btilde.View(); auto Btilde_v = Btilde.View();
thread_loop((int ss = 0; ss < tmp.Grid()->oSites(); ss++) ,{ thread_for(ss,tmp.Grid()->oSites(),{
for (int s = 0; s < Ls; s++) { for (int s = 0; s < Ls; s++) {
int sF = s + Ls * ss; int sF = s + Ls * ss;
auto ttmp = traceIndex<SpinIndex>(outerProduct(Btilde_v[sF], Atilde_v[sF])); auto ttmp = traceIndex<SpinIndex>(outerProduct(Btilde_v[sF], Atilde_v[sF]));