Tritset::Reference Tritset::operator[](const int i) {

  if (i >= sizetrits) {
    int tmp = (int) (ceil(i * 2 / (float) 8 / (float) sizeof(unsigned)));

    if (tmp <= sizearr) {
      return Reference(memory, i);
    } else {
      unsigned *newMemory;
      newMemory=new unsigned[tmp];
      memset(newMemory, 0, tmp * sizeof(unsigned));
      memcpy(newMemory, memory, sizearr * sizeof(unsigned));
      sizearr = tmp;
      sizetrits = i;
      delete[] memory;
      memory = newMemory;
      return Reference(newMemory, i);
    }
  }

  return Reference(memory, i);
}