Changeset 1683 in Sophya for trunk/SophyaProg/PrgMap

Timestamp:

Oct 11, 2001, 5:22:03 PM (24 years ago)

Author:

lemeur

Message:

methode iterative pour analyse harmonique

File:

• trunk/SophyaProg/PrgMap/map2cl.cc (modified) (8 diffs)

trunk/SophyaProg/PrgMap/map2cl.cc

@@ -50 +50 @@
   else {
     cout << " map2cl : Spherical harmonics analysis - HEALPix map -> Power spectrum C_l \n" 
-         << " Usage: map2cl [-float/-r4] [-lmax lval] [-thetacut dtdeg] \n" 
+         << " Usage: map2cl [-float/-r4] [-lmax lval] [-thetacut dtdeg] [-iter_order lval]\n" 
          << "        [-fitsin] [-fitsout] InFileName OutFileName \n" 
          << "   -float (-r4): single precision C_l and map (default = double) \n"
@@ -56 +56 @@
          << "   -thetacut dtdeg : Symmetric delta-theta cut (in degree) along equator \n" 
          << "                    (default=0 -> no cut)\n"
+         << "   -iter_order lval : 1,2,3,4... order of an iterative analysis , 3rd order is usually optimal (default=0 -> standard analysis)\n"
          << "   -fitsout: Select the FITS format for the output map (default PPF format) \n"
          << "   -fitsin : Select the FITS format for the input vector (default PPF format) \n"
@@ -69 +70 @@
 public :
 static void ComputeCl(string & infile, string & outfile, int lmax, double tcut, 
-                      bool fgfitsin, bool fgfitsout)
+                      int iterationOrder, bool fgfitsin, bool fgfitsout)
 {
   double deg2rad =  M_PI/180.;
   double minute2rad =  M_PI/(180.*60.);
-
+ 
   SphereHEALPix<T> sph;
   if (fgfitsin) {
@@ -84 +85 @@
     PInPersist ppi(infile);
     ppi >> sph;
-  }
+  } 
 
   cout << " Input map : NbPixels= " <<  sph.NbPixels() << " NSide= " 
@@ -92 +93 @@
   double ctcut = (tcut < 1.e-19) ? 0. : cos((90.-tcut)*deg2rad);
   cout << "--- Calling  DecomposeToCl() (lmax= " << lmax 
-       << " cos_theta_cut= " << ctcut << ") theta_cut=" << tcut << " deg" << endl;
+       << " cos_theta_cut= " << ctcut << " iter_order= " << iterationOrder << ") theta_cut=" << tcut << " deg"  << endl;
   // Decomposition de la carte en C_l
   SphericalTransformServer<T> sphtr;
-  TVector<T> clvec = sphtr.DecomposeToCl(sph, lmax, ctcut);
+  TVector<T> clvec = sphtr.DecomposeToCl(sph, lmax, ctcut, iterationOrder);
  
   T min, max;
@@ -125 +126 @@
   
   int lmax = 255;
+  int iterationOrder = 0;
   double tcut = 0.;
   string infile;
@@ -142 +144 @@
       if (k == narg-1) Usage(true);  // -thetacut est suivi d'un argument 
       tcut = atof(arg[k+1]);  k++;       // k++ pour sauter au suivant
+    }
+    else if (strcmp(arg[k], "-iter_order") == 0)  {
+      if (k == narg-1) Usage(true);  // -iter_order est suivi d'un argument 
+      iterationOrder = atof(arg[k+1]);  k++;   // k++ pour sauter au suivant
     }
     else if (strcmp(arg[k], "-fitsin") == 0) {
@@ -166 +172 @@
     if (fgr4) {
       cout << " SphereHEALPix<r_4>  --> Power spectrum C_l<r_4> (float)" << endl;
-      _Map2Cl<r_4>::ComputeCl(infile, outfile, lmax, tcut, fgfitsin, fgfitsout);
+      _Map2Cl<r_4>::ComputeCl(infile, outfile, lmax, tcut, iterationOrder, fgfitsin, fgfitsout);
     }
     else {
       cout << " SphereHEALPix<r_8>  --> Power spectrum C_l<r_8> (double)" << endl;
-      _Map2Cl<r_8>::ComputeCl(infile, outfile, lmax, tcut, fgfitsin, fgfitsout);
+      _Map2Cl<r_8>::ComputeCl(infile, outfile, lmax, tcut, iterationOrder,fgfitsin, fgfitsout);
     }
   }