Changeset 3613 in Sophya for trunk/SophyaProg

Timestamp:

Apr 30, 2009, 7:07:46 PM (16 years ago)

Author:

ansari

Message:

Adaptation a l'introduction de la suite des classes RandomGeneratorInterface et Cie , Reza 30/04/2009

Location:

trunk/SophyaProg

Files:

• PrgMap/cl2map.cc (modified) (2 diffs)
• Tests/tmtrnd.cc (modified) (1 diff)
• Tests/tmtrnd3.cc (modified) (9 diffs)

trunk/SophyaProg/PrgMap/cl2map.cc

@@ -14 +14 @@
 #include "fitsspherehealpix.h"
 #include "fitstarray.h"
+
+#include "randr48.h"
 
 /*!
@@ -115 +117 @@
        << "' msph= " << msph << ") beam=" << beam*minute2rad << " rad" << endl;
   // on cree la carte 
-  RandomGenerator rg;
+  ThSDR48RandGen rg;
   if (fgair) {
     cout << " Setting random number generator seed (using time) " << endl; 

trunk/SophyaProg/Tests/tmtrnd.cc

@@ -28 +28 @@
 
 #include "srandgen.h"
-#include "stsrand.h"
 #include "randr48.h"
 
 
-// Choix du nom de la classe thread-safe a tests 
+// Choix du nom de la classe thread-safe a tester 
 #define TC_RandomGenerator ThSDR48RandGen
 

trunk/SophyaProg/Tests/tmtrnd3.cc

@@ -17 +17 @@
 #include "timing.h"
 
-#include "stsrand.h"
+#include "randr48.h"
 
 long Check_Redondances(TVector<r_8>& v);
@@ -55 +55 @@
   cout<<"MTRndTest::run() - Nom= "<<nom_<<" vv.Size()= "<<vv_.Size()
       <<" with threadsafe="<<thsafe_<<endl;
-  if(thsafe_) {  // On appele RandomGenerator Thread-Safe
-    RandomGenerator rgen(nbuff_,true);
+  if(thsafe_) {  // On appele ThSDR48RandGen Thread-Safe
+    ThSDR48RandGen rgen(nbuff_,true);
     for(sa_size_t k=0; k<vv_.Size(); k++) vv_(k) = rgen.Flat01();
-  } else {  // On appele RandomGenerator NO Thread-Safe
+  } else {  // On appele ThSDR48RandGen NO Thread-Safe
     // ATTENTION: le comportement est different selon qu'un thread
     // se termine avant ou apres que le suivant commence:
@@ -65 +65 @@
     // cela donne un comportement "tread_safe" APPARENT mais
     // dans la realite drand48() N'EST PAS thread_safe !!!!
-    RandomGenerator rgen(0,false);
+    ThSDR48RandGen rgen(0,false);
     for(sa_size_t k=0; k<vv_.Size(); k++) vv_(k) = rgen.Flat01();
   }
@@ -83 +83 @@
   if(NTH<=0) NTH = 1;
   cout<<"tmtrnd/arguments, VSZ= "<<VSZ<<" NTH= "<<NTH<<" NBUFF="<<NBUFF<<endl;
-  // Il ne faut pas que les buffers de RandomGenerator soient remplis
+  // Il ne faut pas que les buffers de ThSDR48RandGen soient remplis
   // au dela de ce qui est necessaire si on veut faire
   // la comparaison avec le mono thread
@@ -100 +100 @@
     //--- Init et set de l'aleatoire
     //---
-    RandomGenerator rgen_nothsf(0,false);
+    ThSDR48RandGen rgen_nothsf(0,false);
     unsigned short seed_16v[3];
     {
     cout<<"\n>>>>>>> Init de l'aleatoire"<<endl;
     rgen_nothsf.AutoInit(5);
-    rgen_nothsf.GetSeed(seed_16v,5);
+    rgen_nothsf.GetSeed(seed_16v);
     }
 
@@ -116 +116 @@
     {
     cout<<"\n>>>>>>> Remplissage en mono-thread avec drand48()"<<endl;
-    rgen_nothsf.Init(seed_16v,5);
+    rgen_nothsf.SetSeed(seed_16v);
     for(int i=0;i<VSZ;i++) DATA0(i) = drand48();
     PrtTim(" ");
@@ -122 +122 @@
 
     //---
-    //--- Remplissage par threads avec RandomGenerator thread safe
+    //--- Remplissage par threads avec ThSDR48RandGen thread safe
     //---
     TVector<r_8> DATA1(VSZ);
     {
-    cout<<"\n>>>>>>> Remplissage par threads avec RandomGenerator thread safe"<<endl;
-    rgen_nothsf.Init(seed_16v,5);
+    cout<<"\n>>>>>>> Remplissage par threads avec ThSDR48RandGen thread safe"<<endl;
+    rgen_nothsf.SetSeed(seed_16v);
     vector<MTRndTest *> vth1;
     cout<<"...tmtrnd/creating threads "<<endl;
@@ -147 +147 @@
 
     //---
-    //--- Remplissage par threads avec RandomGenerator SANS thread safe
+    //--- Remplissage par threads avec ThSDR48RandGen SANS thread safe
     //---
     TVector<r_8> DATA2(VSZ);
     {
-    cout<<"\n>>>>>>> Remplissage par threads avec RandomGenerator SANS thread safe"<<endl;
-    rgen_nothsf.Init(seed_16v,5);
+    cout<<"\n>>>>>>> Remplissage par threads avec ThSDR48RandGen SANS thread safe"<<endl;
+    rgen_nothsf.SetSeed(seed_16v);
     vector<MTRndTest *> vth2;
     cout << "......tmtrnd/creating threads " << endl;
@@ -172 +172 @@
 
     //---
-    //--- remplissage en mono-thread avec RandomGenerator
+    //--- remplissage en mono-thread avec ThSDR48RandGen
     //---
     TVector<r_8> DATA0R(VSZ);
     {
-    cout<<"\n>>>>>>> Remplissage en mono-thread avec RandomGenerator"<<endl;
-    rgen_nothsf.Init(seed_16v,5);
+    cout<<"\n>>>>>>> Remplissage en mono-thread avec ThSDR48RandGen"<<endl;
+    rgen_nothsf.SetSeed(seed_16v);
     for(int i=0;i<VSZ;i++) DATA0R(i) = rgen_nothsf.Flat01();
     PrtTim(" ");