CompbineHoppingVectors(tmp,psi,0) ; for(int s(0);s<N5;s++) { D.apply(chi[s],tmp[s],isign,0); chi[s][rb[0]] *= (-0.5);
switch ( isign ){ case PLUS: CompbineHoppingVectors(tmp,psi,isign,1) ; for(int s(0);s<N5;s++) { D.apply(chi[s],tmp[s],isign,0); chi[s][rb[0]] *= (-0.5); } break; case MINUS: for(int s(0);s<N5;s++){ D.apply(tmp[s],psi[s],isign,0); tmp[s][rb[0]] *= (-0.5) ;
  void evenOddLinOp(multi1d<LatticeFermion>& chi, 		    const multi1d<LatticeFermion>& psi, 		    enum PlusMinus isign) const    {      multi1d<LatticeFermion> tmp(N5) ;      CompbineHoppingVectors(tmp,psi,0) ;      for(int s(0);s<N5;s++)	{	  D.apply(chi[s],tmp[s],isign,0);	  chi[s][rb[0]] *= (-0.5);	}    }