Changeset 20453 for trunk/psphot/src/psphotFitSourcesLinear.c
- Timestamp:
- Oct 28, 2008, 3:05:05 PM (18 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotFitSourcesLinear.c (modified) (11 diffs)
Legend:
- Unmodified
- Added
- Removed
-
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;
Note:
See TracChangeset
for help on using the changeset viewer.
