Changeset 3606 in Sophya for trunk/SophyaProg/Tests/tmtrnd.cc

Timestamp:

Apr 29, 2009, 12:16:11 PM (16 years ago)

Author:

ansari

Message:

Modif tmtrnd.cc pour utilisation classe ThSDR48RandGen , Reza 29/04/2009

File:

• trunk/SophyaProg/Tests/tmtrnd.cc (modified) (6 diffs)

trunk/SophyaProg/Tests/tmtrnd.cc

@@ -29 +29 @@
 #include "srandgen.h"
 #include "stsrand.h"
-
+#include "randr48.h"
+
+
+// Choix du nom de la classe thread-safe a tests 
+#define TC_RandomGenerator ThSDR48RandGen
 
 //---------------------------------------------------------------------------------
@@ -62 +66 @@
   Timer tm(nom_.c_str());
   cout << "MTRnd::run() - Nom= " << nom_ << " vv.Size()= " << vv_.Size() << endl;
-  RandomGenerator rgen;
+  TC_RandomGenerator rgen;
   if (fgg_) 
     for(sa_size_t k=0; k<vv_.Size(); k++) vv_(k) = rgen.Gaussian();
@@ -95 +99 @@
   cout << "[1] f1_tmtrnd/starting, VSZ= " << VSZ << " SeqSz=" << seqsz << endl;
   TVector<r_8> DATA(VSZ);
-  RandomGenerator rg(seqsz, (seqsz==0)?false:true);  
+  TC_RandomGenerator rg(seqsz, (seqsz==0)?false:true);  
   unsigned short seed[3];
-  rg.GetSeed(seed, 2);  // pour imprimer l'etat du generateur
+//  rg.GetSeed(seed,2);  // pour imprimer l'etat du generateur
   {
     Timer tm("f0_tmtrnd-GauRnd/drand01");
@@ -111 +115 @@
   POutPersist porg("rg1.ppf");
   porg << rg;
-  rg.GetSeed(seed, 2);  // pour imprimer l'etat du generateur
+//  rg.GetSeed(seed);  // pour imprimer l'etat du generateur
   cout << "[4] f1_tmtrnd/  RandGen saved to file rg1.ppf " << endl;
   return 0; 
@@ -120 +124 @@
   cout << "[1] f2_tmtrnd/starting, VSZ= " << VSZ << " InPPFName=" << inppf << endl;
   TVector<r_8> DATA(VSZ);
-  RandomGenerator rg;
+  TC_RandomGenerator rg;
   PInPersist pirg(inppf);
   pirg >> rg;
   cout << "[1.b] f2_tmtrnd/ RandGen read from InPPF OK " << endl;
   unsigned short seed[3];
-  rg.GetSeed(seed, 2);  // pour imprimer l'etat du generateur
+//  rg.GetSeed(seed,2);  // pour imprimer l'etat du generateur
   {
     Timer tm("f0_tmtrnd-GauRnd/drand01");
@@ -139 +143 @@
   POutPersist porg("rg2.ppf");
   porg << rg;
-  rg.GetSeed(seed, 2);  // pour imprimer l'etat du generateur
+//  rg.GetSeed(seed,2);  // pour imprimer l'etat du generateur
   cout << "[4] f1_tmtrnd/  RandGen saved to file rg2.ppf " << endl;
   return 0;