Index: trunk/Eval/JET/tjet.cc
===================================================================
--- trunk/Eval/JET/tjet.cc	(revision 2939)
+++ trunk/Eval/JET/tjet.cc	(revision 4054)
@@ -20,4 +20,5 @@
 //   Test d'expression template sur tableaux 
 //       Reza - Avril 2003
+//       revisite : Mars 2012 
 //--------------------------------------------------------
 
@@ -81,6 +82,6 @@
     SimpleMatrix<r_8> m5(nrow, ncol);
     for(k=0; k<N; k++) {
-      double xxr = frand01();
-      double yyr = 2.*frand01();
+      double xxr = drand01();
+      double yyr = 2.*drand01();
       for(i=0; i<nrow; i++)
 	for(j=0; j<ncol; j++) {
@@ -95,7 +96,7 @@
     // Calcul m1*c1 + m2*c2 + m3*c3;
     for(k=0; k<N; k++) {
-      double c1 = frand01() + 1.2;
-      double c2 = frand01() + 3.5;
-      double c3 = frand01() + 6.7;
+      double c1 = drand01() + 1.2;
+      double c2 = drand01() + 3.5;
+      double c3 = drand01() + 6.7;
       m5 = m1;  m5.MulCst(c1); 
       m4 = m2;  m4.MulCst(c2);
@@ -108,7 +109,7 @@
 
     for(k=0; k<N; k++) {
-      double c1 = frand01() + 1.2;
-      double c2 = frand01() + 3.5;
-      double c3 = frand01() + 6.7;
+      double c1 = drand01() + 1.2;
+      double c2 = drand01() + 3.5;
+      double c3 = drand01() + 6.7;
       m5 = m1*c1+m2*c2+m3*c3;
     }
@@ -117,7 +118,7 @@
     // Calcul m1*c1 + m1*(m2+c2) + m2*c3; 
     for(k=0; k<N; k++) {
-      double c1 = frand01() + 1.2;
-      double c2 = frand01() + 3.5;
-      double c3 = frand01() + 6.7;
+      double c1 = drand01() + 1.2;
+      double c2 = drand01() + 3.5;
+      double c3 = drand01() + 6.7;
       m3 = m1;  m3.MulCst(c1); 
       m4 = m2;  m4.AddCst(c2); m4.MulElt(m1);
@@ -129,7 +130,7 @@
 
     for(k=0; k<N; k++) {
-      double c1 = frand01() + 1.2;
-      double c2 = frand01() + 3.5;
-      double c3 = frand01() + 6.7;
+      double c1 = drand01() + 1.2;
+      double c2 = drand01() + 3.5;
+      double c3 = drand01() + 6.7;
       //      m4 = SMExprAdd< r_8, SMExprMtx<r_8>, SMExprMtx<r_8> >( SMExprMtx<r_8>(m1), SMExprMtx<r_8> (m2) );  
     // Calcul m1*c1 + m1*(m2+c2) + m2*c3; 
@@ -140,7 +141,7 @@
 
     for(k=0; k<N; k++) {
-      double c1 = frand01() + 1.2;
-      double c2 = frand01() + 3.5;
-      double c3 = frand01() + 6.7;
+      double c1 = drand01() + 1.2;
+      double c2 = drand01() + 3.5;
+      double c3 = drand01() + 6.7;
     // Calcul m1*c1 + m1*(m2+c2) + m2*c3; 
       m5 = m1*c1 + m2 + c2*Sin(c3*m3+m1);
@@ -161,6 +162,6 @@
     TMatrix<r_8> m5(nrow, ncol);
     for(k=0; k<N; k++) {
-      double xxr = frand01();
-      double yyr = 2.*frand01();
+      double xxr = drand01();
+      double yyr = 2.*drand01();
       for(i=0; i<nrow; i++)
 	for(j=0; j<ncol; j++) {
@@ -174,7 +175,7 @@
     // Calcul m1*c1 + m2*c2 + m3*c3;
     for(k=0; k<N; k++) {
-      double c1 = frand01() + 1.2;
-      double c2 = frand01() + 3.5;
-      double c3 = frand01() + 6.7;
+      double c1 = drand01() + 1.2;
+      double c2 = drand01() + 3.5;
+      double c3 = drand01() + 6.7;
       /* Changement interface SOPHYA  : 2004-2005  
 	 m5 = m1;  m5.Mul(c1); 
@@ -194,7 +195,7 @@
 
     for(k=0; k<N; k++) {
-      double c1 = frand01() + 1.2;
-      double c2 = frand01() + 3.5;
-      double c3 = frand01() + 6.7;
+      double c1 = drand01() + 1.2;
+      double c2 = drand01() + 3.5;
+      double c3 = drand01() + 6.7;
       m5 = m1*c1+m2*c2+m3*c3;
     }
@@ -203,7 +204,7 @@
 
     for(k=0; k<N; k++) {
-      double c1 = frand01() + 1.2;
-      double c2 = frand01() + 3.5;
-      double c3 = frand01() + 6.7;
+      double c1 = drand01() + 1.2;
+      double c2 = drand01() + 3.5;
+      double c3 = drand01() + 6.7;
     // Calcul m1*c1 + m1*(m2+c2) + m2*c3; 
       m5 = m1*c1 + m1 && (m2+c2) + m2*c3;
@@ -213,7 +214,7 @@
 
     for(k=0; k<N; k++) {
-      double c1 = frand01() + 1.2;
-      double c2 = frand01() + 3.5;
-      double c3 = frand01() + 6.7;
+      double c1 = drand01() + 1.2;
+      double c2 = drand01() + 3.5;
+      double c3 = drand01() + 6.7;
       m5 = m1*c1 + m2 + c2*Sin(c3*m3+m1);
     }
@@ -231,6 +232,6 @@
     r_8 * m5 = new r_8[nrow*ncol];
     for(k=0; k<N; k++) {
-      double xxr = frand01();
-      double yyr = 2.*frand01();
+      double xxr = drand01();
+      double yyr = 2.*drand01();
       for(i=0; i<nrow*ncol; i++) {
 	  m1[i] = k*300+10.*i+j+xxr;	
@@ -243,7 +244,7 @@
     // Calcul m1*c1 + m2*c2 + m3*c3;
     for(k=0; k<N; k++) {
-      double c1 = frand01() + 1.2;
-      double c2 = frand01() + 3.5;
-      double c3 = frand01() + 6.7;
+      double c1 = drand01() + 1.2;
+      double c2 = drand01() + 3.5;
+      double c3 = drand01() + 6.7;
       for(i=0; i<nrow*ncol; i++) 
 	m5[i] = c1*m1[i]+c2*m2[i]+c3*m3[i];
