diff --git a/Grid/qcd/action/pseudofermion/TwoFlavourRatio.h b/Grid/qcd/action/pseudofermion/TwoFlavourRatio.h index e0e474bf..c4bb2fd4 100644 --- a/Grid/qcd/action/pseudofermion/TwoFlavourRatio.h +++ b/Grid/qcd/action/pseudofermion/TwoFlavourRatio.h @@ -97,21 +97,9 @@ public: tmp = Zero(); ActionSolver(MdagMOp,Phi,tmp); // (VdagV)^-1 Mdag eta = V^-1 Vdag^-1 Mdag eta NumOp.M(tmp,Phi); // Vdag^-1 Mdag eta -#define FILTER -#ifdef FILTER - Integer OrthogDir=0; - Integer plane=0; - if ( getenv("DIR") ) OrthogDir = atoi(getenv("DIR")); - if ( getenv("COOR") ) plane = atoi(getenv("COOR")); - std::cout << " *** PseudoFermion FILTER DIR " < > coor(NumOp.FermionGrid()); - LatticeCoordinate(coor,OrthogDir); - tmp = Zero(); - Phi = where(coor==plane,Phi,tmp); -#endif - + Phi=Phi*scale; - + std::cout << GridLogMessage<<" TwoFlavourRatio Expect action to be "< > coor(NumOp.GaugeGrid()); - - LatticeCoordinate(coor,mu); - int L = NumOp.GaugeGrid()->FullDimensions()[mu]; - for (Integer p=0;p