Changeset 3496 in Sophya for trunk

Timestamp:

Apr 30, 2008, 3:29:47 PM (17 years ago)

Author:

ansari

Message:

MAJ prog test fits suite ajout methodes lecture-ecriture tableaux par morceaux - Reza 30/04/2008

File:

• trunk/SophyaProg/Tests/tnfits.cc (modified) (1 diff)

trunk/SophyaProg/Tests/tnfits.cc

@@ -168 +168 @@
     }
 
+    {
+      cout << "5/=== Ecriture tableau par morceaux dans tamor.fits  " << endl;
+      FitsInOutFile fiosc("!tamor.fits", FitsInOutFile::Fits_Create);
+      cout << "5.a/ TArray<int_4> iam(6,4) = RegularSequence(0., 1.) " << endl;
+      TArray<int_4> iam(6,4);      
+      iam = RegularSequence(0., 1.);
+      LONGLONG naxes[3] = {6,4,3};
+      int_4 x;
+      fiosc.CreateImageHDU(FitsTypes::ImageType(x), 3, naxes);
+      iam.Info()["Info"] = "Write test of 3 x TArray<r_4>(6,4)";
+      fiosc.WriteHeaderRecords(iam.Info());
+      sa_size_t off[3] = {0,0,0};
+      FitsArrayHandler<int_4> fah(iam);
+      for(int_4 kk=0; kk<3; kk++) {
+        iam += 100;
+        off[2] = kk;
+        cout << "5.b/ Appel fah.WriteAtOffset(fiosc, off) Avec off[2]=" 
+             << off[2] << endl;
+        fah.WriteAtOffset(fiosc, off);
+      }
+    }
+
+    {
+      cout << "6/=== lecture tableau par morceaux depuis tamor.fits  " << endl;
+      FitsInOutFile fiosc("tamor.fits", FitsInOutFile::Fits_RO);
+      TArray<int_4> iam(6,4), iamck(6,4), diff(6,4);      
+      iamck = RegularSequence(0., 1.);
+      int rc = 0;
+      int_4 min,max;
+      sa_size_t off[3] = {0,0,0};
+      int rcf[3] = {7,70,700};
+      FitsArrayHandler<int_4> fah(iam);
+      for(int_4 kk=0; kk<3; kk++) {
+        iamck += 100;
+        off[2] = kk;
+        cout << "6.b/ Appel fah.ReadAtOffset(fiosc, off) Avec off[2]=" 
+             << off[2] << endl;
+        fah.ReadAtOffset(fiosc, off);
+        cout << iam;
+        diff = iam-iamck;
+        diff.MinMax(min, max);
+        if ((min != 0) || (max != 0))  rc += rcf[kk];
+      }
+      cout << "6.b/ fah.ReadAtOffset(tamor.fits ...) Rc= " << rc << endl;
+    }
+
   }
   catch (PThrowable & exc) {