source: Sophya/trunk/SophyaLib/BaseTools/ppfbinstream.cc@ 2476

#include "machdefs.h"
#include <stdio.h>
#include <sys/types.h>
#include <time.h>
#include "ppfbinstream.h"
#include "pexceptions.h"
#include <iostream>


/*
//  strptime n'est pas defini sous Linux  - Reza Mars 2000
#if defined(OS_LINUX) || defined(OS_MACOSX)
extern "C" {
char *strptime(const char *buf, const char *format, struct tm *tm);
}
#endif
*/

#define MAXTAGLEN_V2 255


static inline void bswap8(void* p)
{
  uint_8 tmp = *(uint_8*)p;
  *(uint_8*)p = ((tmp >> (7*8)) & 0x000000FF) | 
                ((tmp >> (5*8)) & 0x0000FF00) |
                ((tmp >> (3*8)) & 0x00FF0000) |
                ((tmp >> (1*8)) & 0xFF000000) |
                ((tmp & 0xFF000000) << (1*8))  |
                ((tmp & 0x00FF0000) << (3*8))  |
                ((tmp & 0x0000FF00) << (5*8))  |
                ((tmp & 0x000000FF) << (7*8));
}

static inline void bswap4(void* p)
{
  uint_4 tmp = *(uint_4*)p;
  *(uint_4*)p = ((tmp >> 24) & 0x000000FF) | 
                ((tmp >> 8)  & 0x0000FF00) |
                ((tmp & 0x0000FF00) << 8)  |
                ((tmp & 0x000000FF) << 24);
}

static inline void bswap2(void* p)
{
  uint_2 tmp = *(uint_2*)p;
  *(uint_2*)p = ((tmp >> 8) & 0x00FF) | 
                ((tmp & 0x00FF) << 8);
}


PPFBinaryInputStream::PPFBinaryInputStream(RawInOutStream * is, bool ad, bool scan)
{
  s = is;
  _ads = ad;
  Init(scan);
}

//++
PPFBinaryInputStream::PPFBinaryInputStream(string const& flnm, bool scan)
//
//      Constructeur. Ouvre le fichier.
//--
{
  s = new RawInFileStream(flnm.c_str());
  _ads = true;
  Init(scan);
}

PPFBinaryInputStream::~PPFBinaryInputStream()
{
  if (_ads && (s!=NULL) ) delete s;
}

void
PPFBinaryInputStream::Init(bool scan)
{
  // Read and check header

  char rbuf[36];
  GetRawBytes(rbuf, 32);
  if (strncmp(rbuf,"SOS-SOPHYA-PPersistFile", 23) != 0)  {
    throw FileFormatExc("PPFBinaryInputStream::PPFBinaryInputStream  bad header");
  }
  rbuf[32] = '\0';
  version = atoi(rbuf+25);
  if (version < 2) {
    cerr << "PPFBinaryInputStream::PPFBinaryInputStream(" << FileName() << ") Version(=" << version 
         << ") < 2 not supported !" << endl;
    throw FileFormatExc("PPFBinaryInputStream::PPFBinaryInputStream() - Unsupported (Old) Version");
  }
  // read endianness
  GetRawBytes(rbuf, 32);
  if (strncmp(rbuf,"BIG-ENDIAN",10) == 0)
    bigEndian = true;
  else if (strncmp(rbuf,"LITTLE-ENDIAN",13) == 0)
    bigEndian = false;
  else {
    throw FileFormatExc("PPFBinaryInputStream::PPFBinaryInputStream  bad header - endianness");
  }

  // read creation date
  GetRawBytes(rbuf, 32);
  rbuf[32] = '\0'; 
  struct tm tm;
  /*  #if !(defined(__MWERKS__) || defined(OS_MACOSX))  RZ-DEL */
  strptime(rbuf,"%d/%m/%Y %H:%M:%S GMT",&tm);
  /*  #else
  sscanf(rbuf,"%2d/%2d/%4d %2d:%2d:%2d GMT",&tm.tm_mday,&tm.tm_mon,&tm.tm_year,
                                            &tm.tm_hour,&tm.tm_min,&tm.tm_sec);

  tm.tm_mon  --;
  tm.tm_year -= 1900;
  #endif    RZ-DEL */

  creationdate = mktime(&tm);
  seqread = true;   // To flag non sequential reads
  if (scan && s->isSeekable()) { 
    if (Version() >= 3) ReadNameTagTable();
    else ReadNameTagTableV2();
  }
}





string 
PPFBinaryInputStream::CreationDateStr()
{
  time_t cdt = CreationDate();
  string cdate = ctime(&cdt);
  return(cdate);
}

void
PPFBinaryInputStream::ReadNameTagTable()
{
  unsigned char ppstype;
  int_8 debut;
  debut = s->tellg();

  s->seekg(-(sizeof(int_8)+1), ios::end);

  GetTypeTag(ppstype);
  if (ppstype != PPS_EOF)
    throw FileFormatExc("PPFBinaryInputStream::ReadNameTagTable()  Corrupted file, no EOF tag at end of file");

  // Lecture position NameTagTable
  int_8 pos;
  GetRawI8(pos);
  if (pos < 0) {  // no tags
    s->seekg(debut);
    return;
  }
  
  // On se positionne au debut du NameTagTable
  s->seekg(pos);
  GetTypeTag(ppstype);
  if (ppstype != PPS_NAMETAG_TABLE)
    throw FileFormatExc("PPFBinaryInputStream::ReadNameTagTable()  Corrupted file PPS_NAMETAG_TABLE not found");
  uint_8 ttsz,it;
  GetU8(ttsz);
  for(it=0; it<ttsz; it++) {
    int_8 tpos;
    string tname;
    GetI8(tpos);
    GetStr(tname);
    tags[tname] = tpos;
  }
  // On revient au debut du float, juste apres l'entete
  s->seekg(debut);
}

