cell.cpp

Go to the documentation of this file.

#include <cell.h>

cell::~cell()
{
  delete[] vi;
  delete[] fi;
  delete[] count;
  delete[] index;
}

cell::cell()
  : N(0)
{
  W = WW = 0;
  count = 0;
  index = 0;
  vi = 0;
  fi = 0;
}

void cell::reset(uintc N_)
{
  if (count)
    delete[] count;
  if (vi)
    delete[] vi;
  if (index)
    delete[] index;
  if (fi)
    delete[] fi;

  N = N_;
  W = (uint)( sqrt((double)N) )+1;
  WW = W*W;
  count = new uint[WW];
  index = new uint[WW+1];
  index[WW] = N; // Sentinal value
  //  suchthat  index[k+1]-index[k] = count[k]
  vi = new uint[N];

  fi = new uint[N];

  for (uint i=0; i<WW; ++i)
    count[i]=0;
}

cell::cell(uintc N_)
{
  reset(N_);
}