const LatticeFermion& x = static_cast<const LatticeFermion &>(rhs.expression().left()); const Real& a = static_cast<const Real &> (rhs.expression().right());  REAL ar = a.elem().elem().elem().elem().elem();
const OLattice< TVec >& x = static_cast<const OLattice< TVec > &>(rhs.expression().left()); const OScalar< TScal >& a = static_cast<const OScalar< TScal > &> (rhs.expression().right());   REAL ar = -( a.elem().elem().elem().elem().elem());
void evaluate(OLattice< TVec >& d, 	      const OpAddAssign& op, 	      const QDPExpr<BinaryNode<OpMultiply, 	      Reference< QDPType< TVec, OLattice< TVec > > >,	      Reference< QDPType< TScal, OScalar < TScal > > > >,	      OLattice< TVec > > &rhs,	      const OrderedSubset& s){#ifdef DEBUG_BLAS  QDPIO::cout << "y += x*a" << endl;#endif  const LatticeFermion& x = static_cast<const LatticeFermion &>(rhs.expression().left());  const Real& a = static_cast<const Real &> (rhs.expression().right());    REAL ar = a.elem().elem().elem().elem().elem();  REAL* aptr = &ar;  REAL* xptr = (REAL *)&(x.elem(s.start()).elem(0).elem(0).real().elem());  REAL* yptr = &(d.elem(s.start()).elem(0).elem(0).real().elem());  // cout << "Specialised axpy a ="<< ar << endl;    int n_3vec = (s.end()-s.start()+1)*Ns;  vaxpy3(yptr, aptr, xptr, yptr, n_3vec);}