cust_lifetime.H

/* cust_lifetime.H */
/* $Id: cust_lifetime.H,v 1.5 2007/10/10 11:52:45 vlad Exp $ */
#ifdef __CubeStorage_H


/***********************************************************************
 * Create 2D cube - matrix.
 */
template<class T>
CubeStorage<T>::CubeStorage (int nR, int nC, item_order_t itor)
  : pDim(NULL), /*iCur(NULL),*/ pStorage(NULL),
    fpLog(NULL), szFormat(NULL), szStreamDir(NULL)
{
  int   dims[2];
  dims[0] = nR;
  dims[1] = nC;

  new_size(2, dims, itor);

  fpLog = stdout;
  szFormat = NULL;
}


/***********************************************************************
 * Create any dimension supercube.
 */
template<class T>
CubeStorage<T>::CubeStorage (int dim_n, int* dims, item_order_t itor)
  : pDim(NULL), /*iCur(NULL),*/ pStorage(NULL),
    fpLog(NULL), szFormat(NULL), szStreamDir(NULL)
{
  new_size(dim_n, dims, itor);

  fpLog = stdout;
  szFormat = NULL;
}


/***********************************************************************
 * Create cube as a copy of another one.
 */
template<class T>
CubeStorage<T>::CubeStorage (const CubeStorage<T>& cs, item_order_t itor)
  : pDim(NULL), /*iCur(NULL),*/ pStorage(NULL),
    fpLog(NULL), szFormat(NULL), szStreamDir(NULL)
{
  operator=(cs);

  fpLog = stdout;
  szFormat = NULL;
}


/***********************************************************************
 * Destroy object.
 */
template<class T>
CubeStorage<T>::~CubeStorage ()
{
  /*delete[] iCur;*/
  delete[] pDim;
  delete[] pStorage;
  delete[] szFormat;
  delete[] szStreamDir;
}


/***********************************************************************
 * Copy cs to this supercube erasing previous contents.
 */
template<class T>
CubeStorage<T>&
CubeStorage<T>::operator= (const CubeStorage<T>& cs)
{
  new_size(cs.nDim, cs.pDim, eItemOrder);

  if(cs.eItemOrder == eItemOrder)
    memcpy(pDim, cs.pDim, sizeof(T));
  else /* copy item by item */
    {
      int   dim_i, *i = new int[nDim];
      bool  good = true;

      for(dim_i = 0; dim_i < nDim; ++dim_i)
    i[dim_i] = 0;

      while(good)
    {
      /* copy itself */
      pStorage[lin_index(i)] = cs.pStorage[lin_index(i)];

      /* flush flag of allowed index */
      good = FALSE;

      /* increment index */
      for(dim_i = 0; dim_i < nDim; ++dim_i)
        if(i[dim_i] < pDim[dim_i] - 1)
          {
        ++i[dim_i];
        good = TRUE;
        break;
          }
        else
          i[dim_i] = 0;
    }
    }
}


#endif /* cust_lifetime.H */