Changeset 20943 for trunk/magic/remove/src/streaksremove.c

Timestamp:

Dec 8, 2008, 3:27:32 PM (17 years ago)

Author:

bills

Message:

add some timers

File:

• trunk/magic/remove/src/streaksremove.c (modified) (8 diffs)

trunk/magic/remove/src/streaksremove.c

@@ -16 +16 @@
 
     psLibInit(NULL);
+    psTimerStart("STREAKSREMOVE");
 
     pmConfig *config = parseArguments(argc, argv);
@@ -61 +62 @@
         // if no skycells are provided sfiles->exciseAll is set to true
 
+        psTimerStart("COMPUTE_WARPED_PIXELS");
         if (! computeWarpedPixels(sfiles) ) {
             // we have no choice to excise all pixels
             exciseAll = true;
         }
+        psF64 cwp_t = psTimerClear("COMPUTE_WARPED_PIXELS");
+        printf("time to compute warped pixels: %f\n", cwp_t);
     }
     
@@ -81 +85 @@
     do {
         bool exciseImageCube = false;
+
+        psTimerStart("PROCESS_COMPONENT");
 
         // read the images and copy the data from the inputs to the outputs
@@ -95 +101 @@
             // Identify pixels to mask because of streaks
 
+            psTimerStart("GET_STREAK_PIXELS");
+
             StreakPixels *pixels = getStreakPixels(sfiles, streaks);
 
-            // XXX: 
-            // 
-            // PFS: Need to get mask weight for PM_MASK_DETECTOR.  Is this stored
-            // in a config somewhere?  Not sure how to properly set the maskValue.
-
+            printf("time to get streak pixels: %f\n", psTimerClear("GET_STREAK_PIXELS"));
 
             if (checkNonWarpedPixels) {
+                psTimerStart("EXCISE_NON_WARPED");
+
                 exciseNonWarpedPixels(sfiles, maskStreak);
+
+                printf("time to excise non warped pixels: %f\n", psTimerClear("EXCISE_NON_WARPED"));
             }
             
             if (sfiles->inImage->image) {
+                psTimerStart("REMOVE_STREAKS");
                 for (int i = 0; i < psArrayLength (pixels); ++i) {
                     PixelPos *pixelPos = psArrayGet (pixels, i);
@@ -120 +129 @@
                     }
                 }
+                printf("time to remove streak pixels: %f\n", psTimerClear("REMOVE_STREAKS"));
             } else { 
                 // this component contains an image cube, excise it completely
@@ -130 +140 @@
         writeImages(sfiles, exciseImageCube);
 
+        printf("time to process component %d: %f\n", sfiles->extnum, psTimerClear("PROCESS_COMPONENT"));
     } while (streakFilesNextExtension(sfiles));
 
+
+    psTimerStart("CLOSE_IMAGES");
     // close all files
     closeImages(sfiles);
+    printf("time to close images: %f\n", psTimerClear("CLOSE_IMAGES"));
 
     // NOTE: from here on we can't just quit if something goes wrong.
@@ -175 +189 @@
     }
 #endif  // REPLACE, REMOVE
+    printf("time to run streaksremove: %f\n", psTimerClear("STREAKSREMOVE"));
+
     // nebServerFree(ourNebServer);
     psFree(config);
@@ -183 +199 @@
     //  PAU
     fprintf(stderr, "Found %d leaks at %s\n", psMemCheckLeaks (0, NULL, NULL, false), "streaksremove");
+
 
     return 0;