d.elem(i) = l.elem() * r.elem(i);
{ d.elem(i) = l.elem(i,0) * r.elem(0); for(int j=1; j < N; ++j) d.elem(i) += l.elem(i,j) * r.elem(j); }
operator*(const PScalar<T1>& l, const PVector<T2,N,C>& r){  typename BinaryReturn<PScalar<T1>, PVector<T2,N,C>, OpMultiply>::Type_t  d;  for(int i=0; i < N; ++i)    d.elem(i) = l.elem() * r.elem(i);  return d;}