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