const OpAddAssign &op,
const OpSubtractAssign &op,
void evaluate( OLattice< TVec > &d,	       const OpAddAssign &op,	       const QDPExpr< 	       UnaryNode<OpIdentity,	       Reference< QDPType< TVec, OLattice< TVec > > > >,	       OLattice< TVec > > &rhs,	       const OrderedSubset& s){  const OLattice< TVec >& x = static_cast< const OLattice<TVec >&>(rhs.expression().child()); #ifdef DEBUG_BLAS  QDPIO::cout << "BJ: v += v" << endl;#endif  int n_3vec = (s.end() - s.start()+1)*Ns;  REAL *xptr = (REAL *)(&x.elem(s.start()).elem(0).elem(0).real().elem());  REAL *yptr = (REAL *)(&d.elem(s.start()).elem(0).elem(0).real().elem());  REAL one = 1;  vaxpy3(yptr, &one, yptr, xptr,n_3vec);}