Changeset 3512 in Sophya for trunk

Timestamp:

Aug 8, 2008, 3:12:51 PM (17 years ago)

Author:

ansari

Message:

Correction,adaptation suite modifSphericalTransformServer , Reza 08/08/2008

File:

• trunk/SophyaProg/PrgMap/cl2map.cc (modified) (3 diffs)

trunk/SophyaProg/PrgMap/cl2map.cc

@@ -85 +85 @@
 public :
 static void BuildMap(string & infile, string & outfile, int msph, double beam, 
-                     bool fgfitsin, bool fgfitsout)
+                     bool fgfitsin, bool fgfitsout, bool fgair)
 {
   double deg2rad =  M_PI/180.;
@@ -115 +115 @@
        << "' msph= " << msph << ") beam=" << beam*minute2rad << " rad" << endl;
   // on cree la carte 
-  SphericalTransformServer<T> sphtr;
-  sphtr.GenerateFromCl(sph, msph, clvec, beam*minute2rad);
+  RandomGenerator rg;
+  if (fgair) {
+    cout << " Setting random number generator seed (using time) " << endl; 
+    rg.AutoInit(); 
+  }
+  SphericalTransformServer<T> sphtr(rg);
+  sphtr.GenerateFromCl(sph, msph, clvec, beam);
 
   cout << " Output map : NbPixels= " <<  sph.NbPixels() << " NSide= " 
@@ -186 +191 @@
     InitTim();
     SophyaInit();
-    if (fginirand) {
-      cout << " Setting random number generator seed (using time) " << endl; 
-      Auto_Ini_Ranf(1);
-    }
     if (fgr4) {
       cout << " Power spectrum C_l (r_4) --> SphereHEALPix<r_4> (float) " << endl;
-      _Cl2Map<r_4>::BuildMap(infile, outfile, msph, beam, fgfitsin, fgfitsout);
+      _Cl2Map<r_4>::BuildMap(infile, outfile, msph, beam, fgfitsin, fgfitsout, fginirand);
     }
     else {
       cout << " Power spectrum C_l (r_8) --> SphereHEALPix<r_8> (double) " << endl;
-      _Cl2Map<r_8>::BuildMap(infile, outfile, msph, beam, fgfitsin, fgfitsout);
+      _Cl2Map<r_8>::BuildMap(infile, outfile, msph, beam, fgfitsin, fgfitsout, fginirand);
     }
   }