multi1d<LatticeFermion> tmp(N5) ;  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) ; } CompbineHoppingVectors(chi,tmp,isign,0) ; break ; }
applyOffDiag(chi, psi, isign, 0);
  void evenOddLinOp(multi1d<LatticeFermion>& chi, 		    const multi1d<LatticeFermion>& psi, 		    enum PlusMinus isign) const    {      multi1d<LatticeFermion> tmp(N5) ;      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) ;	}	CompbineHoppingVectors(chi,tmp,isign,0) ;	break ;      }    }