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