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

Looping construct changes

Browse Source
This commit is contained in:
Peter Boyle 2019-06-15 08:35:57 +01:00
parent 09720c40cd
commit e8b78f596e
1 changed files with 7 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
Grid/qcd/utils/SUn.h
View File

@ -225,7 +225,9 @@ public:
auto subgroup_v = subgroup.View(); auto subgroup_v = subgroup.View();
auto source_v = source.View(); auto source_v = source.View();
auto Determinant_v = Determinant.View(); auto Determinant_v = Determinant.View();
thread_loop( (int ss = 0; ss < grid->oSites(); ss++) ,{
thread_for(ss, grid->oSites(), {
subgroup_v[ss]()()(0, 0) = source_v[ss]()()(i0, i0); subgroup_v[ss]()()(0, 0) = source_v[ss]()()(i0, i0);
subgroup_v[ss]()()(0, 1) = source_v[ss]()()(i0, i1); subgroup_v[ss]()()(0, 1) = source_v[ss]()()(i0, i1);
subgroup_v[ss]()()(1, 0) = source_v[ss]()()(i1, i0); subgroup_v[ss]()()(1, 0) = source_v[ss]()()(i1, i0);
@ -257,7 +259,8 @@ public:
dest = 1.0; // start out with identity dest = 1.0; // start out with identity
auto dest_v = dest.View(); auto dest_v = dest.View();
auto subgroup_v = subgroup.View(); auto subgroup_v = subgroup.View();
thread_loop( (int ss = 0; ss < grid->oSites(); ss++) ,{ thread_for(ss, grid->oSites(),
{
dest_v[ss]()()(i0, i0) = subgroup_v[ss]()()(0, 0); dest_v[ss]()()(i0, i0) = subgroup_v[ss]()()(0, 0);
dest_v[ss]()()(i0, i1) = subgroup_v[ss]()()(0, 1); dest_v[ss]()()(i0, i1) = subgroup_v[ss]()()(0, 1);
dest_v[ss]()()(i1, i0) = subgroup_v[ss]()()(1, 0); dest_v[ss]()()(i1, i0) = subgroup_v[ss]()()(1, 0);
@ -517,8 +520,8 @@ public:
a[2] = a123mag * sin_theta * sin(phi); a[2] = a123mag * sin_theta * sin(phi);
a[3] = a123mag * cos_theta; a[3] = a123mag * cos_theta;
ua = toComplex(a[0]) * ident + toComplex(a[1]) * pauli1 + ua = toComplex(a[0]) * ident + toComplex(a[1]) * pauli1 +
toComplex(a[2]) * pauli2 + toComplex(a[3]) * pauli3; toComplex(a[2]) * pauli2 + toComplex(a[3]) * pauli3;
b = 1.0; b = 1.0;
b = where(wheremask, uinv * ua, b); b = where(wheremask, uinv * ua, b);