elem(i,j) = rhs.elem(i,j);
if (i == j) elem(i,j) = rhs.elem(); else zero_rep(elem(i,j));
  CC& assign(const C<T1,N>& rhs)     {      for(int i=0; i < N; ++i)	for(int j=0; j < N; ++j)	  elem(i,j) = rhs.elem(i,j);      return static_cast<CC&>(*this);    }