PSpinMatrix& operator=(const PSpinMatrix<T1,N>& rhs)
PSpinMatrix& operator=(const PScalar<T1>& rhs)
  PSpinMatrix& operator=(const PSpinMatrix<T1,N>& rhs)     {      assign(rhs);      return *this;    }