typedef BinaryReturn<ILatticeFloat, ILatticeFloat, OpMultiply>::Type_t  Ret_t;
operator*(const PMatrix<RComplexFloat,3,PColorMatrix>& l, 	  const PMatrix<RComplexFloat,3,PColorMatrix>& r){//  cout << "M*M" << endl;  BinaryReturn<PMatrix<RComplexFloat,3,PColorMatrix>,     PMatrix<RComplexFloat,3,PColorMatrix>, OpMultiply>::Type_t  d;  REAL32 *dd = (REAL32*)&d;  REAL32 *ll = (REAL32*)&l;  REAL32 *rr = (REAL32*)&r;  _inline_ssevec_mult_su3_nn(dd,ll,rr,0);  _inline_ssevec_mult_su3_nn(dd,ll,rr,1);  _inline_ssevec_mult_su3_nn(dd,ll,rr,2);  return d;}