source: Sophya/trunk/SophyaExt/FitsIOServer/fitsxntuple.cc@ 1047

#include "pexceptions.h"
#include "fitsxntuple.h"
///////////////////////////////////////////////////////////
//   Les objets delegues pour la gestion de persistance sur fichiers fits 
//    pout XNTuple
///////////////////////////////////////////////////////////


#define LONNOM 31 


FITS_XNTuple::FITS_XNTuple()
{
  dobj_ = new XNTuple;
  InitNull();
  ownobj_ = true;
}

FITS_XNTuple::FITS_XNTuple(char inputfile[],int hdunum)
{
  dobj_ = new XNTuple;
  InitNull();
  ownobj_ = true; 

  ReadF(inputfile,hdunum);
}


FITS_XNTuple::FITS_XNTuple(const XNTuple & obj) 
{ 
  dobj_ = new XNTuple(obj);
  InitNull();
  ownobj_ = true; 
}
FITS_XNTuple::FITS_XNTuple(XNTuple* obj) 
{ 
  dobj_ = obj;
  InitNull();
  ownobj_ = false; 
}
FITS_XNTuple::~FITS_XNTuple()
{
  Clean();
}

void FITS_XNTuple::Clean()
{
  if (dcolumn_ != NULL) delete [] dcolumn_;
  if (fcolumn_ != NULL) delete [] fcolumn_;
  if (icolumn_ != NULL) delete [] icolumn_;
  if (ccolumn_ != NULL)
    {
      if (dobj_ != NULL)
        {
          for (int k=0; k<dobj_->NEntry(); k++) delete [] ccolumn_[k];
          delete [] ccolumn_;
        }
      else
        {
          cout << "FITS_XNTuple, destructeur: bizarre, ccolumn non vide, sans objet?" << endl;;         }
    }
  if (ownobj_ && dobj_ != NULL) delete dobj_;
}

void FITS_XNTuple::Read(char inputfile[],int hdunum)
{
  ReadF(inputfile,hdunum);
}

void FITS_XNTuple::ReadLines(char inputfile[],int firstLine, int numberOfLines,int hdunum)
{
  fistLineToBeRead_ = firstLine;
  numberOfLinesToBeRead_ = numberOfLines;
  ReadF(inputfile,hdunum);
  //  return dobj_;
}


void FITS_XNTuple::Write(char outputfile[], bool OldFile)
{       
  WriteF(outputfile, OldFile);
} 
 