void
PPFBinaryInputStream::ReadNameTagTableV2()
{
  // On cherche la liste des tags, a la fin du fichier

  unsigned char ppstype;
  int_8 debut;
  debut = s->tellg();

  s->seekg(-(sizeof(int_8)+1), ios::end);

  GetTypeTag(ppstype);
  if (ppstype != PPS_EOF)
    throw FileFormatExc("PPFBinaryInputStream::ReadNameTagTableV2() Corrupted file, no eof entry at end of file");

  int_8 pos;
  GetRawI8(pos);
  if (pos < 0) {  // no tags
    s->seekg(debut);
    return;
  }

  char buffer[MAXTAGLEN_V2+1];
  s->seekg(pos);
  while (true) {
    GetTypeTag(ppstype);
    if (ppstype == PPS_EOF) break;
    
    if (ppstype != PPS_NAMETAG_TABLE)
      throw FileFormatExc("PPFBinaryInputStream::ReadNameTagTableV2() Corrupted file, bad tag entry");

    GetRawI8(pos);
    int_4 len;
    GetRawI4(len);
    if (len > MAXTAGLEN_V2) 
      throw FileFormatExc("PPFBinaryInputStream::ReadNameTagTableV2() Corrupted file, tag name too long");
    GetRawBytes(buffer, len);
    buffer[len] = '\0';

    tags[buffer] = pos;
  }
  s->seekg(debut);
}


int
PPFBinaryInputStream::NbNameTags()
{
  return tags.size();
}

bool 
PPFBinaryInputStream::GotoPositionTag(int_8 pos)
{
  s->seekg(pos);
  int_8 tpos;
  GetRawI8(tpos);
  if (tpos != pos) 
    throw FileFormatExc("PPFBinaryInputStream::GotoPositionTag() - Wrong tag position!");
  return true;
}

bool 
PPFBinaryInputStream::GotoNameTag(string const& name)
{
  map<string, int_8>::iterator i = tags.find(name);
  if (i == tags.end()) 
    return false;
      //    throw NotFoundExc("PPFBinaryInputStream::GotoNameTag   tag not found");
  s->seekg((*i).second);
  seqread = false;
  //  objList.clear(); $CHECK$ EA 171199 Ne pas faire ? Reza 03/2000 ?
  return(true); 
}


bool
PPFBinaryInputStream::GotoNameTagNum(int itag)
{
  if (itag<0 || itag >= (int)tags.size()) return false;
  map<string, int_8>::iterator i = tags.begin();
  for (int j=0; j<itag; j++) i++;
  s->seekg((*i).second);
  seqread = false;
  // objList.clear();  $CHECK$ EA 171199  Ne pas faire ? Reza 03/2000 ?
  return(true);
}

string
PPFBinaryInputStream::GetTagName(int itag)
{
  if (itag<0 || itag >= (int)tags.size()) return "";
  map<string, int_8>::iterator i = tags.begin();
  for (int j=0; j<itag; j++) i++;
  return((*i).first);
}


static vector<string> * ret_tag_names = NULL;
vector<string> const &
PPFBinaryInputStream::GetNameTags()
{
if (ret_tag_names) delete ret_tag_names;
ret_tag_names = new vector<string> ;
map<string, int_8>::iterator i;
for(i=tags.begin(); i!=tags.end(); i++) ret_tag_names->push_back((*i).first);
return(*ret_tag_names);  
}

bool
PPFBinaryInputStream::SkipToNextObject()
{
  // A FAIRE NOV 2003 - REZA
  return true;
}

bool
PPFBinaryInputStream::SkipNextItem()
{
  // A FAIRE NOV 2003 - REZA
  return true;
}

char
PPFBinaryInputStream::NextItemTag(short datasz, size_t sz)
{
  datasz = 0;
  sz = 0;
  return 0;
}

void 
PPFBinaryInputStream::SkipItem(bool fgrdt, char itag)
{
  return;
}

//++
// void PPFBinaryInputStream::GetByte(char& c)
// void PPFBinaryInputStream::GetBytes(void* ptr, size_t bytes)
// void PPFBinaryInputStream::GetR4   (r_4& result)
// void PPFBinaryInputStream::GetR4s  (r_4* tab, size_t n)
// void PPFBinaryInputStream::GetR8   (r_8& result)
// void PPFBinaryInputStream::GetR8s  (r_8* tab, size_t n)
// void PPFBinaryInputStream::GetI2   (int_2& result)
// void PPFBinaryInputStream::GetI2s  (int_2* tab, size_t n)
// void PPFBinaryInputStream::GetU2   (uint_2& result)
// void PPFBinaryInputStream::GetU2s  (uint_2* tab, size_t n)
// void PPFBinaryInputStream::GetI4   (int_4& result)
// void PPFBinaryInputStream::GetI4s  (int_4* tab, size_t n)
// void PPFBinaryInputStream::GetU4   (uint_4& result)
// void PPFBinaryInputStream::GetU4s  (uint_4* tab, size_t n)
// void PPFBinaryInputStream::GetI8   (int_8& result)
// void PPFBinaryInputStream::GetI8s  (int_8* tab, size_t n)
// void PPFBinaryInputStream::GetU8   (uint_8& result)
// void PPFBinaryInputStream::GetU8s  (uint_8* tab, size_t n)
//      Lecture de données portables depuis le fichier PPersist. Pour chaque type
//      de données, on peut lire une valeur, ou un tableau de valeurs.
// void PPFBinaryInputStream::GetLine(char* ptr, size_t len)
//      Lecture d'une ligne de texte depuis le fichier PPersist.
//--


void
PPFBinaryInputStream::GetTypeTag(unsigned char& c) 
{
  int_8 tpos;
  c = PPS_NAMETAG_MARK;
  while ( (c == PPS_NAMETAG_MARK) || (c == PPS_POSTAG_MARK) ) {
    GetRawUByte(c);
    if (c == PPS_POSTAG_MARK)  GetRawI8(tpos);
  } 
  //  while (c == PPS_NAMETAG_MARK) {    Il ne faut plus faire ca !
  //    objList.clear();             $CHECK$ Reza 03/2000
  //    GetRawByte(c);
  //  }
}


void
PPFBinaryInputStream::GetRawByte(char& c)
{
  GetRawBytes(&c, 1);
}

void
PPFBinaryInputStream::GetRawUByte(unsigned char& c)
{
  GetRawBytes(&c, 1);
}

void
PPFBinaryInputStream::GetRawBytes(void* ptr, size_t bytes)
{
  s->read((char*)ptr, bytes);
}

void
PPFBinaryInputStream::GetRawI2   (int_2& result)
{
  GetRawBytes(&result, sizeof(int_2));
  if (bigEndian != IS_BIG_ENDIAN)
    bswap2(&result);
}

void
PPFBinaryInputStream::GetRawI4   (int_4& result)
{
  GetRawBytes(&result, sizeof(int_4));
  if (bigEndian != IS_BIG_ENDIAN)
    bswap4(&result);
}

