Changeset 2587 in Sophya for trunk/SophyaLib/TArray

Timestamp:

Jul 30, 2004, 3:31:39 PM (21 years ago)

Author:

ansari

Message:

Diverses petites optimisations ds les methodes operations/affectations de TArray - Reza 30 Juillet 2004

File:

• trunk/SophyaLib/TArray/tarray.cc (modified) (11 diffs)

trunk/SophyaLib/TArray/tarray.cc

@@ -533 +533 @@
     pe = Data();
     per = res.Data();
-    for(k=0; k<maxx; k++) { *per = (*pe)+x;  per++;  pe++; }
+    for(k=0; k<maxx; k++)  *per++ = *pe++ + x;  
   }
   else {    // Non regular data spacing ...
@@ -577 +577 @@
     per = res.Data();
     if (!fginv) 
-      for(k=0; k<maxx; k++) { *per = (*pe)-x;  per++;  pe++; }
+      for(k=0; k<maxx; k++)  *per++ = *pe++ - x;  
     else 
-      for(k=0; k<maxx; k++) { *per = x-(*pe);  per++;  pe++; }
+      for(k=0; k<maxx; k++)  *per++ = x - *pe++; 
   }
   else {    // Non regular data spacing ...
@@ -624 +624 @@
     pe = Data();
     per = res.Data();
-    for(k=0; k<maxx; k++) { *per = (*pe)*x;  per++;  pe++; }
+    for(k=0; k<maxx; k++)  *per++ = *pe++ * x;  
   }
   else {    // Non regular data spacing ...
@@ -670 +670 @@
     per = res.Data();
     if (!fginv) 
-      for(k=0; k<maxx; k++) { *per = (*pe)/x;  per++;  pe++; }
+      for(k=0; k<maxx; k++)  *per++ = *pe++ / x;  
     else 
-      for(k=0; k<maxx; k++) { *per = x/(*pe);  per++;  pe++; }
+      for(k=0; k<maxx; k++) *per++ = x / *pe++; 
   }
   else {    // Non regular data spacing ...
@@ -715 +715 @@
     pe = Data();
     per = res.Data();
-    for(k=0; k<maxx; k++) { *per = -(*pe);  per++;  pe++; }
+    for(k=0; k<maxx; k++)  *per++ = -(*pe++);  
   }
   else {    // Non regular data spacing ...
@@ -771 +771 @@
     pea = a.Data();
     per = res.Data();
-    for(k=0; k<maxx;  k++, pe++, pea++, per++)  *per = *pe + *pea ;
+    //    for(k=0; k<maxx;  k++, pe++, pea++, per++)  *per = *pe + *pea ;
+    for(k=0; k<maxx;  k++)  *per++ = *pe++ + *pea++ ;
   }
   else {    // Non regular data spacing ...
@@ -839 +840 @@
     per = res.Data();
     if (!fginv) 
-      for(k=0; k<maxx;  k++, pe++, pea++, per++)  *per = *pe - *pea ;
+      for(k=0; k<maxx;  k++)  *per++ = *pe++ - *pea++ ;
     else 
-      for(k=0; k<maxx;  k++, pe++, pea++, per++)  *per = *pea - *pe ;
+      for(k=0; k<maxx;  k++)  *per++ = *pea++ - *pe++ ;
   }
   else {    // Non regular data spacing ...
@@ -910 +911 @@
     pea = a.Data();
     per = res.Data();
-    for(k=0; k<maxx;  k++, pe++, pea++, per++ )  *per = (*pe) * (*pea) ;
+    for(k=0; k<maxx;  k++)  *per++ = *pe++ * *pea++ ;
   }
   else {    // Non regular data spacing ...
@@ -981 +982 @@
     if(divzero) {
       if (!fginv) 
-        for(k=0; k<maxx;  k++, pe++, pea++, per++ )  
-          if (*pea==(T)0) *per = (T)0;  else   *per = *pe / *pea ;
+        for(k=0; k<maxx;  k++)  
+          if (*pea==(T)0) *per = (T)0;  else   *per++ = *pe++ / *pea++ ;
       else 
-        for(k=0; k<maxx;  k++, pe++, pea++, per++ )  
-          if (*pe==(T)0)  *per = (T)0;  else   *per = *pea / *pe ;
+        for(k=0; k<maxx;  k++)  
+          if (*pe==(T)0)  *per = (T)0;  else   *per++ = *pea++ / *pe++ ;
     }
     else {
       if (!fginv) 
-        for(k=0; k<maxx;  k++, pe++, pea++, per++ )  *per = *pe / *pea ;
+        for(k=0; k<maxx;  k++)  *per++ = *pe++ / *pea++ ;
       else
-        for(k=0; k<maxx;  k++, pe++, pea++, per++ )  *per = *pea / *pe ;
+        for(k=0; k<maxx;  k++)  *per = *pea++ / *pe++ ;
     }
   }
@@ -1052 +1053 @@
   sa_size_t j,k;
   if (smo && (AvgStep() > 0) && (a.AvgStep() > 0) )   {  // regularly spaced elements
-    sa_size_t step = AvgStep();
-    sa_size_t stepa = a.AvgStep();
-    sa_size_t maxx = totsize_*step;
-    pe = Data();
-    pea = a.Data();
-    for(k=0;  k<maxx;  k+=step, pe+=step, pea+=stepa )  *pe = *pea ;
+    if (IsPacked() && a.IsPacked())  memcpy(Data(), a.Data(), totsize_*sizeof(T)); // Packed arrays
+    else {
+      sa_size_t step = AvgStep();
+      sa_size_t stepa = a.AvgStep();
+      sa_size_t maxx = totsize_*step;
+      pe = Data();
+      pea = a.Data();
+      for(k=0;  k<maxx;  k+=step, pe+=step, pea+=stepa )  *pe = *pea ;
+    }
   }
   else {    // Non regular data spacing ...
@@ -1129 +1133 @@
     pe = Data();
     pea = a.Data();
-    for(k=0; k<maxx; k++, pe++, pea++)  res += (*pe)*(*pea);
+    for(k=0; k<maxx; k++)  res += *pe++ * *pea++;
   }
   else {    // Non regular data spacing ...