Changeset 2907 in Sophya for trunk/SophyaExt/FitsIOServer/fitsntuple.cc

Timestamp:

Jan 17, 2006, 10:09:29 AM (20 years ago)

Author:

ansari

Message:

debug/check fits : gestion NTuple colonnes double + init hdunum_ ds FitsFile(FitsInOutFile&) + correction gestion TArray/TMatrix/TVector ds fitsarrhand.h - Reza 17/01/2006

File:

• trunk/SophyaExt/FitsIOServer/fitsntuple.cc (modified) (5 diffs)

trunk/SophyaExt/FitsIOServer/fitsntuple.cc

@@ -104 +104 @@
       ColName[k][LONNOM] =  '\0';
     }
+  // Reza Dec 2006 : Colonnes du NTuple de type double par defaut 
   for (k=0; k<nbcols;k++)
     {
       FitsFile::FitsDataType ss= is.ColTypeFromFits(k);
       string type;
-      if (ss !=  FitsFile::FitsDataType_float) 
+      if ( (ss !=  FitsFile::FitsDataType_double) && (ss != FitsFile::FitsDataType_float) )
         {
-          if (ss == FitsFile::FitsDataType_double) type= string("double");
+          if (ss == FitsFile::FitsDataType_int) type= string("integer");
           else
-            if (ss == FitsFile::FitsDataType_int) type= string("integer");
+            if (ss == FitsFile::FitsDataType_char) type = string("char*");
             else
-              if (ss == FitsFile::FitsDataType_char) type = string("char*");
-              else
-                type = string("unknown");
+              type = string("unknown");
           cout << " WARNING: the column " << k << " on fits file is not float but : " << type << endl;
         } 
     }
+  bool fgdouble = true;
+  if ( is.ColTypeFromFits(0) == FitsFile::FitsDataType_float ) fgdouble = false;
   if(dobj_ == NULL) 
     { 
-      dobj_= new NTuple(nbcols,ColName);
+      dobj_= new NTuple(nbcols,ColName,512,fgdouble);
       ownobj_ = true;      
     }
@@ -128 +129 @@
     {
       dobj_->Clean();
-      (*dobj_) = NTuple(nbcols,ColName);
+      (*dobj_) = NTuple(nbcols,ColName,512,fgdouble);
     }
   for (k=0; k<nbcols;k++) 
@@ -136 +137 @@
   delete [] ColName;
 
-  float* ligne = new float[nbcols];
+  // Reza Dec 2006 : Colonnes du NTuple de type double par defaut 
+  double* dligne = NULL;
+  float* fligne = NULL;
+  if (fgdouble) dligne = new double[nbcols];
+  else fligne = new float[nbcols];
 
   int firstln, lastln;
@@ -150 +155 @@
     }
   int numLigne;
-  for (numLigne=firstln; numLigne < lastln; numLigne++)
-    {
-      is.GetBinTabLine(numLigne, ligne);
-      dobj_->Fill((r_4*)ligne);
-    }
-  delete [] ligne;
+  if (fgdouble) 
+    for (numLigne=firstln; numLigne < lastln; numLigne++) {
+      is.GetBinTabLine(numLigne, dligne);
+      dobj_->Fill((r_8*)dligne);
+    }
+  else 
+    for (numLigne=firstln; numLigne < lastln; numLigne++) {
+      is.GetBinTabLine(numLigne, fligne);
+      dobj_->Fill((r_4*)fligne);
+    }
+  
+  if (dligne) delete [] dligne;
+  if (fligne) delete [] fligne;
 
   dobj_->Info()=is.DVListFromFits();
@@ -191 +203 @@
     }
 
-  string type(ncols, 'E');
-  vector<int> dummy;
-
-  os.makeHeaderBntblOnFits(type,Noms, nentries, ncols, &dvl, extname, dummy);
-   float* column =  new float[nentries];
-   for (k=0; k<ncols;k++)
-     {
-       int j;
-       for(j = 0; j < nentries; j++) column[j]= dobj_->GetVal(j,k);
-       os.PutColToFits(k, nentries, column);
-     }
-   delete [] column;
-
-}
-
+  // Gestion NTuple en double ou float / Reza Jan 2006
+  if (dobj_->mFgDouble) {    // NTuple avec colonnes double
+    string type(ncols, 'D');
+    vector<int> dummy;
+    
+    os.makeHeaderBntblOnFits(type,Noms, nentries, ncols, &dvl, extname, dummy);
+    double* column =  new double[nentries];
+    for (k=0; k<ncols;k++)
+      {
+        int j;
+        for(j = 0; j < nentries; j++) column[j]= dobj_->GetCell(j,k);
+        os.PutColToFits(k, nentries, column);
+      }
+    delete [] column;
+  }
+  else {   // NTuple avec colonnes float
+    string type(ncols, 'E');
+    vector<int> dummy;
+    
+    os.makeHeaderBntblOnFits(type,Noms, nentries, ncols, &dvl, extname, dummy);
+    float* column =  new float[nentries];
+    for (k=0; k<ncols;k++)
+      {
+        int j;
+        for(j = 0; j < nentries; j++) column[j]= dobj_->GetVal(j,k);
+        os.PutColToFits(k, nentries, column);
+      }
+    delete [] column;
+  }
+}
+