Changeset 1994 in Sophya for trunk/ArchTOIPipe/Kernel/fitstoirdr.cc

Timestamp:

May 13, 2002, 3:11:32 PM (23 years ago)

Author:

ansari

Message:

1/ Ajout classe CGT (Compagnie Generale des Tuyaux) - Assistance a l'assemblage
des TOIProcessors .
2/ Corrections mineures ds FitsTOIReader, correction de smkmflib

Reza 13/5/2002

File:

• trunk/ArchTOIPipe/Kernel/fitstoirdr.cc (modified) (9 diffs)

trunk/ArchTOIPipe/Kernel/fitstoirdr.cc

@@ -3 +3 @@
 //                               Christophe Magneville
 //                               Reza Ansari
-// $Id: fitstoirdr.cc,v 1.26 2001-12-21 17:00:40 aubourg Exp $
+// $Id: fitstoirdr.cc,v 1.27 2002-05-13 13:11:32 ansari Exp $
 
 #include "fitstoirdr.h"
@@ -27 +27 @@
 }
 
-void FITSTOIReader::setImplicitSN(int snStart) {
+void FITSTOIReader::setImplicitSN(long snStart) {
   implicitSN = true;
   implicitSNStart = snStart;
+}
+
+void FITSTOIReader::setBufferSize(int buffsz)
+{
+  Buff_Sz = (buffsz>0) ? buffsz: 1024;
+  return;
 }
 
@@ -144 +150 @@
   for (vector<string>::iterator i=allfn.begin(); i!=allfn.end(); i++) {
     openFile(*i);
-    //run1();
-    run2();
+    if (Buff_Sz > 1) run2();  // Lecture bufferise
+    else run1();              // Lecture un echantillon a la fois
   }
 }
 
 // run 1 : deprecated. NON MAINTENU. Incompatible avec implicit SN.
+//                     ^^^^^^^^^^^^  Reza , 13/5/2002 , Je viens de rajouter 
+//                                   implicitSNStart et switvh run1/run2 
+//                                   suivant buffersize
+  
 void FITSTOIReader::run1() {
   // Il faudrait optimiser en fonction de ce qui a ete demande comme samplenum,
@@ -165 +175 @@
     double y;
     fits_lock();
-    fits_read_col_dbl(fptr,1,i+1,1,1,0,&y,&anyNul,&fstatus);
+    long sn;
+    if (implicitSN) {
+      fits_read_col_dbl(fptr,1,i+1,1,1,0,&y,&anyNul,&fstatus);
+      sn = (long) (y+.1); 
+    }
+    else sn = implicitSNStart+i;
     //fits_unlock();
-    long sn = (long) (y+.1); 
     TOIManager* mgr = TOIManager::getManager();
     if (sn > mgr->getRequestedEnd()) {fits_unlock(); break;}
@@ -228 +242 @@
  long ideb=-1,ifin=-1;
  for(long i=0; i<nrows; i+=Buff_Sz) {
-
+   fits_lock();   // lock en debut de boucle de lecture - Reza 13/05/2002
    // faut-il lire dans le fichier fits ?
    if(i<ideb || i>ifin) { // Toujours vrai avec le += Buff_Sz
@@ -244 +258 @@
          samplenum[j] = j+ideb+implicitSNStart;
        }
-       fits_lock();
+       //---- Pas la peine, fait en debut de boucle - Reza 13/05/2002  fits_lock();
      } else {
-       fits_lock();
+       //---- Pas la peine, fait en debut de boucle - Reza 13/05/2002 fits_lock();
        fits_read_col_dbl(fptr,1,ideb+1,1,n,0,samplenum,&anyNul,&fstatus);
        if (samplenum[0] > mgr->getRequestedEnd() 
            || (forcedMaxIn > 0 && samplenum[0] > forcedMaxIn)) {
-         fits_unlock();
+         fits_unlock();  // unlock avant de casser la boucle - Reza 13/05/2002
          break;
        }
        if (samplenum[n-1] < mgr->getRequestedBegin()
            || (forcedMinIn > 0 && samplenum[n-1] < forcedMinIn)) {
-         fits_unlock();
+         //---- Pas la peine, fait en fin de boucle - Reza 13/05/2002 fits_unlock();
          continue;
        }
@@ -268 +282 @@
      if(fstatus!=0) {
        fits_report_error(stderr,fstatus);
+       fits_unlock(); // On supprime le lock avant le throw - Reza 13/05/2002
        throw RangeCheckError("FITSTOIReader::run2: Error Reading Fits file\n");
      }
-     fits_unlock();
    }
+   fits_unlock();   // unlock en debut de boucle de lecture - Reza 13/05/2002
 
    long ip = i-ideb; // pointeurs dans les buffers
@@ -288 +303 @@
      }
    }
-   totnscount++;
+   totnscount += (ifin-ideb+1);
  }
 
@@ -300 +315 @@
  colval.resize(0); colflg.resize(0);
 
- cout << "reader done reading... " << pthread_self() << endl;
-}
+ cout << "reader (buffered) done reading... " << pthread_self() << endl;
+}