utmatrix_t utmatrix_t::operator+(const utmatrix_t& o) {
  utmatrix_t product(size());

  for (int i = 0; i < size()*(size()+1)/2; i++)
  product.data[i]=data[i]+o.data[i];

  return product;
}