void
PPFBinaryInputStream::GetRawI8   (int_8& result)
{
  GetRawBytes(&result, sizeof(int_8));
  if (bigEndian != IS_BIG_ENDIAN)
    bswap8(&result);
}

void
PPFBinaryInputStream::GetRawU8   (uint_8& result)
{
  GetRawBytes(&result, sizeof(uint_8));
  if (bigEndian != IS_BIG_ENDIAN)
    bswap8(&result);
}

void
PPFBinaryInputStream::CheckTag(short datasz, short datatype)
// datatype = PPS_DATATYPE_CHAR or PPS_DATATYPE_FLOAT or PPS_DATATYPE_INTEGER or PPS_DATATYPE_UNSIGNED
{
  unsigned char ppstype;
  GetTypeTag(ppstype);
  if (ppstype != PPS_SIMPLE + datasz + datatype)
    throw FileFormatExc("PPFBinaryInputStream::CheckTag   bad type in ppersist file");
}

void
PPFBinaryInputStream::CheckArrayTag(short datasz, size_t sz, short datatype)
// datatype = PPS_DATATYPE_CHAR or PPS_DATATYPE_FLOAT or PPS_DATATYPE_INTEGER or PPS_DATATYPE_UNSIGNED
{
  unsigned char ppstype;
  GetTypeTag(ppstype);
  size_t filesz;
  if (sz <= 0x7fffffff) {
    if (ppstype != PPS_SIMPLE_ARRAY4 + datasz + datatype)
      throw FileFormatExc("PPFBinaryInputStream::CheckArrayTag   bad type in ppersist file");
    int_4 ff;
    GetRawI4(ff); filesz=ff;
  } else {
    if (ppstype != PPS_SIMPLE_ARRAY8 + datasz + datatype)
      throw FileFormatExc("PPFBinaryInputStream::CheckArrayTag   bad type in ppersist file");
    uint_8 ff;
    GetRawU8(ff); filesz=ff;
  }
  if (filesz != sz)
    throw FileFormatExc("PPFBinaryInputStream::CheckArrayTag   bad array size in ppersist file");
}

void
PPFBinaryInputStream::GetByte(char& c)
{
  CheckTag(1,PPS_DATATYPE_CHAR);
  GetRawBytes(&c, 1);
}


void
PPFBinaryInputStream::GetBytes(void* ptr, size_t bytes)
{
  CheckArrayTag(1, bytes, PPS_DATATYPE_CHAR);
  GetRawBytes(ptr, bytes);
}
void
PPFBinaryInputStream::GetR4   (r_4& result)
{
  CheckTag(4,PPS_DATATYPE_FLOAT);
  GetRawBytes(&result, sizeof(r_4));
  if (bigEndian != IS_BIG_ENDIAN)
    bswap4(&result);
}


void
PPFBinaryInputStream::GetR4s  (r_4* tab, size_t n)
{
  CheckArrayTag(4,n,PPS_DATATYPE_FLOAT);
  GetRawBytes(tab, n*sizeof(r_4));
  if (bigEndian == IS_BIG_ENDIAN) return;

  for (unsigned int i=0; i<n; i++)
    bswap4(tab+i);

  return;
}

void
PPFBinaryInputStream::GetR8   (r_8& result)
{
  CheckTag(8,PPS_DATATYPE_FLOAT);
  GetRawBytes(&result, sizeof(r_8));
  if (bigEndian != IS_BIG_ENDIAN)
    bswap8(&result);
}

void
PPFBinaryInputStream::GetR8s  (r_8* tab, size_t n)
{
  CheckArrayTag(8,n,PPS_DATATYPE_FLOAT);
  GetRawBytes(tab, n*sizeof(r_8));
  if (bigEndian == IS_BIG_ENDIAN) return;

  for (unsigned int i=0; i<n; i++)
    bswap8(tab+i);

  return;
}

void
PPFBinaryInputStream::GetI1   (int_1& result)
{
  CheckTag(1,PPS_DATATYPE_INTEGER);
  GetRawBytes(&result, sizeof(int_1));
}

void
PPFBinaryInputStream::GetI1s  (int_1* tab, size_t n)
{
  CheckArrayTag(1,n,PPS_DATATYPE_INTEGER);
  GetRawBytes(tab, n*sizeof(int_1));
}

void
PPFBinaryInputStream::GetU1   (uint_1& result)
{
  CheckTag(1,PPS_DATATYPE_UNSIGNED);
  GetRawBytes(&result, sizeof(uint_1));
}

void
PPFBinaryInputStream::GetU1s  (uint_1* tab, size_t n)
{
  CheckArrayTag(1,n,PPS_DATATYPE_UNSIGNED);
  GetRawBytes(tab, n*sizeof(uint_1));
}

void
PPFBinaryInputStream::GetI2   (int_2& result)
{
  CheckTag(2,PPS_DATATYPE_INTEGER);
  GetRawBytes(&result, sizeof(int_2));
  if (bigEndian != IS_BIG_ENDIAN)
    bswap2(&result);
}

void
PPFBinaryInputStream::GetI2s  (int_2* tab, size_t n)
{
  CheckArrayTag(2,n,PPS_DATATYPE_INTEGER);
  GetRawBytes(tab, n*sizeof(int_2));
  if (bigEndian == IS_BIG_ENDIAN) return;

  for (unsigned int i=0; i<n; i++)
    bswap2(tab+i);

  return;
}

void
PPFBinaryInputStream::GetU2   (uint_2& result)
{
  CheckTag(2,PPS_DATATYPE_UNSIGNED);
  GetRawBytes(&result, sizeof(uint_2));
  if (bigEndian != IS_BIG_ENDIAN)
    bswap2(&result);
}

void
PPFBinaryInputStream::GetU2s  (uint_2* tab, size_t n)
{
  CheckArrayTag(2,n,PPS_DATATYPE_UNSIGNED);
  GetRawBytes(tab, n*sizeof(uint_2));
  if (bigEndian == IS_BIG_ENDIAN) return;

  for (unsigned int i=0; i<n; i++)
    bswap2(tab+i);

  return;
}

void
PPFBinaryInputStream::GetI4   (int_4& result)
{
  CheckTag(4,PPS_DATATYPE_INTEGER);
  GetRawBytes(&result, sizeof(int_4));
  if (bigEndian != IS_BIG_ENDIAN)
    bswap4(&result);
}

