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