Changeset 3089 in Sophya for trunk/SophyaProg/Tests/tswfits.cc

Timestamp:

Oct 11, 2006, 12:58:10 PM (19 years ago)

Author:

ansari

Message:

prog tests: ajout test perfo multi-fft (tfft.cc) et amelioration test SwFits (tswfits.cc), Reza 11/10/2006

File:

• trunk/SophyaProg/Tests/tswfits.cc (modified) (6 diffs)

trunk/SophyaProg/Tests/tswfits.cc

@@ -11 +11 @@
 #include "swfitsdtable.h"
 
+#include "datatable.h"
+#include "fitshdtable.h"
+
 void test_fdtable() ; 
+static long NROW = 64;
+static long NPRT = 10;
 
 int main(int narg, char *arg[])
@@ -17 +22 @@
   SophyaInit();
   FitsIOServerInit();
+  if (narg < 3) {
+    cout << "Erreur arg/Usage: tswfits NRowTable NPRT" << endl;
+    return 1;
+  }
+  NROW = atoi(arg[1]);
+  NPRT = atoi(arg[2]);
   try {
-    /*
-    if (*arg[1] == 'd')  test_dvl();
-    else if (*arg[1] == 'n')  test_ntuple();
-    //    else if (*arg[1] == 'x')  test_xntuple();
-    else if (*arg[1] == 'X')  test_Xntuple();
-    */
     test_fdtable();
   }
@@ -35 +40 @@
     cerr << "tswfits-main() , Catched ... ! " << endl;
   }
+  
 }
 
@@ -40 +46 @@
 void test_fdtable() 
 {
-  cout << "======= test_fdtable:  simple DataTable+FITS test ======= " << endl;
+  cout << "======= test_fdtable:  simple SwapOnFits test ======= " << endl;
     int SWSEGSZ = 4;
     vector<int_8> tags;
     {
-      cout << " Creation FITSDataSwapper<double> swout(so) avec swseg.fits " << endl;
+      cout << "1/ Creation FITSDataSwapper<double> swout(so) avec swseg.fits " << endl;
       vector<string> ttype, tform, tunit;
       ttype.push_back("D-Col1");
@@ -74 +80 @@
     }
     {
-      cout << "6/ Creation PPFDataSwapper<float> avec lecture depuis swseg.ppf " << endl;
+      cout << "2/ Creation FITSDataSwapper<double> avec lecture depuis swseg.fits" << endl;
       FitsInOutFile si("swseg.fits", FitsInOutFile::Fits_RO);
       si.MoveAbsToHDU(2);
       FITSDataSwapper<double> swin(si, 1);
       SwSegDataBlock<double> swsd(swin, tags, SWSEGSZ);
-      cout <<"7/ Recopie ds SegDataBlock<float> depuis SwSegDataBlock<float> et print" << endl;
+      cout <<"2.b/ Recopie ds SegDataBlock<double> depuis SwSegDataBlock<double> et print" << endl;
       SegDataBlock<double> segd;
       segd = swsd;
@@ -86 +92 @@
 
     {
-      cout << "8/ Creation SwFitsDataTable sur myswtable.fits " << endl;
+      cout << " ------------------------------------------------ " << endl;
+      cout << "3/ Creation SwFitsDataTable sur myswtable.fits " << endl;
       FitsInOutFile so("!myswtable.fits", FitsInOutFile::Fits_Create);
-      SwFitsDataTable dt(so, 16);
+      SwFitsDataTable dt(so, 45);
       // define table columns
       dt.AddFloatColumn("X0_f");
       dt.AddFloatColumn("X1_f");
       dt.AddDoubleColumn("X0X0pX1X1_d");
-      cout << "8.b/ start fill " << endl;
+      cout << "3.b/ start fill " << endl;
       // Fill the table
       r_8 x[5];
-      for(int i=0; i<63; i++) {
-        if (i < 5)  cout << " DBG - i= " << i << endl;
+      int pmod = NROW/50;
+      if (pmod < 10) pmod = 10;
+      for(int i=0; i<NROW; i++) {
+        //      if (i%pmod == 0)  cout << " DBG - i= " << i << endl;
         x[0] = (i%9)-4.;  x[1] = (i/9)-3.;  x[2] = x[0]*x[0]+x[1]*x[1];
         dt.AddLine(x);
       }
-      cout << "8.c/ End fill " << endl;
+      cout << "3.c/ End fill " << endl;
       
       // Printing table info
       cout << dt ;
+      for(int k = 0; k<NPRT; k++) 
+        cout << "Line[" << k << "] " << dt.LineToString(k) << endl ; 
+     
     }
+    {
+      cout << " ------------------------------------------------ " << endl;
+      cout << "4/ Lecture SwFitsDataTable depuis myswtable.fits " << endl;
+      FitsInOutFile si("myswtable.fits", FitsInOutFile::Fits_RO);
+      si.SkipEmptyFirstHDU();
+      SwFitsDataTable dt(si,40,false);
+      //      si >> dt;
+      // Printing table info
+      cout << dt ;
+      for(int k = 0; k<NPRT; k++) 
+        cout << "Line[" << k << "] " << dt.LineToString(k) << endl ; 
+    }
+
   cout << "============ FIN  test_fdtable ======== ======= " << endl;  
 }