Changeset 813 in Sophya for trunk/SophyaLib/TArray/matharr.cc

Timestamp:

Apr 5, 2000, 5:44:19 PM (25 years ago)

Author:

ansari

Message:

Correction bug/amelioarions TArray,TMatrix,TVector - Reza 5/4/2000

File:

• trunk/SophyaLib/TArray/matharr.cc (modified) (3 diffs)

trunk/SophyaLib/TArray/matharr.cc

@@ -27 +27 @@
     uint_8 step = a.Step(ka);
     uint_8 gpas = a.Size(ka)*step;
-    for(j=0; j<a.Size(); j += a.Size(ka))  {
-      pe = a.DataBlock().Begin()+a.Offset(j);
+    uint_8 naxa = a.Size()/a.Size(ka);
+    for(j=0; j<naxa; j++)  {
+      pe = a.DataBlock().Begin()+a.Offset(ka,j);
       for(k=0; k<gpas; k+=step)  pe[k] = (T)(f((double)pe[k]));
     }
@@ -52 +53 @@
     uint_8 step = a.Step(ka);
     uint_8 gpas = a.Size(ka)*step;
-    for(j=0; j<a.Size(); j += a.Size(ka))  {
-      pe = a.DataBlock().Begin()+a.Offset(j);
+    uint_8 naxa = a.Size()/a.Size(ka);
+    for(j=0; j<naxa; j++)  {
+      pe = a.DataBlock().Begin()+a.Offset(ka,j);
       for(k=0; k<gpas; k+=step)  pe[k] = (T)(f((float)pe[k]));
     }
@@ -102 +104 @@
     uint_8 step = a.Step(ka);
     uint_8 gpas = a.Size(ka)*step;
-    for(j=0; j<a.Size(); j += a.Size(ka))  {
-      pe = a.DataBlock().Begin()+a.Offset(j);
+    uint_8 naxa = a.Size()/a.Size(ka);
+    for(j=0; j<naxa; j++)  {
+      pe = a.DataBlock().Begin()+a.Offset(ka,j);
       for(k=0; k<gpas; k+=step) { 
         valok = (double) pe[k];