#include using namespace Grid; using namespace Grid::QCD; int main (int argc, char ** argv) { Grid_init(&argc,&argv); std::vector latt_size = GridDefaultLatt(); int Nd = latt_size.size(); std::vector simd_layout = GridDefaultSimd(Nd,vComplex::Nsimd()); std::vector mpi_layout = GridDefaultMpi(); std::vector mask(Nd,1); mask[0]=0; GridCartesian Fine (latt_size,simd_layout,mpi_layout); GridRedBlackCartesian RBFine(latt_size,simd_layout,mpi_layout,mask,1); GridParallelRNG FineRNG(&Fine); FineRNG.SeedRandomDevice(); LatticeComplex U(&Fine); LatticeComplex ShiftU(&Fine); LatticeComplex rbShiftU(&Fine); LatticeComplex Ue(&RBFine); LatticeComplex Uo(&RBFine); LatticeComplex ShiftUe(&RBFine); LatticeComplex ShiftUo(&RBFine); LatticeComplex lex(&Fine); lex=zero; Integer stride =1; { double nrm; LatticeComplex coor(&Fine); for(int d=0;d coor(4); std::cout << "Checking the non-checkerboard shift"< scoor(coor); scoor[dir] = (scoor[dir]+shift)%latt_size[dir]; Integer slex = scoor[0] + latt_size[0]*scoor[1] + latt_size[0]*latt_size[1]*scoor[2] + latt_size[0]*latt_size[1]*latt_size[2]*scoor[3]; Complex scm(slex); double nrm = abs(scm-cm()()()); std::vector peer(4); Complex ctmp = cm; Integer index=real(ctmp); Fine.CoorFromIndex(peer,index,latt_size); if (nrm > 0){ std::cerr<<"FAIL shift "<< shift<<" in dir "<< dir <<" ["< scoor(coor); scoor[dir] = (scoor[dir]+shift)%latt_size[dir]; Integer slex = scoor[0] + latt_size[0]*scoor[1] + latt_size[0]*latt_size[1]*scoor[2] + latt_size[0]*latt_size[1]*latt_size[2]*scoor[3]; Complex scm(slex); nrm = abs(scm-cm()()()); std::vector peer(4); Complex ctmp=cm; Integer index=real(ctmp); Fine.CoorFromIndex(peer,index,latt_size); if (nrm > 0){ std::cerr<<"FAIL shift "<< shift<<" in dir "<< dir <<" ["<