void
PPFBinaryInputStream::GetI4s  (int_4* tab, size_t n)
{
  CheckArrayTag(4,n,PPS_DATATYPE_INTEGER);
  GetRawBytes(tab, n*sizeof(int_4));
  if (bigEndian == IS_BIG_ENDIAN) return;

  for (unsigned int i=0; i<n; i++)
    bswap4(tab+i);

  return;
}

void
PPFBinaryInputStream::GetU4   (uint_4& result)
{
  CheckTag(4,PPS_DATATYPE_UNSIGNED);
  GetRawBytes(&result, sizeof(uint_4));
  if (bigEndian != IS_BIG_ENDIAN)
    bswap4(&result);
}

void
PPFBinaryInputStream::GetU4s  (uint_4* tab, size_t n)
{
  CheckArrayTag(4,n,PPS_DATATYPE_UNSIGNED);
  GetRawBytes(tab, n*sizeof(uint_4));
  if (bigEndian == IS_BIG_ENDIAN) return;

  for (unsigned int i=0; i<n; i++)
    bswap4(tab+i);

  return;
}


void
PPFBinaryInputStream::GetI8   (int_8& result)
{
  CheckTag(8,PPS_DATATYPE_INTEGER);
  GetRawBytes(&result, sizeof(int_8));
  if (bigEndian != IS_BIG_ENDIAN)
    bswap8(&result);
}

void
PPFBinaryInputStream::GetI8s  (int_8* tab, size_t n)
{
  CheckArrayTag(8,n,PPS_DATATYPE_INTEGER);
  GetRawBytes(tab, n*sizeof(int_8));
  if (bigEndian == IS_BIG_ENDIAN) return;

  for (unsigned int i=0; i<n; i++)
    bswap8(tab+i);

  return;
}

void
PPFBinaryInputStream::GetU8   (uint_8& result)
{
  CheckTag(8,PPS_DATATYPE_UNSIGNED);
  GetRawBytes(&result, sizeof(uint_8));
  if (bigEndian != IS_BIG_ENDIAN)
    bswap8(&result);
}

void
PPFBinaryInputStream::GetU8s  (uint_8* tab, size_t n)
{
  CheckArrayTag(8,n,PPS_DATATYPE_UNSIGNED);
  GetRawBytes(tab, n*sizeof(uint_8));
  if (bigEndian == IS_BIG_ENDIAN) return;

  for (unsigned int i=0; i<n; i++)
    bswap8(tab+i);

  return;
}


void
PPFBinaryInputStream::GetLine(char* ptr, size_t len)
{
  string str;
  GetStr(str);
  strncpy(ptr, str.c_str(), len);
  ptr[len] = '\0';
}

void
PPFBinaryInputStream::GetStr(string& str)
{
  unsigned char ppstype;
  GetTypeTag(ppstype);
  if (ppstype != PPS_STRING) 
    throw FileFormatExc("PPFBinaryInputStream::GetStr   bad type in ppersist file");
  int_4 len;
  GetRawI4(len);
  char * buff = new char[len+1];
  GetRawBytes(buff, len);
  buff[len] = '\0';
  str = buff;
  delete[] buff;
}

void
PPFBinaryInputStream::GetZ4   (complex<r_4>& result)
{
  CheckTag(4,PPS_DATATYPE_COMPLEX);
  r_4 reim[2];
  GetRawBytes(reim, 2*sizeof(r_4));
  if (bigEndian != IS_BIG_ENDIAN) {
    bswap4(reim);
    bswap4(reim+1);
  }
  result = complex<r_4>(reim[0], reim[1]);
}

void
PPFBinaryInputStream::GetZ4s  (complex<r_4>* tab, size_t n)
{  
  CheckArrayTag(4,n,PPS_DATATYPE_COMPLEX);
  GetRawBytes(tab, n*2*sizeof(r_4));
  if (bigEndian == IS_BIG_ENDIAN) return;

  r_4 * p = (r_4 *)tab;
  for (unsigned int i=0; i<n; i++) {
    bswap4(p);  p++;
    bswap4(p);  p++;
  }
  return;
}

void
PPFBinaryInputStream::GetZ8   (complex<r_8>& result)
{
  CheckTag(8,PPS_DATATYPE_COMPLEX);
  r_8 reim[2];
  GetRawBytes(reim, 2*sizeof(r_8));
  if (bigEndian != IS_BIG_ENDIAN) {
    bswap8(reim);
    bswap8(reim+1);
  }
  result = complex<r_8>(reim[0], reim[1]);
}

void
PPFBinaryInputStream::GetZ8s  (complex<r_8>* tab, size_t n)
{
  CheckArrayTag(8,n,PPS_DATATYPE_COMPLEX);
  GetRawBytes(tab, n*2*sizeof(r_8));
  if (bigEndian == IS_BIG_ENDIAN) return;

  r_8 * p = (r_8 *)tab;
  for (unsigned int i=0; i<n; i++) {
    bswap8(p);  p++;
    bswap8(p);  p++;
  }
  return;
}


void
PPFBinaryInputStream::GetPosTagTable(int_8* ptab, size_t sz)
{
  unsigned char ppstype;
  GetTypeTag(ppstype);
  if (ppstype != PPS_POSTAG_TABLE) 
    throw FileFormatExc("PPFBinaryInputStream::GetPosTagTable  bad type in ppersist stream");
  int_4 tsz;
  GetRawI4(tsz);
  if (tsz != sz)
    throw FileFormatExc("PPFBinaryInputStream::GetPosTagTable Size mismatch ");
  for(int kk=0; kk<tsz; kk++)
    GetRawI8(ptab[kk]);
  return;
}

void
PPFBinaryInputStream::GetPosTagTable(vector<int_8>& ptab)
{
  unsigned char ppstype;
  GetTypeTag(ppstype);
  if (ppstype != PPS_POSTAG_TABLE) 
    throw FileFormatExc("PPFBinaryInputStream::GetPosTagTable  bad type in ppersist stream");
  ptab.clear();
  int_4 tsz;
  GetRawI4(tsz);
  int_8 tpos;
  for(int kk=0; kk<tsz; kk++) {
    GetRawI8(tpos);
    ptab.push_back(tpos);
  }
  return;
}


