diff --git a/Grid/qcd/action/momentum/DDHMCfilter.h b/Grid/qcd/action/momentum/DDHMCfilter.h index c96351e2..db2ac156 100644 --- a/Grid/qcd/action/momentum/DDHMCfilter.h +++ b/Grid/qcd/action/momentum/DDHMCfilter.h @@ -70,32 +70,25 @@ struct DDHMCFilter: public MomentumFilterBase Coordinate Global=grid->GlobalDimensions(); LatticeInteger coor(grid); for(int mu=0;mu0) ) { LatticeCoordinate(coor,mu); auto P_mu = PeekIndex(P, mu); P_mu = where(mod(coor,Block[mu])==Integer(Block[mu]-1),zz,P_mu); - /* if(Width>=2) { - P_mu = where(mod(coor,Block[mu])==Integer(0),zz,P_mu); - P_mu = where(mod(coor,Block[mu])==Integer(Block[mu]-2),zz,P_mu); - }*/ PokeIndex(P, P_mu, mu); for(int nu=0;nu(P, nu); - if ( Width>=1 ){ - P_nu = where(mod(coor,Block[mu])==Integer(Block[mu]-1),zz,P_nu); - P_nu = where(mod(coor,Block[mu])==Integer(0) ,zz,P_nu); - } - /* if ( Width>=2 ){ - P_nu = where(mod(coor,Block[mu])==Integer(Block[mu]-2),zz,P_nu); - P_nu = where(mod(coor,Block[mu])==Integer(1) ,zz,P_nu); - }*/ + P_nu = where(mod(coor,Block[mu])==Integer(Block[mu]-1),zz,P_nu); + P_nu = where(mod(coor,Block[mu])==Integer(0) ,zz,P_nu); PokeIndex(P, P_nu, nu); } } + } + } } };