From c18025c0b8b7be8ffa916518ff1d676a2298f0a9 Mon Sep 17 00:00:00 2001
From: Quadro <Peter Boyle>
Date: Tue, 18 May 2021 16:35:39 -0400
Subject: [PATCH] Luscher's filter - will tune this later

---
 Grid/qcd/action/momentum/DDHMCfilter.h | 21 +++++++--------------
 1 file changed, 7 insertions(+), 14 deletions(-)

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<MomentaField>
     Coordinate Global=grid->GlobalDimensions();
     LatticeInteger coor(grid); 
     for(int mu=0;mu<Nd;mu++) {
-      if ( Block[mu] <= Global[mu] ) {
+
+      if ( (Block[mu] <= Global[mu])&&(Block[mu]>0) ) {
 
 	LatticeCoordinate(coor,mu);
 	auto P_mu = PeekIndex<LorentzIndex>(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<LorentzIndex>(P, P_mu, mu);
 
 	for(int nu=0;nu<Nd;nu++){
-	  if ( (mu!=nu) && Block[mu] <= Global[mu]){
+	  if ( mu!=nu ){
 	    auto P_nu = PeekIndex<LorentzIndex>(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<LorentzIndex>(P, P_nu, nu);
 	  }
 	}
+
       }
+      
     }
   }
 };