Changeset 42392 for trunk/extsrc/gpcsw/gpcsrc/fits/burntool/persistfix.c

Timestamp:

Feb 14, 2023, 12:48:54 PM (3 years ago)

Author:

tdeboer

Message:

new version of the burntool programme

File:

• trunk/extsrc/gpcsw/gpcsrc/fits/burntool/persistfix.c (modified) (3 diffs)

trunk/extsrc/gpcsw/gpcsrc/fits/burntool/persistfix.c

@@ -23 +23 @@
 
 /* Merge all overlapping streaks */
-   persist_merge(cell);
+  persist_merge(cell);
+
 
 /* Fix up all the persistence streaks */
    for(k=0; k<cell->npersist; k++) {
-
-/* Is this just a blasted area being carried for IPP? */
-      if((cell->persist)[k].func == BURN_BLASTED) {
-         if(cell->time - (cell->persist)[k].time > EXPIRE_TRAIL_TIME) {
+     /* Is this just a blasted area being carried for IPP? */
+     if((cell->persist)[k].func == BURN_BLASTED) {
+           if(cell->time - (cell->persist)[k].time > EXPIRE_TRAIL_TIME) {
             (cell->persist)[k].fiterr = FIT_EXPIRED;
-         }
-         continue;
-      }
+           }
+           continue;
+     }
 
 /* This had a significant positive slope, what do we do now? */
-      if( (cell->persist)[k].func == BURN_POSSLOPE) {
+     if( (cell->persist)[k].func == BURN_POSSLOPE) {
 /* try again, so let it slide through... */
-      }
+     }
 
 /* Fit the trail */
-      err = fit_trail(nx, ny, NX, imbuf, mbuf, cell->persist+k, 0, 
+     err = fit_trail(nx, ny, NX, imbuf, mbuf, cell->persist+k, 0, 
              cell->sky+cell->bias, cell->rms, BURN_EXP);
 
 /* Subtract out the fit */
-      if(!cell->persist[k].fiterr) 
-         err = sub_fit(nx, ny, NX, buf, cell->persist+k, 1);
+     if(!cell->persist[k].fiterr) 
+           err = sub_fit(nx, ny, NX, buf, cell->persist+k, 1);
    }
 
@@ -65 +65 @@
    for(k=0; k<cell->npersist; k++) {
       i = (cell->persist[k].ex - cell->persist[k].sx + 1) / 2;
-      ymid = (cell->persist[k].y0m + cell->persist[k].y1m +
-              cell->persist[k].y0p + cell->persist[k].y1p + 2) / 4;
-      printf("%3d %5d %3d %5d %5d %3d %5d %3d %5d %2d %8.5f %4d\n", 
+      /*ymid = (cell->persist[k].y0m + cell->persist[k].y1m +
+              cell->persist[k].y0p + cell->persist[k].y1p + 2) / 4;*/
+        ymid = (cell->persist[k].sy + cell->persist[k].ey + 1) / 2;
+      printf("%3d %5d %3d %5d %5d %3d %5d %3d %5d %2d %8.5f %8.5f\n", 
              k, cell->persist[k].cx, cell->persist[k].cy,
              cell->persist[k].max, 
@@ -73 +74 @@
              cell->persist[k].ex, cell->persist[k].ey, ymid,
              cell->persist[k].func, 
-             cell->persist[k].slope, NINT((cell->persist[k]).zero[i]));
+             cell->persist[k].slope, (cell->persist[k]).zero[i]);
    }
    return(0);