void FITS_XNTuple::ReadFromFits()
{   
  //  if (!fn.IsFitsTable())
  if (!IsFitsTable())
    {
      throw PException("ReadFromFits: the fits file seems not to be a bintable nor ASCII table");
    }
  int nbcols, nbentries;
  //  nbcols = fn.NbColsFromFits();
  nbcols = NbColsFromFits();
  nbentries = 0;
  int k;
  //  for (k=0; k<nbcols; k++) nbentries=max( nbentries, fn.NentriesFromFits(k) );
  for (k=0; k<nbcols; k++) nbentries=max( nbentries, NentriesFromFits(k) );

  //
  // pour mettre les colonnes dans l'ordre double, float, int, char :
  // tableau de correspondance 
  // DfitsCol(j)= numero dans le fichier fits de la  jeme variable double du 
  // xntuple;
  // FfitsCol(j)= numero dans le fichier fits de la  jeme variable float du 
  // xntuple;
  // etc.
  vector<int> DfitsCol;
  vector<int> FfitsCol;
  vector<int> IfitsCol;
  vector<int> SfitsCol;
  for (k=0; k<nbcols;k++)
    {
      //      char ss= fn.ColTypeFromFits(k);
      char ss= ColTypeFromFits(k);
      if (ss == 'D') DfitsCol.push_back(k);
      else if (ss == 'E') FfitsCol.push_back(k);
      else if (ss == 'I') IfitsCol.push_back(k);
      else if (ss == 'S') SfitsCol.push_back(k);
      else {
        cout << " FITS_XNTuple: colonne fits " << k << " type= " << ss << endl;
        throw  IOExc("type de champ inconnu");
      }
    }
  char ** ColName = new char*[nbcols];
  int compt=0;
  for (k=0; k<DfitsCol.size(); k++)
    {
      ColName[compt] = new char[LONNOM+1]; 
      //      strncpy(ColName[compt], fn.ColNameFromFits(DfitsCol[k]).c_str(), LONNOM);
      strncpy(ColName[compt], ColNameFromFits(DfitsCol[k]).c_str(), LONNOM);
      ColName[compt++][ LONNOM] =  '\0';
    }
  for (k=0; k<FfitsCol.size(); k++)
    {
      ColName[compt] = new char[LONNOM+1]; 
      //      strncpy(ColName[compt], fn.ColNameFromFits(FfitsCol[k]).c_str(), LONNOM);
      strncpy(ColName[compt], ColNameFromFits(FfitsCol[k]).c_str(), LONNOM);
      ColName[compt++][ LONNOM] =  '\0';
    }
  for (k=0; k<IfitsCol.size(); k++)
    {
      ColName[compt] = new char[LONNOM+1]; 
      //      strncpy(ColName[compt], fn.ColNameFromFits(IfitsCol[k]).c_str(), LONNOM);
      strncpy(ColName[compt], ColNameFromFits(IfitsCol[k]).c_str(), LONNOM);
      ColName[compt++][ LONNOM] =  '\0';
    }
  for (k=0; k<SfitsCol.size(); k++)
    {
      ColName[compt] = new char[LONNOM+1]; 
      //      strncpy(ColName[compt], fn.ColNameFromFits(SfitsCol[k]).c_str(), LONNOM);
      strncpy(ColName[compt], ColNameFromFits(SfitsCol[k]).c_str(), LONNOM);
      ColName[compt++][LONNOM] =  '\0';
    }

  if(dobj_ == NULL) 
    { 
      dobj_= new XNTuple(DfitsCol.size(), FfitsCol.size(), IfitsCol.size(), SfitsCol.size(),ColName);
      ownobj_ = true;      
    }
  else 
    {
      if (ownobj_)
        {
      (*dobj_)= XNTuple(DfitsCol.size(), FfitsCol.size(), IfitsCol.size(), SfitsCol.size(),ColName);
        }
      else
        {
        if (DfitsCol.size() != dobj_->NDVar() || FfitsCol.size() != dobj_->NFVar() ||  IfitsCol.size() != dobj_->NIVar() || SfitsCol.size() != dobj_->NSVar())

          throw  SzMismatchError("FITS_XNTuple : structure incorrecte du ntuple");
        }
    }

  for (k=0; k<nbcols;k++) 
    {
      delete [] ColName[k];
    }
  delete [] ColName;

    r_8* dligne;
    r_4* fligne;
    int_4* iligne;
    char** cligne;

  if (DfitsCol.size()>0) dligne = new r_8[DfitsCol.size()];
  else dligne=NULL;
  if (FfitsCol.size()>0) fligne = new r_4[FfitsCol.size()];
  else fligne=NULL;
  if (IfitsCol.size()>0) iligne = new int_4[IfitsCol.size()];
  else iligne=NULL;
  if (SfitsCol.size()>0) 
    {
      cligne = new char*[SfitsCol.size()];
      int  taille_des_chaines=0;
      //      for (k=0; k< SfitsCol.size(); k++)  taille_des_chaines = max( taille_des_chaines, fn.ColStringLengthFromFits(SfitsCol[k]) );
      for (k=0; k< SfitsCol.size(); k++)  taille_des_chaines = max( taille_des_chaines, ColStringLengthFromFits(SfitsCol[k]) );
      for (k=0; k<SfitsCol.size(); k++) cligne[k]=new char[taille_des_chaines+1];
    }
  else cligne=NULL;
  int firstln, lastln;
  if (numberOfLinesToBeRead_ > 0) 
    {
      firstln = fistLineToBeRead_;
      lastln = firstln + numberOfLinesToBeRead_;
    }
  else
    {
      firstln = 0;
      lastln  = nbentries;
    }
  int numLigne;
  for (numLigne=firstln; numLigne < lastln; numLigne++)
    {
      //      fn.GetBinTabLine(numLigne, dligne, fligne, iligne, cligne );
      GetBinTabLine(numLigne, dligne, fligne, iligne, cligne );
      dobj_->Fill(dligne, fligne, iligne, cligne);    
    }
    delete [] dligne;
    delete [] fligne;
    delete [] iligne;
    for (k=0; k< SfitsCol.size(); k++) delete []  cligne[k];
    delete [] cligne; 
    //  dobj_->Info()=fn.DVListFromFits();
  dobj_->Info()=DVListFromFits();

}
void FITS_XNTuple::WriteToFits()
{
  if(dobj_ == NULL) 
    {
      cout << " WriteToFits:: dobj_= null " << endl;
      return;
    }

  // table will have 'ncols'  columns
  int ncols = dobj_->NVar(); 
  // table will have 'nrows' rows
  int nrows = dobj_->NEntry();
  // get names and values from the join DVList object
  DVList dvl= dobj_->Info();
  // extension name
  char* extname = "XNTuple_Binary_tbl"; 
  dvl.Print();
  char** Noms = new char*[ncols];   
  int k;
  for (k=0; k< ncols; k++)
    {
      Noms[k]= new char[LONNOM+1];
      strncpy(Noms[k],dobj_->NomIndex(k).c_str(),LONNOM+1);
    }
   char* types=new char[ncols+1];
  int compt=0;
  for (k=0; k<dobj_->NDVar();k++)
    {
      types[compt++]='D'; 
    }
  for (k=0; k<dobj_->NFVar();k++)
    {
      types[compt++]='E'; 
    }
  for (k=0; k<dobj_->NIVar();k++)
    {
      types[compt++]='I'; 
    }
  for (k=0; k<dobj_->NSVar();k++)
    {
      types[compt++]='A'; 
    }            
  types[ncols]='\0';
  vector<int> StringSizes(dobj_->NSVar());
  for (k=0; k< StringSizes.size(); k++) StringSizes[k]=dobj_->mStrSz;
  // la librairie fitsio ecrit colonne par colonne 
  //    fn.makeHeaderBntblOnFits(types, Noms, nrows, ncols, dvl, extname,StringSizes); 
    makeHeaderBntblOnFits(types, Noms, nrows, ncols, dvl, extname,StringSizes); 
   for (k=0; k< ncols; k++)
    {
      delete [] Noms[k];
    }
   delete [] Noms;
   delete [] types;
  compt=0;  
  for (k=0; k<dobj_->NDVar();k++)
    {
      putColToFits(compt, dobj_->NEntry(), getColDFromObj(compt));
      compt++;
    }
  for (k=0; k<dobj_->NFVar();k++)
    {
      putColToFits(compt, dobj_->NEntry(), getColFFromObj(compt));
      compt++;
    }
  for (k=0; k<dobj_->NIVar();k++)
    {
      putColToFits(compt, dobj_->NEntry(), getColIFromObj(compt));
      compt++;
    }
  for (k=0; k<dobj_->NSVar();k++)
    {
      putColToFits(compt, dobj_->NEntry(), getColSFromObj(compt));
      compt++;
    }
}
double* FITS_XNTuple::getColDFromObj(int colNr)
{
  if (dcolumn_ != NULL)
    {
      delete [] dcolumn_;
      dcolumn_ = NULL;
    }
  dcolumn_ = new double[dobj_->NEntry()];
  for(int j = 0; j < dobj_->NEntry(); j++) dcolumn_[j]= dobj_->GetDVal(j,colNr);
  return dcolumn_;
}

