Changeset 13776 for trunk/psModules/src/detrend/pmFringeStats.c

Timestamp:

Jun 12, 2007, 5:42:58 PM (19 years ago)

Author:

magnier

Message:

fixed error copying

File:

• trunk/psModules/src/detrend/pmFringeStats.c (modified) (2 diffs)

trunk/psModules/src/detrend/pmFringeStats.c

@@ -502 +502 @@
         memcpy(&newRegions->mask->data.U8[offset], regions->mask->data.U8, regions->mask->n * sizeof(psU8));
         memcpy(&newStats->f->data.F32[offset], fringe->f->data.F32, fringe->f->n * sizeof(psF32));
-        memcpy(&newStats->df->data.F32[offset], fringe->f->data.F32, fringe->df->n * sizeof(psF32));
+        memcpy(&newStats->df->data.F32[offset], fringe->df->data.F32, fringe->df->n * sizeof(psF32));
         if (x0 && y0) {
             for (long j = offset; j < offset + regions->x->n; j++) {
@@ -952 +952 @@
     }
 
+# if (0)
+    // Write fringe data to file for a test
+    FILE *f = fopen ("fringe.dat", "w");
+    for (int j = 0; j < numRegions; j++) {
+        if (mask->data.U8[j]) continue;
+        fprintf (f, "%d %f %f ", j, science->f->data.F32[j], science->df->data.F32[j]);
+        for (int i = 0; i < fringes->n; i++) {
+            pmFringeStats *fringe = fringes->data[i]; // The fringe of interest
+            fprintf (f, "%f  ", fringe->f->data.F32[j]);
+        }
+        fprintf (f, "\n");
+    }
+    fclose (f);
+# endif
+
     // Get rid of the extreme outliers by assuming most of the points are somewhat clustered
     psVectorStats(median, science->f, NULL, NULL, 0);