void
PPFBinaryInputStream::AnalyseTags(int lev)
{
  unsigned char ppstag=0;
  unsigned char ppst1,ppst2,ppst3,ppst30;
  int_8 cpos,fsize;
  uint_8 ui8,cid,oid;
  int_4 i4;
  int_8 i8;
  char * buff;
  string str;
  
  cout << "\n ---------------------------------------------------------- " << endl;
  cout << " PPFBinaryInputStream::AnalyseTags(Level= " << lev << ")" << endl; 


  cpos = s->tellg();
  s->seekg(0,ios::end);
  fsize = s->tellg();
  s->seekg(cpos,ios::beg);

  cout << "   Version= " << Version() << " FileSize= " << fsize 
       << " Creation Date= " << CreationDateStr() <<  endl;
  
  uint_8 totntags = 0;
  bool eofok = false;

  while ( (ppstag != PPS_EOF) &&  (cpos < fsize) ) {
    cpos = s->tellg();
    GetRawUByte(ppstag);
    totntags++;

    ppst1 = ppstag&0x0f;  // bits 0123
    ppst2 = ppstag&0x30;  // bits     45
    ppst3 = ppstag&0xc0;  // bits       67
    ppst30 = ppstag&0xc1; // bits 0     67
    if ((ppst2 == 0) && (ppst3 == 0) ) {
      switch (ppst1) {

      case PPS_NULL :
        if (lev > 1)  cout << "<PPS_NULL> tag at position " << hex << cpos << dec << endl;
        break;

      case PPS_STRING :
        GetRawI4(i4);
        if (lev > 1)  cout << "<PPS_STRING> tag at position " << hex << cpos << dec 
                           << " Length=" << i4 << endl;
        s->seekg(i4,ios::cur);
        break;

      case PPS_OBJECT :
        GetRawU8(cid);
        GetRawU8(oid);
        cout << "<PPS_OBJECT> tag at position " << hex << cpos << " ClassId= " << cid 
             << "  ObjectId= " << oid << dec << endl;
        break;

      case PPS_REFERENCE :
        GetRawU8(oid);
        GetRawI8(i8);
        cout << "<PPS_REFERENCE> tag at position " << hex << cpos << "  ObjectId= " 
             << oid << "  OrigPos=" << i8 << dec << endl;
        break;

      case PPS_NAMETAG_MARK :
        cout << "<PPS_NAMETAG_MARK> tag at position " << hex << cpos << dec << endl; 
        break;

      case PPS_POSTAG_MARK :
        GetRawI8(i8);   
        cout << "<PPS_POSTAG_MARK> tag at position " << hex << cpos 
             << " TPos=" << i8 << dec << endl; 
        break;
     
      case PPS_ENDOBJECT :
        GetRawU8(oid);
        cout << "<PPS_ENDOBJECT> tag at position " << hex << cpos << "  ObjectId= " 
             << oid << dec << endl;
        break;

      case PPS_POSTAG_TABLE :
        GetRawI4(i4);
        for(int kkt=0; kkt<i4; kkt++)  GetRawI8(i8);    
        cout << "<PPS_POSTAG_TABLE> tag at position " << hex << cpos << dec << endl; 
        break;

      case PPS_NAMETAG_TABLE :
        GetRawI8(i8);
        GetRawI4(i4);
        buff = new char[i4+1];
        GetRawBytes(buff, i4); 
        buff[i4] = '\0';  str = buff;
        delete[] buff;
        cout << "<PPS_TAG> tag at position " << hex << cpos << dec 
             << " Name= " << str << endl; 
        break;

      case PPS_EOF :
        GetRawI8(i8);   
        cout << "<PPS_EOF> tag at position " << hex << cpos  
             << " TagPos=" << i8 << dec << endl;
        eofok = true;
        break;

      default :
        cerr << " ERROR : Unexpected tag value " << hex << ppstag 
             << " At position" << cpos << dec << endl;
        throw FileFormatExc("PPFBinaryInputStream::AnalyseTags() - Unexpected tag value !");
      }
    }
    else {
      string dtype = "???? x";
      if (ppst30 == PPS_DATATYPE_COMPLEX) dtype = "COMPLEX x";
      else if (ppst3 == PPS_DATATYPE_CHAR) dtype = "CHAR x";
      else if (ppst3 == PPS_DATATYPE_FLOAT) dtype = "FLOAT x";
      else if (ppst3 == PPS_DATATYPE_INTEGER) dtype = "INTEGER x";
      else if (ppst3 == PPS_DATATYPE_UNSIGNED) dtype = "UNSIGNED x";
      int_4 dsize = ppst1;
      int_4 dsizeskip = dsize;
      if (ppst30 == PPS_DATATYPE_COMPLEX) { 
        dsize--;
        dsizeskip = 2*dsize;
      }
      char sb[16];
      sprintf(sb, "%d", dsize);
      dtype += sb;

      switch (ppst2) {

      case PPS_SIMPLE :
        if (lev > 2)  cout << "<PPS_SIMPLE> tag at position " << hex << cpos << dec 
                           << " DataType=" << dtype << endl;
        s->seekg(dsizeskip, ios::cur);
        break;

      case PPS_SIMPLE_ARRAY4 :
        GetRawI4(i4);
        if (lev > 0)  cout << "<PPS_SIMPLE_ARRAY4> tag at position " << hex << cpos << dec 
                           << " DataType=" << dtype << " NElts= " << i4 << endl;
        s->seekg((int_8)dsizeskip*(int_8)i4, ios::cur);
        break;

      case PPS_SIMPLE_ARRAY8 :
        GetRawU8(ui8);
        if (lev > 0)  cout << "<PPS_SIMPLE_ARRAY8> tag at position " << hex << cpos << dec 
                           << " DataType=" << dtype << " NElts= " << ui8 << endl;
        s->seekg((int_8)dsizeskip*ui8, ios::cur);
        break;
      }
    }
  }
  if (!eofok) 
    throw FileFormatExc("PPFBinaryInputStream::AnalyseTags() - Not found <PPS_EOF> tag ");

  cout << " PPFBinaryInputStream::AnalyseTags() - End - Total Number of Tags= " << totntags << endl; 
  cout << " ---------------------------------------------------------- \n" << endl;
  return;
}


//++
// Class        POutPersist
// Lib          Outils++
// include      ppersist.h
//
//      Fichier d'objets persistants, en écriture.
//--


