Changeset 880 in Sophya for trunk/SophyaProg/PMixer/skymixer.cc

Timestamp:

Apr 11, 2000, 3:34:55 PM (25 years ago)

Author:

ansari

Message:

Adaptation SphereGorski -> SphereHEALPix apres modifs Sophie - Reza 11/4/2000

File:

• trunk/SophyaProg/PMixer/skymixer.cc (modified) (16 diffs)

trunk/SophyaProg/PMixer/skymixer.cc

@@ -15 +15 @@
 #include "fitsioserver.h"
 
-#include "spheregorski.h"
+#include "spherehealpix.h"
+#include "fiospherehealpix.h"
 
 #include "radspecvector.h"
@@ -51 +52 @@
 //
 template <class T>
-void addComponentBeta(SphereGorski<T>& finalMap, 
-                      SphereGorski<T>& mapToAdd,SpectralResponse& sr, 
-                      SphereGorski<T>& betaMap, double normFreq, double K);
+void addComponentBeta(SphereHEALPix<T>& finalMap, 
+                      SphereHEALPix<T>& mapToAdd,SpectralResponse& sr, 
+                      SphereHEALPix<T>& betaMap, double normFreq, double K);
 //
 template <class T>
 void integratedMap(SpectralResponse&  sr,   
-                   SphereGorski<T>& betaMap, double normFreq, SphereGorski<T>& intBetaMap);
+                   SphereHEALPix<T>& betaMap, double normFreq, SphereHEALPix<T>& intBetaMap);
   //template <class T> 
   //void integratedMap(SpectralResponse&  sr,   
@@ -64 +65 @@
   //                          PixelMap<T>& outMap);
 template <class T>
-void addComponentBeta(SphereGorski<T>& finalMap, 
-                      SphereGorski<T>& mapToAdd,
-                      SphereGorski<T>& intBetaMap, double K);
+void addComponentBeta(SphereHEALPix<T>& finalMap, 
+                      SphereHEALPix<T>& mapToAdd,
+                      SphereHEALPix<T>& intBetaMap, double K);
 //
 // -----------------------------------------------------------------
@@ -119 +120 @@
 
 
-cout << " skymix/Info : NComp = " <<  nskycomp << " SphereGorski_NSide= " << hp_nside << endl;
+cout << " skymix/Info : NComp = " <<  nskycomp << " SphereHEALPix_NSide= " << hp_nside << endl;
 cout << "  ... MapPath = " << (string)mapPath << "  DebugLev= " << debuglev 
      << "  PrintLev= " << printlev << endl;
@@ -126 +127 @@
 if (debuglev > 0) so = new POutPersist("skymixdbg.ppf");
 
-SphereGorski<float> outgs(hp_nside);
+SphereHEALPix<float> outgs(hp_nside);
 bool okout = false;
 
@@ -156 +157 @@
   string key;
 
-  SphereGorski<float> ings(hp_nside);  // The input map
+  SphereHEALPix<float> ings(hp_nside);  // The input map
   double K = 1.;
   double freqOfMap = 1.;
@@ -182 +183 @@
     fios.load(ings, flnm, 2);
     if (debuglev > 4) {  // Writing the input map to the outppf
-      FIO_SphereGorski<float> fiog(ings);
+      FIO_SphereHEALPix<float> fiog(ings);
       fiog.Write(*so, key);
     }
@@ -199 +200 @@
       {
         key = (string)"BETAFITSFILE"+ buff;
-        SphereGorski<float> betaMap;
+        SphereHEALPix<float> betaMap;
         flnm = BuildFITSFileName(dc.SParam(key, 0));
         double normFreq = dc.DParam(key, 1, 1.);
@@ -215 +216 @@
           {
             // integrate the betamap over the SpectralResponse
-            SphereGorski<float> intBetaMap(nSideForInt);
+            SphereHEALPix<float> intBetaMap(nSideForInt);
             integratedMap(*sr, betaMap, normFreq, intBetaMap);
             betaMap.Resize(8);
@@ -248 +249 @@
   FitsIoServer fios;
   fios.save(outgs, arg[2]);
-  cout << "Output Map (SphereGorski<float>) written to FITS file " 
+  cout << "Output Map (SphereHEALPix<float>) written to FITS file " 
        << (string)(arg[2]) << endl;
   PrtTim("End of WriteFITS ");
@@ -254 +255 @@
   if (narg > 3) {
    POutPersist s(arg[3]); 
-   FIO_SphereGorski<float> fiog(&outgs) ;
+   FIO_SphereHEALPix<float> fiog(&outgs) ;
    fiog.Write(s);
-   cout << "Output Map (SphereGorski<float>) written to POutPersist file "  
+   cout << "Output Map (SphereHEALPix<float>) written to POutPersist file "  
         << (string)(arg[3]) << endl;
    PrtTim("End of WritePPF ");
@@ -503 +504 @@
 /* Nouvelle-Fonction */
 template <class T> 
-void addComponentBeta(SphereGorski<T>& finalMap, 
-                      SphereGorski<T>& mapToAdd,SpectralResponse& sr, 
-                      SphereGorski<T>& betaMap, double normFreq, double K)
+void addComponentBeta(SphereHEALPix<T>& finalMap, 
+                      SphereHEALPix<T>& mapToAdd,SpectralResponse& sr, 
+                      SphereHEALPix<T>& betaMap, double normFreq, double K)
 {
   // finalMap = finalMap + coeff* mapToAdd
@@ -520 +521 @@
       cout << "nb pixels: " << finalMap.NbPixels() << endl;
     }
-  SphereGorski<T> bigBetaMap(sqrt(nbpix/12));
+  SphereHEALPix<T> bigBetaMap(sqrt(nbpix/12));
   if(nbpix != betaMap.NbPixels())
     {
@@ -537 +538 @@
 template <class T>
 void integratedMap(SpectralResponse&  sr,   
-                   SphereGorski<T>& betaMap, 
+                   SphereHEALPix<T>& betaMap, 
                    double normFreq, 
-                   SphereGorski<T>& intBetaMap)
+                   SphereHEALPix<T>& intBetaMap)
 {
   //  int nbpix = intBetaMap.NbPixels();
-  //  SphereGorski<T> newMap(sqrt(nbpix/12));
+  //  SphereHEALPix<T> newMap(sqrt(nbpix/12));
   Sph2Sph(betaMap,intBetaMap);
   for(int i=0; i<intBetaMap.NbPixels(); i++)
@@ -554 +555 @@
 
 template <class T>
-void addComponentBeta(SphereGorski<T>& finalMap, 
-                      SphereGorski<T>& mapToAdd,SphereGorski<T>& intBetaMap, double K)
+void addComponentBeta(SphereHEALPix<T>& finalMap, 
+                      SphereHEALPix<T>& mapToAdd,SphereHEALPix<T>& intBetaMap, double K)
 {
   // finalMap = finalMap + coeff* mapToAdd
@@ -575 +576 @@
     {
       cout << "nbpix != intBetaMap.NbPixels()" << endl;
-      SphereGorski<T> bigBetaMap(sqrt(nbpix/12));
+      SphereHEALPix<T> bigBetaMap(sqrt(nbpix/12));
       cout << "new map with size corresponding to mapToAdd" <<  sqrt(nbpix/12) << endl;
       Sph2Sph(intBetaMap,bigBetaMap);