Changeset 9576 for trunk/psphot/src/psphotSourceFits.c
- Timestamp:
- Oct 13, 2006, 5:13:31 PM (20 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotSourceFits.c (modified) (10 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotSourceFits.c
r9270 r9576 12 12 return status; 13 13 } 14 psTrace ("psphot .blend", 5, "trying blend...\n");14 psTrace ("psphot", 5, "trying blend...\n"); 15 15 16 16 // save the PSF model from the Ensemble fit … … 79 79 // if this one failed, skip it 80 80 if (!psphotEvalPSF (blend, model)) { 81 psTrace ("psphot .blend", 5, "failed on blend %d of %ld\n", i, modelSet->n);81 psTrace ("psphot", 5, "failed on blend %d of %ld\n", i, modelSet->n); 82 82 continue; 83 83 } … … 86 86 psFree(blend->modelPSF); 87 87 blend->modelPSF = psMemIncrRefCounter (model); 88 psTrace ("psphot .blend", 5, "fitted blend as PSF\n");88 psTrace ("psphot", 5, "fitted blend as PSF\n"); 89 89 pmModelSub (source->pixels, source->mask, model, false, false); 90 90 blend->mode |= PM_SOURCE_MODE_SUBTRACTED; … … 96 96 // evaluate the primary object 97 97 if (!psphotEvalPSF (source, PSF)) { 98 psTrace ("psphot .blend", 5, "failed on blend 0 of %ld\n", modelSet->n);98 psTrace ("psphot", 5, "failed on blend 0 of %ld\n", modelSet->n); 99 99 psFree (PSF); 100 100 return false; 101 101 } 102 102 103 psTrace ("psphot .blend", 5, "fitted primary as PSF\n");103 psTrace ("psphot", 5, "fitted primary as PSF\n"); 104 104 pmModelSub (source->pixels, source->mask, PSF, false, false); 105 105 psFree (source->modelPSF); … … 140 140 } 141 141 142 psTrace ("psphot .blend", 5, "fitted as PSF\n");142 psTrace ("psphot", 5, "fitted as PSF\n"); 143 143 pmModelSub (source->pixels, source->mask, PSF, false, false); 144 144 … … 189 189 if (!pmSourceMoments (source, EXT_MOMENTS_RAD)) return false; 190 190 191 psTrace ("psphot .blend", 5, "trying blob...\n");191 psTrace ("psphot", 5, "trying blob...\n"); 192 192 193 193 // this temporary source is used as a place-holder by the psphotEval functions below … … 219 219 220 220 if (okEXT && okDBL) { 221 psTrace ("psphot .blend", 5, "blob chisq: %f vs %f\n", chiEXT, chiDBL);221 psTrace ("psphot", 5, "blob chisq: %f vs %f\n", chiEXT, chiDBL); 222 222 // XXX EAM : a bogus bias: need to examine this better 223 223 if (3*chiEXT > chiDBL) goto keepDBL; … … 237 237 psFree (DBL); 238 238 pmModelSub (source->pixels, source->mask, EXT, false, false); 239 psTrace ("psphot .blend", 5, "blob as EXT: %f %f\n", EXT->params->data.F32[2], EXT->params->data.F32[3]);239 psTrace ("psphot", 5, "blob as EXT: %f %f\n", EXT->params->data.F32[2], EXT->params->data.F32[3]); 240 240 241 241 // save new model … … 250 250 pmModelSub (source->pixels, source->mask, (pmModel *) DBL->data[0], false, false); 251 251 pmModelSub (source->pixels, source->mask, (pmModel *) DBL->data[1], false, false); 252 psTrace ("psphot .blend", 5, "blob as DBL: %f %f\n", ONE->params->data.F32[2], ONE->params->data.F32[3]);252 psTrace ("psphot", 5, "blob as DBL: %f %f\n", ONE->params->data.F32[2], ONE->params->data.F32[3]); 253 253 254 254 // drop old model, save new second model... … … 339 339 340 340 if ((source->moments->Sx < 1e-3) || (source->moments->Sx < 1e-3)) { 341 psTrace ("psphot SourceFits", 5, "problem source: moments: %f %f\n", source->moments->Sx, source->moments->Sy);341 psTrace ("psphot", 5, "problem source: moments: %f %f\n", source->moments->Sx, source->moments->Sy); 342 342 } 343 343
Note:
See TracChangeset
for help on using the changeset viewer.