//++
//  POutPersist(string const& flnm, int endianness = PPersist::PPS_NATIVE)
//
//      Crée un nouveau fichier ppersist. Par défaut, il est petit=boutien
//      sur machines petit-boutiennes, et gros-boutien sur machines
//      gros-boutiennes. On peut explicitement spécifier PPersist::PPS_LITTLE_ENDIAN
//      ou PPersist::PPS_BIG_ENDIAN.
//--

PPFBinaryOutputStream::PPFBinaryOutputStream(RawInOutStream* os, bool ad, int endianness)
{
  s = os;
  _ads = ad;
  Init(endianness);
}

PPFBinaryOutputStream::PPFBinaryOutputStream(string const& flnm, int endianness)
{
  // Output stream creation
  s = new RawOutFileStream(flnm.c_str());
  _ads = true;
  Init(endianness);
}

PPFBinaryOutputStream::~PPFBinaryOutputStream()
{
  WriteNameTagTable();
  if (_ads && (s != NULL)) delete s;   // Close the stream
}

void
PPFBinaryOutputStream::Init(int endianness)
{
  if (endianness == -1)
    bigEndian = IS_BIG_ENDIAN;
  else
    bigEndian = endianness;

  version = 3;
  // Header
  PutRawBytes("SOS-SOPHYA-PPersistFile V3               ",32);
  PutRawBytes(bigEndian 
           ? "BIG-ENDIAN                             "
           : "LITTLE-ENDIAN                          ",32);

// ---- GMT creation date of the file
  time_t tm = time(NULL);
  char datestring[33];
  int l=strftime(datestring,32,"%d/%m/%Y %H:%M:%S GMT",gmtime(&tm));
  for(int i=l; i<32; i++)  datestring[i] = ' ';
  datestring[32] = '\0'; 
  PutRawBytes(datestring, 32);
}

void
PPFBinaryOutputStream::WriteNameTagTable()
{  
  if (tags.size() == 0) {
    PutRawUByte(PPS_EOF);
    PutRawI8(-1);
  } 
  else {
    int_8 tagPos;
    tagPos = s->tellp();
    PutRawUByte(PPS_NAMETAG_TABLE);                       // This is a tag
    PutU8((uint_8)tags.size());    // Number of tags
    for (map<string,int_8>::iterator i = tags.begin(); i != tags.end(); i++) {
      int_8 pos = (*i).second;
      PutI8(pos);
      PutStr((*i).first);
    }
    PutRawUByte(PPS_EOF);
    PutRawI8(tagPos);
  }
  
}

void
PPFBinaryOutputStream::WriteNameTagTableV2()
{
  if (tags.size() == 0) {
    PutRawUByte(PPS_EOF);
    PutRawI8(-1);
  } else {
    int_8 tagPos;
    tagPos = s->tellp();
    for (map<string,int_8>::iterator i = tags.begin(); i != tags.end(); i++) {
      string name = (*i).first;
      int_8 pos = (*i).second;
      PutRawUByte(PPS_NAMETAG_TABLE);                       // This is a tag
      PutRawI8(pos);                             // position of previous tag
      PutRawI4(name.length());                   // length of the name
      PutRawBytes(name.c_str(), name.length());  // name, without final "0".
    }
    PutRawUByte(PPS_EOF);
    PutRawI8(tagPos);
  }

}


int_8
PPFBinaryOutputStream::WritePositionTag()
{
  int_8 tagpos;
  tagpos = s->tellp();
  PutRawByte(PPS_POSTAG_MARK);
  PutI8(tagpos);
  return(tagpos);
}

void
PPFBinaryOutputStream::WriteNameTag(string const& name)
{
  //  if (name.length() > MAXTAGLEN_V2) 
  //  throw ParmError("PPFBinaryOutputStream::WriteNameTag   tag name too long");

  if (tags.find(name) != tags.end())
    throw DuplicateIdExc("PPFBinaryOutputStream::WriteNameTag   duplicate tag name");

  // Get current file position
  int_8 tagPos;
  tagPos = s->tellp();

  tags[name] = tagPos;
  PutRawUByte(PPS_NAMETAG_MARK);                       // This is a tag
  //  objList.clear(); // $CHECK$ EA 171199  - Ne pas faire ? Reza 03/2000
}

//++
// void PPFBinaryOutputStream::PutByte(char& c)
// void PPFBinaryOutputStream::PutBytes(void const* ptr, size_t bytes)
// void PPFBinaryOutputStream::PutR4   (r_4 result)
// void PPFBinaryOutputStream::PutR4s  (r_4 const* tab, size_t n)
// void PPFBinaryOutputStream::PutR8   (r_8 result)
// void PPFBinaryOutputStream::PutR8s  (r_8 const* tab, size_t n)
// void PPFBinaryOutputStream::PutI2   (int_2 result)
// void PPFBinaryOutputStream::PutI2s  (int_2 const* tab, size_t n)
// void PPFBinaryOutputStream::PutU2   (uint_2 result)
// void PPFBinaryOutputStream::PutU2s  (uint_2 const* tab, size_t n)
// void PPFBinaryOutputStream::PutI4   (int_4 result)
// void PPFBinaryOutputStream::PutI4s  (int_4 const* tab, size_t n)
// void PPFBinaryOutputStream::PutU4   (uint_4 result)
// void PPFBinaryOutputStream::PutU4s  (uint_4 const* tab, size_t n)
// void PPFBinaryOutputStream::PutI8   (int_8 result)
// void PPFBinaryOutputStream::PutI8s  (int_8 const* tab, size_t n)
// void PPFBinaryOutputStream::PutU8   (uint_8 result)
// void PPFBinaryOutputStream::PutU8s  (uint_8 const* tab, size_t n)
// void PPFBinaryOutputStream::PutStr  (string const&)
//      Ecriture de données portables.. Pour chaque type
//      de données, on peut écrire une valeur, ou un tableau de valeurs.
// void PPFBinaryOutputStream::PutLine(char const* ptr, size_t len)
//      Ecriture d'une ligne de texte dans le fichier PPersist.
//--




void
PPFBinaryOutputStream::PutRawBytes(void const* ptr, size_t bytes)
{
  s->write((char const*)ptr, bytes);
}

void
PPFBinaryOutputStream::PutRawByte(char c)
{
  PutRawBytes(&c, 1);
}

void
PPFBinaryOutputStream::PutRawUByte(unsigned char c)
{
  PutRawBytes(&c, 1);
}

