Changeset 41861 for branches/eam_branches/ipp-dev-20210817/psModules/src/detrend/pmPattern.c

Timestamp:

Oct 25, 2021, 5:04:15 PM (5 years ago)

Author:

eugene

Message:

add some asserts to trap memory overflow

File:

• branches/eam_branches/ipp-dev-20210817/psModules/src/detrend/pmPattern.c (modified) (5 diffs)

branches/eam_branches/ipp-dev-20210817/psModules/src/detrend/pmPattern.c

@@ -418 +418 @@
           for (int x = 0; x < numCols; x++) {
             image->data.F32[y][x] += solution->data.F32[y];
-            psTrace("pattern",5,"B: %d %d %g\n",x,y,solution->data.F32[x]);
+            psTrace("pattern",5,"B: %d %d %g\n",x,y,solution->data.F32[y]);
           }
           corr->data.F64[y][0]  -= solution->data.F32[y];
@@ -444 +444 @@
           for (int x = 0; x < numCols; x++) {
             image->data.F32[y][x] += solution->data.F32[y] * indices->data.F32[x];
-            psTrace("pattern",5,"C: %d %d %g %g\n",x,y,solution->data.F32[x],indices->data.F32[x]);
+            // XXX EAM : this was [x] which is wrong
+            psTrace("pattern",5,"C: %d %d %g %g\n",x,y,solution->data.F32[y],indices->data.F32[x]);
           }
           corr->data.F64[y][1]  -= solution->data.F32[y] ;
@@ -703 +704 @@
         }
 
+        psAssert (solution->n == numCols, "oops");
         for (int x = 0; x < numCols; x++) {
             image->data.F32[y][x] -= solution->data.F32[x];
@@ -747 +749 @@
         }
       } else {
+        psAssert (solution->n == numRows, "oops");
         for (int y = 0; y < numRows; y++) {
           for (int x = 0; x < numCols; x++) {
             image->data.F32[y][x] += solution->data.F32[y];
-            psTrace("pattern",5,"B: %d %d %g\n",x,y,solution->data.F32[x]);
+            // XXX EAM : this was [x], which is wrong
+            psTrace("pattern",5,"B: %d %d %g\n",x,y,solution->data.F32[y]);
           }
           corr->data.F64[y][0]  -= solution->data.F32[y];
@@ -771 +775 @@
         psErrorClear();
       } else {
+        psAssert (solution->n == numRows, "oops");
         for (int y = 0; y < numRows; y++) {
           for (int x = 0; x < numCols; x++) {
             image->data.F32[y][x] += solution->data.F32[y] * indices->data.F32[x];
-            psTrace("pattern",5,"C: %d %d %g %g\n",x,y,solution->data.F32[x],indices->data.F32[x]);
+            // XXX EAM : this was set to [x] which is wrong (numCols > numRows)
+            psTrace("pattern",5,"C: %d %d %g %g\n",x,y,solution->data.F32[y],indices->data.F32[x]);
           }
           corr->data.F64[y][1]  -= solution->data.F32[y] ;