return l.elem() * r.elem();
typedef typename BinaryReturn<RComplex<T1>, RComplex<T2>, OpMultiply>::Type_t  Ret_t;  return Ret_t(l.real()*r.real() - l.imag()*r.imag(), l.real()*r.imag() + l.imag()*r.real());
traceMultiply(const RScalar<T1>& l, const RScalar<T2>& r){//  return traceMultiply(l.elem(), r.elem());  /*! NOTE: removed trace here !!!!!  */  return l.elem() * r.elem();}