CubeStorage.H

/* CubeStorage.H */
/* $Id: CubeStorage.H,v 1.13 2008/04/11 11:49:38 vlad Exp $ */
#ifndef __CubeStorage_H
#define __CubeStorage_H

#include <errno.h>
#include <fcntl.h>
#include <unistd.h>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <stdarg.h>
#include <sys/stat.h>
#include <sys/types.h>

#ifndef GEVSP
#include <mix/Exceptions.h>
#include <mix/Types.h>
#else
/* Replacement for UNIVERS exceptions stuff */

#define Throw(x)    fprintf(stderr,"%s:%d: %s\n",__FILE__, __LINE__,x),abort()
#define ERROR__NON_POSITIVE_NUMBER  "ERROR__NON_POSITIVE_NUMBER"
#define ERROR__DA_OUT_OF_RANGE      "ERROR__DA_OUT_OF_RANGE"
#define ERROR__NULL         "ERROR__NULL"
#define ERROR__MISMATCH         "ERROR__MISMATCH"

#endif /* GEVSP */

#include "cust_defs.h"


/*

  #define CuSt_AtomItem if you are sure the T is atomic type (int,
  float, char, pointer etc) to be printable by passing to fprintf(...)

 */


template <class T>
class CubeStorage
{
public:

  CubeStorage (int nR, int nC,
           item_order_t itor = AlaC);
  CubeStorage (int dim_n, int* dims,
           item_order_t itor = AlaC);
  CubeStorage (const CubeStorage& cs,
           item_order_t itor = AlaC);

  virtual ~CubeStorage ();


  CubeStorage&  operator= (const CubeStorage& cs);


  int   dim_number () const;
  int   dim (int dim_i) const;
  void  new_size (int dim_n, int* dims,
          item_order_t itor = AlaC);
  int   lin_index (int i0, ...) const;
  int   lin_index (int* i) const;
  int   items_number (int dim_n, int* dims) const;

  T&    fetch (int* i);
  T&    fetch (int i0, ...);
  T get (int* i) const;
  T get (int i0, ...) const;

  T&    operator() (int* i);
  T&    operator() (int i0, ...);

  T*    c () const;

  T*    f () const;

  void  clean ();

  void  fill (T v);

  virtual void  print_info () const;
  virtual void  print_index (int* i) const;

  virtual void  print_direction (FILE* fp, const char* format);

  /* Dummy methods in case of CuSt_AtomItem is undefined */
  virtual void  print_cell (int* i) const;
  virtual void  print_vector (int d, int* i) const;
  virtual void  print_matrix (int dv, int dh, int* i) const;

#ifndef GEVSP

  virtual void  stream_dir (const char* dir);

  const char*   stream_dir () const;

  virtual void  stream_put (int ii);

  virtual void  stream_get (int ii);
#endif /* GEVSP */

private:

  item_order_t  eItemOrder; 
  T     *pStorage;  
  int       nDim;       
  int       *pDim;      
  /* print setup */
  FILE      *fpLog;
  char      *szFormat;

  char      *szStreamDir;

#ifndef GEVSP

  void  stream_put_header (int fildes);

  void  stream_get_header (int fildes);
#endif /* GEVSP */

};


#include "cust_lifetime.H"
#include "cust_putget.H"
#include "cust_print.H"
#include "cust_info.H"
#ifndef GEVSP
#include "cust_stream.H"
#endif /* GEVSP */
#include "cust_other.H"


#endif /* CubeStorage.H */