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

View related changes

Browse Source
This commit is contained in:
paboyle 2018-03-04 16:33:33 +00:00
parent c03423250f
commit 8824a54269
1 changed files with 7 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
lib/qcd/action/gauge/GaugeImplTypes.h
View File

@ -74,11 +74,11 @@ public:
// Move this elsewhere? FIXME
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 W_v = W.View();
thread_loop ( (auto ss = 0; ss < U.Grid()->oSites(); ss++),{
U[ss]._internal[mu] =
U[ss]._internal[mu] + W[ss]._internal;
U_v[ss](mu) = U_v[ss](mu) + W_v[ss]();
});
}
@ -101,9 +101,11 @@ public:
//static std::chrono::duration<double> diff;
//auto start = std::chrono::high_resolution_clock::now();
auto U_v = U.View();
auto P_v = P.View();
thread_loop( (int ss=0;ss<P.Grid()->oSites();ss++),{
for (int mu = 0; mu < Nd; mu++)
U[ss]._internal[mu] = ProjectOnGroup(Exponentiate(P[ss]._internal[mu], ep, Nexp) * U[ss]._internal[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();