void
PPFBinaryOutputStream::PutRawI2   (int_2 val)
{
  if (bigEndian != IS_BIG_ENDIAN)
    bswap2(&val);

  PutRawBytes(&val, sizeof(int_2));
}

void
PPFBinaryOutputStream::PutRawI4   (int_4 val)
{
  if (bigEndian != IS_BIG_ENDIAN)
    bswap4(&val);

  PutRawBytes(&val, sizeof(int_4));
}

void
PPFBinaryOutputStream::PutRawI8   (int_8 val)
{
  if (bigEndian != IS_BIG_ENDIAN)
    bswap8(&val);

  PutRawBytes(&val, sizeof(int_8));
}

void
PPFBinaryOutputStream::PutRawU8   (uint_8 val)
{
  if (bigEndian != IS_BIG_ENDIAN)
    bswap8(&val);

  PutRawBytes(&val, sizeof(uint_8));
}

void
PPFBinaryOutputStream::PutArrayTag(short datasz, size_t sz, short datatype) 
// datatype = PPS_DATATYPE_CHAR or PPS_DATATYPE_FLOAT or PPS_DATATYPE_INTEGER or PPS_DATATYPE_UNSIGNED
{
  if (sz <= 0x7fffffff) {
    PutRawUByte(PPS_SIMPLE_ARRAY4 + datasz + datatype);
    PutRawI4(sz);
  } else {
    PutRawUByte(PPS_SIMPLE_ARRAY8 + datasz + datatype);
    PutRawU8(sz);
  } 
}

void
PPFBinaryOutputStream::PutByte(char c)
{
  PutRawByte(PPS_SIMPLE + 1 + PPS_DATATYPE_CHAR);
  PutRawBytes(&c, 1);
}



void
PPFBinaryOutputStream::PutBytes(void const* ptr, size_t bytes)
{
  PutArrayTag(1, bytes, PPS_DATATYPE_CHAR);
  PutRawBytes(ptr, bytes);
}

void
PPFBinaryOutputStream::PutR4   (r_4 val)
{
  PutRawUByte(PPS_SIMPLE + 4 + PPS_DATATYPE_FLOAT);
  
  if (bigEndian != IS_BIG_ENDIAN)
    bswap4(&val);

  PutRawBytes(&val, sizeof(r_4));
}

void
PPFBinaryOutputStream::PutR4s  (r_4 const* tab, size_t n)
{
  PutArrayTag(4, n, PPS_DATATYPE_FLOAT);

  if (bigEndian == IS_BIG_ENDIAN) {
    PutRawBytes(tab, n*sizeof(r_4));
  } else {
    for (unsigned int i=0; i<n; i++) {
      r_4 val = tab[i];
      bswap4(&val);
      PutRawBytes(&val, sizeof(r_4));
    }
  }
}

void
PPFBinaryOutputStream::PutR8   (r_8 val)
{
  PutRawUByte(PPS_SIMPLE + 8 + PPS_DATATYPE_FLOAT);

  if (bigEndian != IS_BIG_ENDIAN)
    bswap8(&val);

  PutRawBytes(&val, sizeof(r_8));
}

void
PPFBinaryOutputStream::PutR8s  (r_8 const* tab, size_t n)
{
  PutArrayTag(8, n, PPS_DATATYPE_FLOAT);

  if (bigEndian == IS_BIG_ENDIAN) {
    PutRawBytes(tab, n*sizeof(r_8));
  } else {
    for (unsigned int i=0; i<n; i++) {
      r_8 val = tab[i];
      bswap8(&val);
      PutRawBytes(&val, sizeof(r_8));
    }
  }
}

void
PPFBinaryOutputStream::PutI1   (int_1 val)
{
  PutRawUByte(PPS_SIMPLE + 1 + PPS_DATATYPE_INTEGER);
  PutRawBytes(&val, sizeof(int_1));
}

void
PPFBinaryOutputStream::PutI1s  (int_1 const* tab, size_t n)
{
  PutArrayTag(1, n, PPS_DATATYPE_INTEGER);
  PutRawBytes(tab, n*sizeof(int_1));
}

void
PPFBinaryOutputStream::PutU1   (uint_1 val)
{
  PutRawUByte(PPS_SIMPLE + 1 + PPS_DATATYPE_UNSIGNED);
  PutRawBytes(&val, sizeof(uint_1));
}

void
PPFBinaryOutputStream::PutU1s  (uint_1 const* tab, size_t n)
{
  PutArrayTag(1, n, PPS_DATATYPE_UNSIGNED);
  PutRawBytes(tab, n*sizeof(uint_1));
}

void
PPFBinaryOutputStream::PutI2   (int_2 val)
{
  PutRawUByte(PPS_SIMPLE + 2 + PPS_DATATYPE_INTEGER);

  if (bigEndian != IS_BIG_ENDIAN)
    bswap2(&val);

  PutRawBytes(&val, sizeof(int_2));
}

void
PPFBinaryOutputStream::PutI2s  (int_2 const* tab, size_t n)
{
  PutArrayTag(2, n, PPS_DATATYPE_INTEGER);
  
  if (bigEndian == IS_BIG_ENDIAN) {
    PutRawBytes(tab, n*sizeof(int_2));
  } else {
    for (unsigned int i=0; i<n; i++) {
      int_2 val = tab[i];
      bswap2(&val);
      PutRawBytes(&val, sizeof(int_2));
    }
  }
}

void
PPFBinaryOutputStream::PutU2   (uint_2 val)
{
  PutRawUByte(PPS_SIMPLE + 2 + PPS_DATATYPE_UNSIGNED);

  if (bigEndian != IS_BIG_ENDIAN)
    bswap2(&val);

  PutRawBytes(&val, sizeof(uint_2));
}

void
PPFBinaryOutputStream::PutU2s  (uint_2 const* tab, size_t n)
{
  PutArrayTag(2, n, PPS_DATATYPE_UNSIGNED);

  if (bigEndian == IS_BIG_ENDIAN) {
    PutRawBytes(tab, n*sizeof(uint_2));
  } else {
    for (unsigned int i=0; i<n; i++) {
      uint_2 val = tab[i];
      bswap2(&val);
      PutRawBytes(&val, sizeof(uint_2));
    }
  }
}

void
PPFBinaryOutputStream::PutI4   (int_4 val)
{
  PutRawUByte(PPS_SIMPLE + 4 + PPS_DATATYPE_INTEGER);

  if (bigEndian != IS_BIG_ENDIAN)
    bswap4(&val);

  PutRawBytes(&val, sizeof(int_4));
}

