Changeset 937 in Sophya for trunk/SophyaProg/PMixer

Timestamp:

Apr 14, 2000, 6:11:47 PM (25 years ago)

Author:

ansari

Message:

Sophie: changes to read the SphereHEALPIX with a limited amount
of needed memory

File:

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

trunk/SophyaProg/PMixer/skymixer.cc

@@ -125 +125 @@
       cout << " Reading output HealPix map from FITS file " << (string)ifnm << endl;
       {
-        FITS_SphereHEALPix<float> fios(ifnm);
-        outgs = (SphereHEALPix<float>)fios;
+        FITS_SphereHEALPix<float> fios(&outgs);
+        fios.ReadF(ifnm,2);
       }
       if(printlev>0)
@@ -188 +188 @@
           SphereHEALPix<float> ings(hp_nside);
           {
-            FITS_SphereHEALPix<float> fiosIn(flnm);
-            ings = (SphereHEALPix<float>)fiosIn;
+            FITS_SphereHEALPix<float> fiosIn(&ings);
+            fiosIn.ReadF(flnm,2);
           }
           if (debuglev > 4) {  // Writing the input map to the outppf
@@ -225 +225 @@
               SphereHEALPix<float> betaMap(hp_nside);
               {
-                FITS_SphereHEALPix<float> fiosBM(flnm);
-                betaMap = (SphereHEALPix<float>)fiosBM;
+                FITS_SphereHEALPix<float> fiosBM(&betaMap);
+                fiosBM.ReadF(flnm,2);
+              }
+              if (printlev > 2) {
+                MeanSig(betaMap.DataBlock(), moy, sig );
+                cout << " MeanSig for Beta map - Mean= " << moy << " Sigma= " << sig << endl;
               }
               if (debuglev > 4) {  // Writing the input map to the outppf
@@ -248 +252 @@
                   
                   betaMap.Resize(8);
-                  MeanSig(intBetaMap.DataBlock(), moy, sig );
-                  if (printlev > 4) 
-                    cout << "....BetaFits... MeanSig for intBetaMap - Mean= " 
-                         << moy << " Sigma= " << sig << endl;
+                  if (printlev > 4)
+                    {
+                      MeanSig(intBetaMap.DataBlock(), moy, sig );
+                      cout << "....BetaFits... MeanSig for intBetaMap - Mean= " 
+                           << moy << " Sigma= " << sig << endl;
+                    }
                   // add the integrated beta map 
                   addComponentBeta(outgs,ings,intBetaMap, K);