Changeset 606 in Sophya for trunk/SophyaProg/Tests/tspm2.cc

Timestamp:

Nov 20, 1999, 9:47:42 PM (26 years ago)

Author:

ansari

Message:

Amelioration test spheres - Reza 20/11/99

File:

• trunk/SophyaProg/Tests/tspm2.cc (modified) (5 diffs)

trunk/SophyaProg/Tests/tspm2.cc

@@ -58 +58 @@
   printf("Xsize= %d   Ysize= %d  NPix= %d\n",img.XSize(),img.YSize(),img.XSize()*img.YSize() );
   n = 0;
+  img.Zero();
   for(j=0; j<img.YSize(); j++) {
     yd = (r_4)(j+0.5)/(r_4)img.YSize()-0.5;
@@ -112 +113 @@
    nomobj = "sphg1";
    fiog.Write(s, nomobj);
-   cout << "SphMap written to POutPersist with name "  << nomobj << endl;
-   }
-
+   cout << "SphMap SphereGorski<double> written to POutPersist with name "  << nomobj << endl;
+   }
+
+  FitsIoServer fios;
   // On projete dans un fichier FITS
-  FitsIoServer fios;
   fios.sinus_picture_projection(sph, "gsin1.fits");
   fios.save(sph, "sph1.fits");
@@ -146 +147 @@
   PrtTim("End of Mean-Sig ");
   }
+
   {
-  SphereGorski<double> sph(m);
+  SphereGorski<float> sph(m);
 
   cout << "Filling spherical map2 NPixels= " << sph.NbPixels() << endl;
@@ -164 +166 @@
 
    {
-   FIO_SphereGorski<double> fiog(&sph) ;
+   FIO_SphereGorski<float> fiog(&sph) ;
    nomobj = "sphg2";
    fiog.Write(s, nomobj);
-   cout << "SphMap written to POutPersist with name "  << nomobj << endl;
+   cout << "SphMap SphereGorski<float> written to POutPersist with name "  << nomobj << endl;
    }
 
@@ -173 +175 @@
   // On projete dans un fichier FITS
   FitsIoServer fios;
+   {
+   cout << "Test of Write/Read SphereGorski<float> to FITS (sphg_r4.fits) " << endl;
+   fios.save(sph, "sphg_r4.fits");
+ 
+   SphereGorski<float> sphr(4);
+   fios.load(sphr, "sphg_r4.fits", 2);
+   cout << " Read from file - SphereGorski<float> NPixels= " << sphr.NbPixels() << endl;
+   int ndiff = 0;
+   for(int k=0; k<sphr.NbPixels(); k++) {
+    if ( sphr(k) != sph(k) )   { 
+      ndiff++;
+      if (ndiff < 20)  cout << "!!!Diff: K= " << k << " SPHR= " << sphr(k) << " SPH= " << sph(k) << endl;
+      }
+    }
+   cout << " ReadFrom FITS NDiff = " << ndiff << " (should be zero = 0) " << endl;
+   
+   }
+
+
   fios.sinus_picture_projection(sph, "gsin2.fits");
   fios.save(sph, "sph2.fits");
 
-  TMatrix<double> mtx(3*m, 6*m);
+  TMatrix<float> mtx(3*m, 6*m);
   Project_Mol(sph, mtx);
   fios.save(mtx, "mtx2.fits");