Changeset 4054 in Sophya for trunk/Eval/JET/tjet.cc

Timestamp:

Mar 17, 2012, 11:13:26 PM (14 years ago)

Author:

ansari

Message:

ajout programme telacc.cc pour test rapidite acces aux elements, Reza 17/03/2012

File:

• trunk/Eval/JET/tjet.cc (modified) (14 diffs)

trunk/Eval/JET/tjet.cc

@@ -20 +20 @@
 //   Test d'expression template sur tableaux 
 //       Reza - Avril 2003
+//       revisite : Mars 2012 
 //--------------------------------------------------------
 
@@ -81 +82 @@
     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 +96 @@
     // 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 +109 @@
 
     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 +118 @@
     // 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 +130 @@
 
     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 +141 @@
 
     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 +162 @@
     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 +175 @@
     // 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 +195 @@
 
     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 +204 @@
 
     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 +214 @@
 
     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 +232 @@
     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 +244 @@
     // 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];