Changeset 3451 in Sophya for trunk/SophyaProg

Timestamp:

Jan 31, 2008, 7:25:14 PM (18 years ago)

Author:

ansari

Message:

Init tableaux avec randoms rendu thread-safe ds zthr.cc - Reza , 31/01/2008

File:

• trunk/SophyaProg/Tests/zthr.cc (modified) (4 diffs)

trunk/SophyaProg/Tests/zthr.cc

@@ -10 +10 @@
 #include "matharr.h"
 #include "tarrinit.h"
+#include "stsrand.h"
 
 #include <stdlib.h>
@@ -60 +61 @@
 }
 
+void randgenths(Array& mx, double sig)
+{
+  RandomGenerator rg;
+  for(sa_size_t j=0; j<mx.SizeY(); j++)   
+    for(sa_size_t i=0; i<mx.SizeX(); i++) mx(i,j) = rg.Gaussian(sig);   
+}
+
 // --- fonction de test simple avec calcul matriciel (produit de 2 matrices double)
 void mtx_funzt(void *arg)
@@ -68 +76 @@
   sa_size_t m = za->M;
   Matrix a1(m,m), a2(m,m), mxprod; 
-  a1 = RandomSequence(RandomSequence::Gaussian, 0., 4.);
-  a2 = RandomSequence(RandomSequence::Gaussian, 0., 3.);
+//  a1 = RandomSequence(RandomSequence::Gaussian, 0., 4.);
+//  a2 = RandomSequence(RandomSequence::Gaussian, 0., 3.);
+  randgenths(a1,4.);
+  randgenths(a2,3.);
   char buff[128];
   sprintf(buff, "mtx-funzt(ThId=%d) EndOfInit", za->thid); 
@@ -109 +119 @@
   TVector<r_8> v1(vsz), v2(vsz); 
   TVector<r_8> coeff(EXS); 
-  coeff = RandomSequence();
+//  coeff = RandomSequence();
+  randgenths(coeff, 1.);
   //  v1 = RegularSequence(1.,0.001);  --- ATTENTION , couteux en temps
   //  NDataBlock<r_8> v1(vsz, false), v2(vsz, false);