Changeset 931 in Sophya for trunk/SophyaProg/PMixer/tgsky.cc

Timestamp:

Apr 14, 2000, 3:37:06 PM (25 years ago)

Author:

ansari

Message:

Sophie: adapt to new class MapOperations

File:

• trunk/SophyaProg/PMixer/tgsky.cc (modified) (3 diffs)

trunk/SophyaProg/PMixer/tgsky.cc

@@ -1 @@
-#include "machdefs.h"
-#include <iostream.h>
-#include <math.h>
+#include "pmixer.h"
+#include "mapoperation.h"
 
-#include "sambainit.h"
-#include "spherehealpix.h"
-#include "fiospherehealpix.h"
-#include "fitsioserver.h"
-#include "nbrandom.h"
-#include "bruit.h"
-#include "timing.h"
 /*! \ingroup PMixer
-  \file tgsky.cc
+ * \file tgsky.cc
  * \brief \b PROGRAM \b tgsky  <BR>
  * \l Program which generates different types of random skies
  */
+// ------------- Main program --------------
 
-// Test program to generate different type of random skies
-template <class T>
-void  MeanSig(PixelMap<T> const & map, double& gmoy, double& gsig);
-
-// ------------- Main program --------------
 int main(int narg, char* arg[]) 
 {
@@ -102 +90 @@
 
   // Computing mean and sigma on the sphere 
-  MeanSig(sph, gmoy, gsig);
+  MeanSig(sph.DataBlock(), gmoy, gsig);
   cout << "SphereHEALPix<float> Mean= " << gmoy << "  Sigma = " << gsig << endl;
   PrtTim("End of Mean-Sig ");
@@ -126 +114 @@
 }
 
-/* Nouvelle-Methode */
-template <class T>
-void  MeanSig(PixelMap<T> const & map, double& gmoy, double& gsig)
-
-{
-  gmoy=0.;
-  gsig = 0.; 
-  double valok;
-  for(int k=0; k<map.NbPixels(); k++) {
-    valok = map(k);
-    gmoy += valok;  gsig += valok*valok; 
-  }
-  gmoy /= (double)map.NbPixels();
-  gsig = gsig/(double)map.NbPixels() - gmoy*gmoy;
-  if (gsig >= 0.) gsig = sqrt(gsig);
-}