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

New loop construct

Browse Source
This commit is contained in:
Peter Boyle 2019-06-15 08:28:45 +01:00
parent b285138be4
commit 52456b9ec7
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
Grid/qcd/action/gauge/GaugeImplTypes.h
View File

@ -77,7 +77,7 @@ public:
static inline void AddLink(Field &U, LinkField &W, int mu) { // U[mu] += W static inline void AddLink(Field &U, LinkField &W, int mu) { // U[mu] += W
auto U_v = U.View(); auto U_v = U.View();
auto W_v = W.View(); auto W_v = W.View();
thread_loop ( (auto ss = 0; ss < U.Grid()->oSites(); ss++),{ thread_for( ss, U.Grid()->oSites(), {
U_v[ss](mu) = U_v[ss](mu) + W_v[ss](); U_v[ss](mu) = U_v[ss](mu) + W_v[ss]();
}); });
} }
@ -134,9 +134,10 @@ public:
//auto start = std::chrono::high_resolution_clock::now(); //auto start = std::chrono::high_resolution_clock::now();
auto U_v = U.View(); auto U_v = U.View();
auto P_v = P.View(); auto P_v = P.View();
thread_loop( (int ss=0;ss<P.Grid()->oSites();ss++),{ thread_for(ss, P.Grid()->oSites(),{
for (int mu = 0; mu < Nd; mu++) for (int mu = 0; mu < Nd; mu++) {
U_v[ss](mu) = ProjectOnGroup(Exponentiate(P_v[ss](mu), ep, Nexp) * U_v[ss](mu)); U_v[ss](mu) = ProjectOnGroup(Exponentiate(P_v[ss](mu), ep, Nexp) * U_v[ss](mu));
}
}); });
//auto end = std::chrono::high_resolution_clock::now(); //auto end = std::chrono::high_resolution_clock::now();