float* FITS_XNTuple::getColFFromObj(int colNr)
{
  if (fcolumn_ != NULL)
    {
      delete [] fcolumn_;
      fcolumn_ = NULL;
    }
  fcolumn_ = new float[dobj_->NEntry()];
  for(int j = 0; j < dobj_->NEntry(); j++) fcolumn_[j]= dobj_->GetFVal(j,colNr);
  return fcolumn_;
}
 
int* FITS_XNTuple::getColIFromObj(int colNr)
{
  if (icolumn_ != NULL)
    {
      delete [] icolumn_;
      icolumn_ = NULL;
    }
  icolumn_ = new int[dobj_->NEntry()];
  for(int j = 0; j < dobj_->NEntry(); j++) icolumn_[j]= dobj_->GetIVal(j,colNr);
  return icolumn_;
}
char** FITS_XNTuple::getColSFromObj(int colNr)
{
  if (ccolumn_ != NULL)
    { 
      for (int k=0; k<dobj_->NEntry(); k++) delete [] ccolumn_[k];
      delete [] ccolumn_;
      ccolumn_ = NULL;
    }
  ccolumn_ = new char*[dobj_->NEntry()];
  for(int j = 0; j < dobj_->NEntry(); j++)
    {
      string s= dobj_->GetSVal(j,colNr);
      ccolumn_[j] = new char[dobj_->mStrSz+1];
      strcpy(ccolumn_[j],s.c_str());
    }
  return ccolumn_; 
}