From 52456b9ec7e060b46d3f17a021b9e9257a631782 Mon Sep 17 00:00:00 2001 From: Peter Boyle Date: Sat, 15 Jun 2019 08:28:45 +0100 Subject: [PATCH] New loop construct --- Grid/qcd/action/gauge/GaugeImplTypes.h | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/Grid/qcd/action/gauge/GaugeImplTypes.h b/Grid/qcd/action/gauge/GaugeImplTypes.h index 3635c8c1..dca0c34f 100644 --- a/Grid/qcd/action/gauge/GaugeImplTypes.h +++ b/Grid/qcd/action/gauge/GaugeImplTypes.h @@ -77,7 +77,7 @@ public: static inline void AddLink(Field &U, LinkField &W, int mu) { // U[mu] += W auto U_v = U.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](); }); } @@ -134,9 +134,10 @@ public: //auto start = std::chrono::high_resolution_clock::now(); auto U_v = U.View(); auto P_v = P.View(); - thread_loop( (int ss=0;ssoSites();ss++),{ - for (int mu = 0; mu < Nd; mu++) + thread_for(ss, P.Grid()->oSites(),{ + for (int mu = 0; mu < Nd; 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();