Changeset 35769 for trunk/psphot/src/psphotKronIterate.c

Timestamp:

Jul 3, 2013, 2:43:13 PM (13 years ago)

Author:

eugene

Message:

thread psphotModelBackground, psphotAddOrSubNoise, psphotGuessModel; add recipe options to control selection of extended models fitted; iterate on the trail window a bit

Location:

trunk/psphot

Files:

• . (modified) (1 prop)
• src (modified) (1 prop)
• src/psphotKronIterate.c (modified) (1 diff)

trunk/psphot/src/psphotKronIterate.c

@@ -166 +166 @@
     }
 
+    fprintf (stderr, "--- starting KRON ---\n");
+
     // We measure the Kron Radius on a smoothed copy of the readout image
     psImage *smoothedImage = NULL;
     if (KRON_SMOOTH) {
+        psTimerStart ("psphot.kron.smooth");
+
         // Build the smoothed source image
         // Replace the subtracted sources
         psphotReplaceAllSourcesReadout(config, view, filerule, index, recipe, false);
+        psLogMsg ("psphot.kron", PS_LOG_INFO, "replaced %ld raw sources %f sec\n", sources->n, psTimerMark ("psphot.kron.smooth"));
+
+        // smoothedImage = psImageCopy(NULL, readout->image, PS_TYPE_F32);
+        // psImageSmooth(smoothedImage, KRON_SMOOTH_SIGMA, KRON_SMOOTH_NSIGMA);
+        
         // Copy the image and smooth
         psTimerStart ("psphot.kron.smooth");
-        smoothedImage = psImageCopy(NULL, readout->image, PS_TYPE_F32);
-        psImageSmooth(smoothedImage, KRON_SMOOTH_SIGMA, KRON_SMOOTH_NSIGMA);
+        bool oldThreads = psImageConvolveSetThreads(true); // Old value of threading in psImageConvolve
+        smoothedImage = psImageSmoothNoMask_Threaded (NULL, readout->image, KRON_SMOOTH_SIGMA, KRON_SMOOTH_NSIGMA, 0.25);
+        psImageConvolveSetThreads(oldThreads);
         psLogMsg ("psphot.kron", PS_LOG_INFO, "smoothed image %f sec\n", psTimerMark ("psphot.kron.smooth"));
 
         // remove the sources
+        psTimerStart ("psphot.kron.smooth");
         psphotRemoveAllSourcesReadout( config, view, filerule, index, recipe, false );
+        psLogMsg ("psphot.kron", PS_LOG_INFO, "removed %ld raw sources %f sec\n", sources->n, psTimerMark ("psphot.kron.smooth"));
+
         // Now subtract smooth versions of the sources from the smoothed image
         psTimerStart ("psphot.kron.smooth.sources");