Changeset 3873 in Sophya for trunk/SophyaProg/PMixer/extractRS.cc

Timestamp:

Sep 2, 2010, 5:38:01 PM (15 years ago)

Author:

cmv

Message:

Suite a la disparition de frand01, modification pour compilation:

frand01 -> drand01

+ Modifs pour palier une certaine vision poetique de l'informatique !!!

if( double == int) blabla !!!!

cmv, 02/09/2010

File:

• trunk/SophyaProg/PMixer/extractRS.cc (modified) (5 diffs)

trunk/SophyaProg/PMixer/extractRS.cc

@@ -43 +43 @@
   string msg;
   int rc = 0;
-  RadSpectra * es = NULL;
-  SpectralResponse * sr = NULL;
+  //RadSpectra * es = NULL;
+  //SpectralResponse * sr = NULL;
   double moy, sig;
   
@@ -91 +91 @@
     // Loop over sky frequencies components
     int maxOfSkyComp = nskycomp;
-    int nbOfPix;
+    int nbOfPix = 0;
     int sk;
     Vector freq(maxOfSkyComp);
@@ -138 +138 @@
         if (printlev > 2) cout << " Reading Input FITS map " << nameOfFile << endl;
         ifstream from(nameOfFile);
-        char ch;
         string dum;
         from >> nb_source >> dum;
@@ -209 +208 @@
     for(i=0;i<OutputMap.NbPixels();i++)
       {
-        float rndval = frand01();
-        if(rndval == 1) rndval = frand01();
-        int rndvalPix = rndval*OutputMap.NbPixels();
-        int outpix = rndval*nbOfPix; 
+        double rndval = drand01();
+        //stupid!!!!! if(rndval == 1) rndval = drand01();
+        //unsused!!!!! int rndvalPix = rndval*OutputMap.NbPixels();
+        int outpix = int(rndval*nbOfPix);
+        if(outpix >= nbOfPix) outpix = nbOfPix - 1; // add clean size protection
         OutputMap(i) = outCoeff(outpix);
       }
@@ -234 +234 @@
     for(i=0;i<nb_source;i++)
       {
-        float rndval = frand01();
-        int outPix = rndval*OutputMap.NbPixels(); 
+        double rndval = drand01();
+        int outPix = int(rndval*OutputMap.NbPixels());
+        if(outPix >= OutputMap.NbPixels()) outPix = OutputMap.NbPixels() - 1; // add clean size protection
         SourceMap(outPix) = sourceCoeff(i); 
         OutputMap(outPix) = OutputMap(i)+ SourceMap(outPix);