zero_rep(d.real()); d.imag() = s1.elem();
d.real() = -s1.imag(); d.imag() =  s1.real();
timesI(const RScalar<T>& s1){  typename UnaryReturn<RScalar<T>, FnTimesI>::Type_t  d;  zero_rep(d.real());  d.imag() = s1.elem();  return d;}