Changeset 20453
- Timestamp:
- Oct 28, 2008, 3:05:05 PM (18 years ago)
- Location:
- trunk
- Files:
-
- 1 deleted
- 24 edited
-
psModules/src/objects/pmPSFtry.c (modified) (5 diffs)
-
psphot/src/psphotAddNoise.c (modified) (2 diffs)
-
psphot/src/psphotApResid.c (modified) (2 diffs)
-
psphot/src/psphotBasicDeblend.c (modified) (2 diffs)
-
psphot/src/psphotBlendFit.c (modified) (2 diffs)
-
psphot/src/psphotChoosePSF.c (modified) (2 diffs)
-
psphot/src/psphotDeblendSatstars.c (modified) (2 diffs)
-
psphot/src/psphotFindDetections.c (modified) (1 diff)
-
psphot/src/psphotFindFootprints.c (modified) (4 diffs)
-
psphot/src/psphotFindPeaks.c (modified) (2 diffs)
-
psphot/src/psphotFitSourcesLinear.c (modified) (11 diffs)
-
psphot/src/psphotGuessModels.c (modified) (2 diffs)
-
psphot/src/psphotImageMedian.c (deleted)
-
psphot/src/psphotMagnitudes.c (modified) (2 diffs)
-
psphot/src/psphotMakeFluxScale.c (modified) (2 diffs)
-
psphot/src/psphotMakeGrowthCurve.c (modified) (2 diffs)
-
psphot/src/psphotMakeResiduals.c (modified) (2 diffs)
-
psphot/src/psphotModelBackground.c (modified) (2 diffs)
-
psphot/src/psphotReplaceUnfit.c (modified) (6 diffs)
-
psphot/src/psphotRoughClass.c (modified) (2 diffs)
-
psphot/src/psphotSignificanceImage.c (modified) (4 diffs)
-
psphot/src/psphotSkyReplace.c (modified) (2 diffs)
-
psphot/src/psphotSourceSize.c (modified) (2 diffs)
-
psphot/src/psphotSourceStats.c (modified) (2 diffs)
-
psphot/src/psphotSubtractBackground.c (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psModules/src/objects/pmPSFtry.c
r20250 r20453 5 5 * @author EAM, IfA 6 6 * 7 * @version $Revision: 1.6 4$ $Name: not supported by cvs2svn $8 * @date $Date: 2008-10- 18 00:49:19$7 * @version $Revision: 1.65 $ $Name: not supported by cvs2svn $ 8 * @date $Date: 2008-10-29 01:04:27 $ 9 9 * 10 10 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 168 168 169 169 // stage 1: fit an EXT model to all candidates PSF sources -- this is independent of the modeled 2D variations in the PSF 170 psTimerStart (" fit");170 psTimerStart ("psf.fit"); 171 171 for (int i = 0; i < psfTry->sources->n; i++) { 172 172 … … 196 196 Next ++; 197 197 } 198 psLogMsg ("psphot.psftry", 4, "fit ext: %f sec for %d of %ld sources\n", psTimerMark ("fit"), Next, sources->n);198 psLogMsg ("psphot.psftry", PS_LOG_MINUTIA, "fit ext: %f sec for %d of %ld sources\n", psTimerMark ("psf.fit"), Next, sources->n); 199 199 psTrace ("psModules.object", 3, "keeping %d of %ld PSF candidates (EXT)\n", Next, sources->n); 200 200 … … 213 213 214 214 // stage 3: refit with fixed shape parameters 215 psTimerStart (" fit");215 psTimerStart ("psf.fit"); 216 216 for (int i = 0; i < psfTry->sources->n; i++) { 217 217 … … 267 267 psfTry->psf->nPSFstars = Npsf; 268 268 269 psLogMsg ("psphot.psftry", 4, "fit psf: %f sec for %d of %ld sources\n", psTimerMark ("fit"), Npsf, sources->n);269 psLogMsg ("psphot.psftry", PS_LOG_MINUTIA, "fit psf: %f sec for %d of %ld sources\n", psTimerMark ("psf.fit"), Npsf, sources->n); 270 270 psTrace ("psModules.object", 3, "keeping %d of %ld PSF candidates (PSF)\n", Npsf, sources->n); 271 271 -
trunk/psphot/src/psphotAddNoise.c
r18001 r20453 20 20 PS_ASSERT (readout->parent->concepts, false); 21 21 22 psTimerStart ("psphot ");22 psTimerStart ("psphot.noise"); 23 23 24 24 // user-defined masks to test for good/bad pixels (build from recipe list if not yet set) … … 93 93 } 94 94 if (add) { 95 psLogMsg ("psphot.noise", PS_LOG_INFO, "add noise for %ld objects: %f sec\n", sources->n, psTimerMark ("psphot "));95 psLogMsg ("psphot.noise", PS_LOG_INFO, "add noise for %ld objects: %f sec\n", sources->n, psTimerMark ("psphot.noise")); 96 96 } else { 97 psLogMsg ("psphot.noise", PS_LOG_INFO, "sub noise for %ld objects: %f sec\n", sources->n, psTimerMark ("psphot "));97 psLogMsg ("psphot.noise", PS_LOG_INFO, "sub noise for %ld objects: %f sec\n", sources->n, psTimerMark ("psphot.noise")); 98 98 } 99 99 return true; -
trunk/psphot/src/psphotApResid.c
r20449 r20453 18 18 PS_ASSERT_PTR_NON_NULL(recipe, false); 19 19 20 psTimerStart ("psphot ");20 psTimerStart ("psphot.apresid"); 21 21 22 22 bool measureAptrend = psMetadataLookupBool (&status, recipe, "MEASURE.APTREND"); … … 218 218 psMetadataAdd (recipe, PS_LIST_TAIL, "NAPMIFIT", PS_DATA_S32 | PS_META_REPLACE, "number of apresid stars", psf->nApResid); 219 219 220 // psLogMsg ("psphot.apresid", PS_LOG_INFO, "measure full-frame aperture residuals for %d of %d objects: %f sec\n", Nkeep, Npsf, psTimerMark ("psphot"));221 220 psLogMsg ("psphot.apresid", PS_LOG_DETAIL, "aperture residual: %f +/- %f\n", psf->ApResid, psf->dApResid); 221 psLogMsg ("psphot.apresid", PS_LOG_INFO, "measure full-frame aperture residuals for: %f sec\n", psTimerMark ("psphot.apresid")); 222 222 223 223 psFree (mag); -
trunk/psphot/src/psphotBasicDeblend.c
r13006 r20453 10 10 int Nblend = 0; 11 11 12 psTimerStart ("psphot ");12 psTimerStart ("psphot.deblend.basic"); 13 13 14 14 float FRACTION = psMetadataLookupF32 (&status, recipe, "DEBLEND_PEAK_FRACTION"); … … 122 122 psFree (contour); 123 123 } 124 psLogMsg ("psphot.deblend", PS_LOG_INFO, "identified %d blended objects: %f sec\n", Nblend, psTimerMark ("psphot "));124 psLogMsg ("psphot.deblend", PS_LOG_INFO, "identified %d blended objects: %f sec\n", Nblend, psTimerMark ("psphot.deblend.basic")); 125 125 126 126 psFree (SN); -
trunk/psphot/src/psphotBlendFit.c
r19910 r20453 10 10 bool status; 11 11 12 psTimerStart ("psphot ");12 psTimerStart ("psphot.fit.nonlinear"); 13 13 14 14 // bit-masks to test for good/bad pixels … … 119 119 } 120 120 121 psLogMsg ("psphot.psphotBlendFit", PS_LOG_INFO, "fit models: %f sec for %d objects (%d psf, %d ext, %d failed, %ld skipped)\n", psTimerMark ("psphot "), Nfit, Npsf, Next, Nfail, sources->n - Nfit);121 psLogMsg ("psphot.psphotBlendFit", PS_LOG_INFO, "fit models: %f sec for %d objects (%d psf, %d ext, %d failed, %ld skipped)\n", psTimerMark ("psphot.fit.nonlinear"), Nfit, Npsf, Next, Nfail, sources->n - Nfit); 122 122 123 123 psphotVisualShowResidualImage (readout); -
trunk/psphot/src/psphotChoosePSF.c
r20287 r20453 21 21 bool status; 22 22 23 psTimerStart ("psphot ");23 psTimerStart ("psphot.choose.psf"); 24 24 25 25 // bit-masks to test for good/bad pixels … … 308 308 309 309 char *modelName = pmModelClassGetName (psf->type); 310 psLogMsg ("psphot.pspsf", PS_LOG_INFO, "select psf model: %f sec\n", psTimerMark ("psphot "));310 psLogMsg ("psphot.pspsf", PS_LOG_INFO, "select psf model: %f sec\n", psTimerMark ("psphot.choose.psf")); 311 311 psLogMsg ("psphot.pspsf", PS_LOG_INFO, "psf model %s, ApResid: %f +/- %f\n", modelName, psf->ApResid, psf->dApResid); 312 312 -
trunk/psphot/src/psphotDeblendSatstars.c
r13006 r20453 6 6 pmSource *source; 7 7 8 psTimerStart ("psphot ");8 psTimerStart ("psphot.deblend.sat"); 9 9 10 10 int Nblend = 0; … … 155 155 psFree (index); 156 156 157 psLogMsg ("psphot", PS_LOG_INFO, "found %d satstar blend peaks, leaving %ld sources: %f sec\n", Nblend, sources->n, psTimerMark ("psphot "));157 psLogMsg ("psphot", PS_LOG_INFO, "found %d satstar blend peaks, leaving %ld sources: %f sec\n", Nblend, sources->n, psTimerMark ("psphot.deblend.sat")); 158 158 return true; 159 159 } -
trunk/psphot/src/psphotFindDetections.c
r20177 r20453 8 8 float NSIGMA_PEAK = 25.0; 9 9 int NMAX = 0; 10 11 psTimerStart ("psphot");12 10 13 11 // user-defined masks to test for good/bad pixels (build from recipe list if not yet set) -
trunk/psphot/src/psphotFindFootprints.c
r19869 r20453 5 5 bool status; 6 6 7 psTimerStart (" footprints");7 psTimerStart ("psphot.footprints"); 8 8 9 9 int npixMin = psMetadataLookupS32(&status, recipe, "FOOTPRINT_NPIXMIN"); … … 40 40 psFree (detections->peaks); 41 41 42 psLogMsg ("psphot", PS_LOG_ DETAIL, "found %ld footprints: %f sec\n", footprints->n, psTimerMark ("footprints"));42 psLogMsg ("psphot", PS_LOG_MINUTIA, "found %ld footprints: %f sec\n", footprints->n, psTimerMark ("psphot.footprints")); 43 43 44 44 // optionally grow footprints isotropically by growRadius pixels 45 45 if (growRadius > 0) { 46 46 psArray *tmp = pmFootprintArrayGrow(footprints, growRadius); 47 psLogMsg ("psphot", PS_LOG_ DETAIL, "grow footprint coverage by %d pixels, %ld footprints become %ld footprints: %f sec\n", growRadius, footprints->n, tmp->n, psTimerMark ("footprints"));47 psLogMsg ("psphot", PS_LOG_MINUTIA, "grow footprint coverage by %d pixels, %ld footprints become %ld footprints: %f sec\n", growRadius, footprints->n, tmp->n, psTimerMark ("psphot.footprints")); 48 48 psFree(footprints); 49 49 footprints = tmp; … … 57 57 58 58 psArray *mergedFootprints = pmFootprintArraysMerge(detections->footprints, footprints, includePeaks); 59 psLogMsg ("psphot", PS_LOG_ DETAIL, "merged %ld new footprints with %ld existing ones: %f sec\n", footprints->n, detections->footprints->n, psTimerMark ("footprints"));59 psLogMsg ("psphot", PS_LOG_MINUTIA, "merged %ld new footprints with %ld existing ones: %f sec\n", footprints->n, detections->footprints->n, psTimerMark ("psphot.footprints")); 60 60 61 61 psFree(footprints); … … 68 68 psphotCullPeaks(readout->image, readout->weight, recipe, detections->footprints); 69 69 detections->peaks = pmFootprintArrayToPeaks(detections->footprints); 70 psLogMsg ("psphot", PS_LOG_INFO, "%ld peaks, %ld total footprints: %f sec\n", detections->peaks->n, detections->footprints->n, psTimerMark (" footprints"));70 psLogMsg ("psphot", PS_LOG_INFO, "%ld peaks, %ld total footprints: %f sec\n", detections->peaks->n, detections->footprints->n, psTimerMark ("psphot.footprints")); 71 71 72 72 return detections; -
trunk/psphot/src/psphotFindPeaks.c
r17870 r20453 8 8 bool status = false; 9 9 10 psTimerStart ("p eaks");10 psTimerStart ("psphot.peaks"); 11 11 12 12 // find the peaks in the smoothed image … … 48 48 pmPeaksWriteText (peaks, output); 49 49 } 50 psLogMsg ("psphot", PS_LOG_INFO, "%ld peaks: %f sec\n", peaks->n, psTimerMark ("p eaks"));50 psLogMsg ("psphot", PS_LOG_INFO, "%ld peaks: %f sec\n", peaks->n, psTimerMark ("psphot.peaks")); 51 51 52 52 return peaks; -
trunk/psphot/src/psphotFitSourcesLinear.c
r19910 r20453 20 20 // float r; 21 21 22 psTimerStart ("psphot ");22 psTimerStart ("psphot.linear"); 23 23 24 24 // bit-masks to test for good/bad pixels … … 93 93 psArrayAdd (fitSources, 100, source); 94 94 } 95 psLogMsg ("psphot.ensemble", PS_LOG_MINUTIA, "built fitSources: %f (%ld objects)\n", psTimerMark ("psphot"), sources->n);95 psLogMsg ("psphot.ensemble", PS_LOG_MINUTIA, "built fitSources: %f sec (%ld objects)\n", psTimerMark ("psphot.linear"), sources->n); 96 96 97 97 // vectors to store stats for each object … … 163 163 164 164 psSparseResort (sparse); 165 psLogMsg ("psphot.ensemble", PS_LOG_MINUTIA, "built matrix: %f (%d elements)\n", psTimerMark ("psphot"), sparse->Nelem);165 psLogMsg ("psphot.ensemble", PS_LOG_MINUTIA, "built matrix: %f sec (%d elements)\n", psTimerMark ("psphot.linear"), sparse->Nelem); 166 166 167 167 // set the sky, sky_x, sky_y components of border matrix 168 168 SetBorderMatrixElements (border, readout, fitSources, CONSTANT_PHOTOMETRIC_WEIGHTS, SKY_FIT_ORDER, markVal); 169 psLogMsg ("psphot.ensemble", PS_LOG_MINUTIA, "set border: %f (%d elements)\n", psTimerMark ("psphot"), sparse->Nelem);170 169 171 170 psSparseConstraint constraint; … … 184 183 skyfit = NULL; 185 184 } 186 psLogMsg ("psphot.ensemble", PS_LOG_MINUTIA, "solve matrix: %f (%d elements)\n", psTimerMark ("psphot"), sparse->Nelem);185 psLogMsg ("psphot.ensemble", PS_LOG_MINUTIA, "solve matrix: %f sec (%d elements)\n", psTimerMark ("psphot.linear"), sparse->Nelem); 187 186 188 187 // adjust I0 for fitSources and subtract … … 203 202 source->mode |= PM_SOURCE_MODE_SUBTRACTED; 204 203 } 205 psLogMsg ("psphot.ensemble", PS_LOG_MINUTIA, "sub models: %f (%d elements)\n", psTimerMark ("psphot"), sparse->Nelem);204 psLogMsg ("psphot.ensemble", PS_LOG_MINUTIA, "sub models: %f sec (%d elements)\n", psTimerMark ("psphot.linear"), sparse->Nelem); 206 205 207 206 // measure chisq for each source … … 211 210 pmSourceChisq (model, source->pixels, source->maskObj, source->weight, maskVal); 212 211 } 213 psLogMsg ("psphot.ensemble", PS_LOG_MINUTIA, "get chisqs: %f (%d elements)\n", psTimerMark ("psphot"), sparse->Nelem);212 psLogMsg ("psphot.ensemble", PS_LOG_MINUTIA, "get chisqs: %f sec (%d elements)\n", psTimerMark ("psphot.linear"), sparse->Nelem); 214 213 215 214 // psFree (index); … … 221 220 psFree (border); 222 221 223 psLogMsg ("psphot.ensemble", PS_LOG_INFO, "measure ensemble of PSFs: %f sec\n", psTimerMark ("psphot "));222 psLogMsg ("psphot.ensemble", PS_LOG_INFO, "measure ensemble of PSFs: %f sec\n", psTimerMark ("psphot.linear")); 224 223 225 224 psphotVisualShowResidualImage (readout); … … 241 240 fullArray = psRegionForImage (readout->mask, fullArray); 242 241 psImageMaskRegion (readout->mask, fullArray, "OR", markVal); 243 psLogMsg ("psphot.ensemble", PS_LOG_INFO, "step 1: %f sec\n", psTimerMark ("psphot"));244 242 245 243 // turn off MARK for all object pixels … … 252 250 psImageMaskCircle (source->maskView, x, y, model->radiusFit, "AND", PS_NOT_U8(markVal)); 253 251 } 254 psLogMsg ("psphot.ensemble", PS_LOG_INFO, "step 2: %f sec\n", psTimerMark ("psphot"));255 252 256 253 // accumulate the image statistics from the masked regions … … 289 286 } 290 287 } 291 psLogMsg ("psphot.ensemble", PS_LOG_INFO, "step 3: %f sec\n", psTimerMark ("psphot"));292 288 293 289 // turn off MARK for all image pixels 294 290 psImageMaskRegion (readout->mask, fullArray, "AND", PS_NOT_U8(markVal)); 295 psLogMsg ("psphot.ensemble", PS_LOG_INFO, "step 4: %f sec\n", psTimerMark ("psphot"));296 291 297 292 // set the Border T elements … … 310 305 psSparseBorderElementT (border, 2, 2, y2); 311 306 } 312 psLogMsg ("psphot.ensemble", PS_LOG_INFO, "step 5: %f sec\n", psTimerMark ("psphot"));313 307 314 308 return true; -
trunk/psphot/src/psphotGuessModels.c
r20411 r20453 27 27 bool status; 28 28 29 psTimerStart ("psphot ");29 psTimerStart ("psphot.models"); 30 30 31 31 // select the appropriate recipe information … … 142 142 } 143 143 144 psLogMsg ("psphot.models", 4, "built models for %ld objects: %f sec\n", sources->n, psTimerMark ("psphot "));144 psLogMsg ("psphot.models", 4, "built models for %ld objects: %f sec\n", sources->n, psTimerMark ("psphot.models")); 145 145 return true; 146 146 } -
trunk/psphot/src/psphotMagnitudes.c
r20069 r20453 6 6 int Nap = 0; 7 7 8 psTimerStart ("psphot ");8 psTimerStart ("psphot.mags"); 9 9 10 10 pmReadout *backModel = psphotSelectBackground (config, view); … … 60 60 } 61 61 62 psLogMsg ("psphot.magnitudes", PS_LOG_DETAIL, "measure magnitudes : %f sec for %ld objects (%d with apertures)\n", psTimerMark ("psphot "), sources->n, Nap);62 psLogMsg ("psphot.magnitudes", PS_LOG_DETAIL, "measure magnitudes : %f sec for %ld objects (%d with apertures)\n", psTimerMark ("psphot.mags"), sources->n, Nap); 63 63 return true; 64 64 } -
trunk/psphot/src/psphotMakeFluxScale.c
r20252 r20453 5 5 bool status = false; 6 6 7 psTimerStart (" residuals");7 psTimerStart ("psphot.fluxscale"); 8 8 9 9 int Nx = psMetadataLookupS32(&status, recipe, "PSF.FLUXSCALE.NX"); … … 72 72 psFree(trend); 73 73 74 psLogMsg ("psphot", PS_LOG_MINUTIA, "built flux scale: %f sec\n", psTimerMark ("psphot.fluxscale")); 75 74 76 return success; 75 77 } -
trunk/psphot/src/psphotMakeGrowthCurve.c
r20236 r20453 5 5 bool status; 6 6 7 psTimerStart (" growth");7 psTimerStart ("psphot.growth"); 8 8 9 9 // set limits on the aperture magnitudes … … 34 34 } 35 35 36 psLogMsg ("psphot", PS_LOG_ INFO, "built growth curve: %f sec\n", psTimerMark ("growth"));36 psLogMsg ("psphot", PS_LOG_MINUTIA, "built growth curve: %f sec\n", psTimerMark ("psphot.growth")); 37 37 38 38 return true; -
trunk/psphot/src/psphotMakeResiduals.c
r20451 r20453 7 7 psU8 mflux; 8 8 9 psTimerStart (" residuals");9 psTimerStart ("psphot.residuals"); 10 10 11 11 if (!psMetadataLookupBool(&status, recipe, "PSF.RESIDUALS")) return true; … … 273 273 psFree (B); 274 274 275 psLogMsg ("psphot.pspsf", PS_LOG_ INFO, "generate residuals for %ld objects: %f sec\n", input->n, psTimerMark ("residuals"));275 psLogMsg ("psphot.pspsf", PS_LOG_MINUTIA, "generate residuals for %ld objects: %f sec\n", input->n, psTimerMark ("psphot.residuals")); 276 276 277 277 psFree (xC); -
trunk/psphot/src/psphotModelBackground.c
r18833 r20453 9 9 static char *defaultStatsName = "FITTED_MEAN"; 10 10 11 psTimerStart ("psphot ");11 psTimerStart ("psphot.background"); 12 12 13 13 // find the currently selected readout … … 283 283 } 284 284 285 psLogMsg ("psphot", PS_LOG_ MINUTIA, "build median image: %f sec\n", psTimerMark ("psphot"));285 psLogMsg ("psphot", PS_LOG_INFO, "built background image: %f sec\n", psTimerMark ("psphot.background")); 286 286 287 287 psMetadataAddF32(recipe, PS_LIST_TAIL, "SKY_MEAN", PS_META_REPLACE, "sky mean", Value); -
trunk/psphot/src/psphotReplaceUnfit.c
r19869 r20453 6 6 pmSource *source; 7 7 8 psTimerStart ("psphot ");8 psTimerStart ("psphot.replace"); 9 9 10 10 for (int i = 0; i < sources->n; i++) { … … 19 19 source->mode &= ~PM_SOURCE_MODE_SUBTRACTED; 20 20 } 21 psLogMsg ("psphot.replace", 3, "replace unfitted models: %f sec (%ld objects)\n", psTimerMark ("psphot "), sources->n);21 psLogMsg ("psphot.replace", 3, "replace unfitted models: %f sec (%ld objects)\n", psTimerMark ("psphot.replace"), sources->n); 22 22 return true; 23 23 } … … 28 28 pmSource *source; 29 29 30 psTimerStart ("psphot ");30 psTimerStart ("psphot.replace"); 31 31 32 32 // user-defined masks to test for good/bad pixels (build from recipe list if not yet set) … … 43 43 source->mode &= ~PM_SOURCE_MODE_SUBTRACTED; 44 44 } 45 psLogMsg ("psphot.replace", PS_LOG_INFO, "replaced models for %ld objects: %f sec\n", sources->n, psTimerMark ("psphot "));45 psLogMsg ("psphot.replace", PS_LOG_INFO, "replaced models for %ld objects: %f sec\n", sources->n, psTimerMark ("psphot.replace")); 46 46 return true; 47 47 } … … 52 52 pmSource *source; 53 53 54 psTimerStart ("psphot ");54 psTimerStart ("psphot.replace"); 55 55 56 56 // user-defined masks to test for good/bad pixels (build from recipe list if not yet set) … … 67 67 source->mode |= PM_SOURCE_MODE_SUBTRACTED; 68 68 } 69 psLogMsg ("psphot.replace", PS_LOG_INFO, "replaced models for %ld objects: %f sec\n", sources->n, psTimerMark ("psphot "));69 psLogMsg ("psphot.replace", PS_LOG_INFO, "replaced models for %ld objects: %f sec\n", sources->n, psTimerMark ("psphot.replace")); 70 70 return true; 71 71 } -
trunk/psphot/src/psphotRoughClass.c
r19986 r20453 14 14 bool status; 15 15 16 psTimerStart ("psphot ");16 psTimerStart ("psphot.rough"); 17 17 18 18 // we make this measurement on a NxM grid of regions across the readout … … 41 41 psphotDumpMoments (recipe, sources); 42 42 43 psLogMsg ("psphot.roughclass", PS_LOG_INFO, "rough classification: %f sec\n", psTimerMark ("psphot "));43 psLogMsg ("psphot.roughclass", PS_LOG_INFO, "rough classification: %f sec\n", psTimerMark ("psphot.rough")); 44 44 45 45 psphotVisualPlotMoments (recipe, sources); -
trunk/psphot/src/psphotSignificanceImage.c
r20286 r20453 11 11 12 12 // smooth the image and weight map 13 psTimerStart (" smooth");13 psTimerStart ("psphot.smooth"); 14 14 15 15 // XXX we can a) choose fft to convolve if needed and b) multithread fftw … … 45 45 psImage *smooth_im = psImageCopy(NULL, readout->image, PS_TYPE_F32); 46 46 psImageSmoothMask(smooth_im, smooth_im, readout->mask, maskVal, SIGMA_SMTH, NSIGMA_SMTH, minGauss); 47 psLogMsg("psphot", PS_LOG_MINUTIA, "smooth image: %f sec\n", psTimerMark(" smooth"));47 psLogMsg("psphot", PS_LOG_MINUTIA, "smooth image: %f sec\n", psTimerMark("psphot.smooth")); 48 48 49 49 // Smooth the weight, applying the mask as we go. The variance is smoothed by the PSF^2, … … 57 57 psImageSmoothMask(smooth_wt, smooth_wt, readout->mask, maskVal, SIGMA_SMTH * M_SQRT1_2, 58 58 NSIGMA_SMTH, minGauss); 59 psLogMsg("psphot", PS_LOG_MINUTIA, "smooth weight: %f sec\n", psTimerMark(" smooth"));59 psLogMsg("psphot", PS_LOG_MINUTIA, "smooth weight: %f sec\n", psTimerMark("psphot.smooth")); 60 60 61 61 psImage *mask = readout->mask; … … 116 116 } 117 117 } 118 psLogMsg ("psphot", PS_LOG_INFO, "built smoothed signficance image: %f sec\n", psTimerMark(" smooth"));118 psLogMsg ("psphot", PS_LOG_INFO, "built smoothed signficance image: %f sec\n", psTimerMark("psphot.smooth")); 119 119 120 120 // optionally save example images under trace -
trunk/psphot/src/psphotSkyReplace.c
r13411 r20453 5 5 bool psphotSkyReplace (pmConfig *config, const pmFPAview *view) { 6 6 7 psTimerStart ("psphot ");7 psTimerStart ("psphot.skyreplace"); 8 8 9 9 // find the currently selected readout … … 28 28 } 29 29 } 30 psLogMsg ("psphot.sky", PS_LOG_DETAIL, "replace background flux : %f sec\n", psTimerMark ("psphot "));30 psLogMsg ("psphot.sky", PS_LOG_DETAIL, "replace background flux : %f sec\n", psTimerMark ("psphot.skyreplace")); 31 31 return true; 32 32 } -
trunk/psphot/src/psphotSourceSize.c
r20255 r20453 15 15 bool status; 16 16 17 psTimerStart ("psphot ");17 psTimerStart ("psphot.size"); 18 18 19 19 // user-defined masks to test for good/bad pixels (build from recipe list if not yet set) … … 228 228 229 229 psLogMsg ("psphot.size", PS_LOG_INFO, "measure source sizes for %ld sources: %f sec\n", 230 sources->n - first, psTimerMark ("psphot "));230 sources->n - first, psTimerMark ("psphot.size")); 231 231 232 232 psphotVisualPlotSourceSize (sources); -
trunk/psphot/src/psphotSourceStats.c
r19910 r20453 7 7 float BIG_RADIUS; 8 8 9 psTimerStart ("psphot ");9 psTimerStart ("psphot.stats"); 10 10 11 11 // bit-masks to test for good/bad pixels … … 119 119 } 120 120 121 psLogMsg ("psphot", PS_LOG_INFO, "%ld sources, %d moments, %d failed: %f sec\n", sources->n, Nmoments, Nfail, psTimerMark ("psphot "));121 psLogMsg ("psphot", PS_LOG_INFO, "%ld sources, %d moments, %d failed: %f sec\n", sources->n, Nmoments, Nfail, psTimerMark ("psphot.stats")); 122 122 123 123 psphotVisualShowMoments (sources); -
trunk/psphot/src/psphotSubtractBackground.c
r18836 r20453 10 10 pmReadout *backSub = NULL; 11 11 12 psTimerStart ("psphot ");12 psTimerStart ("psphot.background"); 13 13 14 14 // find the currently selected readout … … 63 63 } 64 64 65 psLogMsg ("psphot", PS_LOG_MINUTIA, "build resampled image: %f sec\n", psTimerMark ("psphot "));65 psLogMsg ("psphot", PS_LOG_MINUTIA, "build resampled image: %f sec\n", psTimerMark ("psphot.background")); 66 66 67 67 // back-sub image pixels, from output background file (don't create if not requested) … … 105 105 psphotSaveImage (NULL, image, name); 106 106 } 107 psLogMsg ("psphot", PS_LOG_INFO, "subtracted background model: %f sec\n", psTimerMark ("psphot "));107 psLogMsg ("psphot", PS_LOG_INFO, "subtracted background model: %f sec\n", psTimerMark ("psphot.background")); 108 108 109 109 // the pmReadout selected in this function are all view on entries in config->files
Note:
See TracChangeset
for help on using the changeset viewer.