void
PPFBinaryOutputStream::PutI4s  (int_4 const* tab, size_t n)
{
  PutArrayTag(4, n, PPS_DATATYPE_INTEGER);

  if (bigEndian == IS_BIG_ENDIAN) {
    PutRawBytes(tab, n*sizeof(int_4));
  } else {
    for (unsigned int i=0; i<n; i++) {
      int_4 val = tab[i];
      bswap4(&val);
      PutRawBytes(&val, sizeof(int_4));
    }
  }
}

void
PPFBinaryOutputStream::PutU4   (uint_4 val)
{
  PutRawUByte(PPS_SIMPLE + 4 + PPS_DATATYPE_UNSIGNED);

  if (bigEndian != IS_BIG_ENDIAN)
    bswap4(&val);

  PutRawBytes(&val, sizeof(uint_4));
}

void
PPFBinaryOutputStream::PutU4s  (uint_4 const* tab, size_t n)
{
  PutArrayTag(4, n, PPS_DATATYPE_UNSIGNED);

  if (bigEndian == IS_BIG_ENDIAN) {
    PutRawBytes(tab, n*sizeof(uint_4));
  } else {
    for (unsigned int i=0; i<n; i++) {
      uint_4 val = tab[i];
      bswap4(&val);
      PutRawBytes(&val, sizeof(uint_4));
    }
  }
}

void
PPFBinaryOutputStream::PutI8   (int_8 val)
{
  PutRawUByte(PPS_SIMPLE + 8 + PPS_DATATYPE_INTEGER);

  if (bigEndian != IS_BIG_ENDIAN)
    bswap8(&val);

  PutRawBytes(&val, sizeof(int_8));
}

void
PPFBinaryOutputStream::PutI8s  (int_8 const* tab, size_t n)
{
  PutArrayTag(8, n, PPS_DATATYPE_INTEGER);

  if (bigEndian == IS_BIG_ENDIAN) {
    PutRawBytes(tab, n*sizeof(int_8));
  } else {
    for (unsigned int i=0; i<n; i++) {
      int_8 val = tab[i];
      bswap8(&val);
      PutRawBytes(&val, sizeof(int_8));
    }
  }
}

void
PPFBinaryOutputStream::PutU8   (uint_8 val)
{
  PutRawUByte(PPS_SIMPLE + 8 + PPS_DATATYPE_UNSIGNED);

  if (bigEndian != IS_BIG_ENDIAN)
    bswap8(&val);

  PutRawBytes(&val, sizeof(uint_8));
}

void
PPFBinaryOutputStream::PutU8s  (uint_8 const* tab, size_t n)
{
  PutArrayTag(8, n, PPS_DATATYPE_UNSIGNED);

  if (bigEndian == IS_BIG_ENDIAN) {
    PutRawBytes(tab, n*sizeof(uint_8));
  } else {
    for (unsigned int i=0; i<n; i++) {
      uint_8 val = tab[i];
      bswap8(&val);
      PutRawBytes(&val, sizeof(uint_8));
    }
  }
}

void
PPFBinaryOutputStream::PutStr(string const& str) 
{
  PutRawUByte(PPS_STRING);
  PutRawI4(str.length());
  PutRawBytes(str.c_str(), str.length());
}

void
PPFBinaryOutputStream::PutLine(char const* ptr, size_t len)
{
  string str = ptr;
  PutStr(str);
}


void
PPFBinaryOutputStream::PutZ4   (complex<r_4> val)
{
  PutRawUByte(PPS_SIMPLE + 4 + PPS_DATATYPE_COMPLEX);
  r_4 reim[2];
  reim[0] = val.real();
  reim[1] = val.imag();
  if (bigEndian != IS_BIG_ENDIAN) {
    bswap4(reim);
    bswap4(reim+1);
  }
  PutRawBytes(reim, 2*sizeof(r_4));
}

void
PPFBinaryOutputStream::PutZ4s  (complex<r_4> const* tab, size_t n)
{
  PutArrayTag(4, n, PPS_DATATYPE_COMPLEX);

  if (bigEndian == IS_BIG_ENDIAN) {
    PutRawBytes(tab, n*2*sizeof(r_4));
  } else {
    for (unsigned int i=0; i<n; i++) {
      r_4 reim[2];
      reim[0] = tab[i].real();
      reim[1] = tab[i].imag();
      bswap4(reim);
      bswap4(reim+1);
      PutRawBytes(reim, 2*sizeof(r_4));
    }
  }
}

void
PPFBinaryOutputStream::PutZ8   (complex<r_8> val)
{
  PutRawUByte(PPS_SIMPLE + 8 + PPS_DATATYPE_COMPLEX);
  r_8 reim[2];
  reim[0] = val.real();
  reim[1] = val.imag();
  if (bigEndian != IS_BIG_ENDIAN) {
    bswap8(reim);
    bswap8(reim+1);
  }
  PutRawBytes(reim, 2*sizeof(r_8));
}

void
PPFBinaryOutputStream::PutZ8s  (complex<r_8> const* tab, size_t n)
{
  PutArrayTag(8, n, PPS_DATATYPE_COMPLEX);

  if (bigEndian == IS_BIG_ENDIAN) {
    PutRawBytes(tab, n*2*sizeof(r_8));
  } else {
    for (unsigned int i=0; i<n; i++) {
      r_8 reim[2];
      reim[0] = tab[i].real();
      reim[1] = tab[i].imag();
      bswap8(reim);
      bswap8(reim+1);
      PutRawBytes(reim, 2*sizeof(r_8));
    }
  }
}

void
PPFBinaryOutputStream::PutPosTagTable(int_8 const * ptab, size_t sz)
{
  PutRawUByte(PPS_POSTAG_TABLE);
  int_4 tsz = sz;
  PutRawI4(tsz);
  for(int kk=0; kk<tsz; kk++)
    PutRawI8(ptab[kk]);
  return;
}

void
PPFBinaryOutputStream::PutPosTagTable(vector<int_8> const& ptab)
{
  PutRawUByte(PPS_POSTAG_TABLE);
  int_4 tsz = ptab.size();
  PutRawI4(tsz);
  for(int kk=0; kk<tsz; kk++)
    PutRawI8(ptab[kk]);
  return;
}