portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-11-04 14:04:32 +00:00
Code Releases Activity

Dont alias lhs and rhs in a where statement

Browse Source
This commit is contained in:
paboyle
2018-03-06 04:14:13 -08:00
parent 3cb1b545d0
commit 40699221e2
1 changed files with 11 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

19
lib/qcd/action/gauge/GaugeImplementations.h
View File

@@ -102,11 +102,13 @@ public:
Lattice<iScalar<vInteger>> coor(grid); Lattice<iScalar<vInteger>> coor(grid);
LatticeCoordinate(coor, mu); LatticeCoordinate(coor, mu);
GaugeLinkField tmp(grid); GaugeLinkField tmp1(grid);
tmp = adj(Link); GaugeLinkField tmp2(grid);
tmp = where(coor == Lmu, conjugate(tmp), tmp); tmp1 = adj(Link);
return Cshift(tmp, mu, -1); // moves towards positive mu tmp2 = where(coor == Lmu, conjugate(tmp1), tmp1);
return Cshift(tmp2, mu, -1); // moves towards positive mu
} }
static inline GaugeLinkField static inline GaugeLinkField
CovShiftIdentityForward(const GaugeLinkField &Link, int mu) { CovShiftIdentityForward(const GaugeLinkField &Link, int mu) {
return Link; return Link;
@@ -119,10 +121,11 @@ public:
Lattice<iScalar<vInteger>> coor(grid); Lattice<iScalar<vInteger>> coor(grid);
LatticeCoordinate(coor, mu); LatticeCoordinate(coor, mu);
GaugeLinkField tmp(grid); GaugeLinkField tmp1(grid);
tmp = Cshift(Link, mu, 1); GaugeLinkField tmp2(grid);
tmp = where(coor == Lmu, conjugate(tmp), tmp); tmp1 = Cshift(Link, mu, 1);
return tmp; tmp2 = where(coor == Lmu, conjugate(tmp1), tmp1);
return tmp2;
} }
static inline bool isPeriodicGaugeField(void) { return false; } static inline bool isPeriodicGaugeField(void) { return false; }