cout << "PSc<M>*adj(PSc<M>)" << endl;
cout << "M*adj(M)" << endl;
multiplyAdj(const PScalar<PColorMatrix<RComplexFloat,3> >& l, 	    const PScalar<PColorMatrix<RComplexFloat,3> >& r){  BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,     PScalar<PColorMatrix<RComplexFloat,3> >, OpMultiplyAdj>::Type_t  d;#if defined(QDP_SCALARSITE_DEBUG)  cout << "PSc<M>*adj(PSc<M>)" << endl;#endif  _inline_generic_mult_su3_na(l.elem(),r.elem(),d.elem());  return d;}