portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-04 19:25:56 +01:00
Code Releases Activity

Revert a change

Browse Source
This commit is contained in:
paboyle 2018-03-07 16:54:11 +00:00
parent 4e3458516a
commit 4154fc6f44
1 changed files with 9 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
lib/qcd/utils/CovariantCshift.h
View File

@ -95,13 +95,11 @@ namespace ConjugateBC {
Lattice<iScalar<vInteger> > coor(grid); LatticeCoordinate(coor,mu); Lattice<iScalar<vInteger> > coor(grid); LatticeCoordinate(coor,mu);
Lattice<covariant> tmp1(grid); Lattice<covariant> field_bc = Cshift(field,mu,1);// moves towards negative mu;
Lattice<covariant> tmp2(grid);
tmp1 = Cshift(field,mu,1);// moves towards negative mu; field_bc = where(coor==Lmu,conjugate(field_bc),field_bc);
tmp2 = where(coor==Lmu,conjugate(tmp1),tmp1); // std::cout<<"Gparity::CovCshiftForward mu="<<mu<<std::endl;
return Link*field_bc;
return Link*tmp2;
} }
template<class covariant,class gauge> Lattice<covariant> CovShiftBackward(const Lattice<gauge> &Link, template<class covariant,class gauge> Lattice<covariant> CovShiftBackward(const Lattice<gauge> &Link,
@ -116,13 +114,12 @@ namespace ConjugateBC {
Lattice<iScalar<vInteger> > coor(grid); LatticeCoordinate(coor,mu); Lattice<iScalar<vInteger> > coor(grid); LatticeCoordinate(coor,mu);
Lattice<covariant> tmp1(grid); Lattice<covariant> tmp(grid);
Lattice<covariant> tmp2(grid);
tmp1 = adj(Link)*field; tmp = adj(Link)*field;
tmp2 = where(coor==Lmu,conjugate(tmp1),tmp1); tmp = where(coor==Lmu,conjugate(tmp),tmp);
// std::cout<<"Gparity::CovCshiftBackward mu="<<mu<<std::endl;
return Cshift(tmp2,mu,-1);// moves towards positive mu return Cshift(tmp,mu,-1);// moves towards positive mu
} }
} }