Changeset 3264 in Sophya for trunk/SophyaProg

Timestamp:

Jun 11, 2007, 9:47:11 PM (18 years ago)

Author:

ansari

Message:

amelioration / optimisation des operations arrdl/arrmf dans zthr.cc , Reza 11/06/2007

File:

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

trunk/SophyaProg/Tests/zthr.cc

@@ -102 +102 @@
   ztarg * za = (ztarg *)arg;
   sa_size_t vsz = za->M*za->M;
-
-  cout << ">>>> arrdl-funzt(ThId=" << za->thid << " VecSz=M*M " << vsz  
-       << " V2=DLO4(V1) , NOp ~= 10*vsz=" << 10*vsz << endl;
+  sa_size_t EXS = 10;
+
+  cout << ">>>> arrdl-funzt(ThId=" << za->thid << " DataBlock_Sz=M*M " << vsz  
+       << " V2=DLO4(V1) , NOp ~= 10*10*vsz=" << 10*EXS*vsz << endl;
   
   TVector<r_8> v1(vsz), v2(vsz); 
-  v1 = RegularSequence(1.,0.001);
+  TVector<r_8> coeff(EXS); 
+  coeff = RandomSequence();
+  //  v1 = RegularSequence(1.,0.001);  --- ATTENTION , couteux en temps
+  //  NDataBlock<r_8> v1(vsz, false), v2(vsz, false); 
+  for(sa_size_t i=0; i<vsz; i++)   v1(i) = i*0.001; 
   char buff[128];
   sprintf(buff, "arrdl-funzt(ThId=%d) EndOfInit", za->thid); 
   PrtTim(buff);
-  for(sa_size_t i=0; i<vsz; i++) {
-    double x = v1(i);
-    v2(i) = x*(1+x*(x/2.+x*(x/6.+x*x/24.)));
+
+  double c2 = 0.5;
+  double c3 = 1./6.;
+  double c4 = 1./24;
+  for(sa_size_t k=0; k<EXS; k++) {
+    for(sa_size_t i=0; i<vsz; i++) {
+      register double x = v1(i)*coeff(k); 
+      v2(i) = x*(1+x*(x*c2+x*(x*c3+x*x*c4)));
+    }
   }
   sprintf(buff, "arrdl-funzt(ThId=%d) EndOfOper", za->thid); 
@@ -126 +137 @@
 
   cout << ">>>> arrmf-funzt(ThId=" << za->thid << " VecSz=M*M " << vsz  
-       << " V2=Sin(V1) , NOp ~= 10*vsz=" << 10*vsz << endl;
+       << " V2=Sin(V1) , NOp ~= 50*vsz=" << 50*vsz << endl;
   
   TVector<r_8> v1(vsz), v2(vsz); 
-  v1 = RegularSequence(1.,0.001);
-  char buff[128];
-  sprintf(buff, "arrmf-funzt(ThId=%d) EndOfInit", za->thid); 
-  PrtTim(buff);
+  //--  v1 = RegularSequence(1.,0.001);  COUTEUX en TCPU
+  for(sa_size_t i=0; i<vsz; i++)  v1(i) = i*0.001; 
+  char buff[128];
+  sprintf(buff, "arrmf-funzt(ThId=%d) EndOfInit", za->thid);
+  PrtTim(buff);
+
   v2 = Sin(v1);
-
+  v2 += Cos(v1);
   sprintf(buff, "arrmf-funzt(ThId=%d) EndOfOper", za->thid); 
   PrtTim(buff);