Changeset 31867

Timestamp:

Jul 12, 2011, 10:52:06 AM (15 years ago)

Author:

eugene

Message:

adding more tracing info on fits

File:

• branches/eam_branches/ipp-20110710/psphot/src/psphotExtendedSourceFits.c (modified) (4 diffs)

branches/eam_branches/ipp-20110710/psphot/src/psphotExtendedSourceFits.c

@@ -401 +401 @@
           // XXX psTraceSetLevel ("psphot.psphotModelWithPSF_LMM", 6);
 
+          // XXX this does not make sense in a threaded context
+          psTimerStart ("psphot.extended.fit");
+
           // fit the model as convolved or not
           pmModel *modelFit = NULL;
@@ -410 +413 @@
                   continue;
               }
-              psTrace ("psphot", 4, "fit psf-conv model for %f, %f : %s chisq = %f\n", source->moments->Mx, source->moments->My, pmModelClassGetName (modelFit->type), modelFit->chisq);
+              psTrace ("psphot", 4, "fit psf-conv model for %f, %f : %s chisq = %f (npix: %d, niter: %d)\n", 
+                       source->moments->Mx, source->moments->My, pmModelClassGetName (modelFit->type), modelFit->chisq, modelFit->nPix, modelFit->nIter);
               Nconvolve ++;
               if (!(modelFit->flags & (PM_MODEL_STATUS_BADARGS | PM_MODEL_STATUS_NONCONVERGE | PM_MODEL_STATUS_OFFIMAGE))) {
@@ -426 +430 @@
               }
               pmSourceCacheModel (source, maskVal);
-              psTrace ("psphot", 4, "fit plain model for %f, %f : %s chisq = %f\n", source->moments->Mx, source->moments->My, pmModelClassGetName (modelFit->type), modelFit->chisq);
+              psTrace ("psphot", 4, "fit plain model for %f, %f : %s chisq = %f (npix: %d, niter: %d)\n", 
+                       source->moments->Mx, source->moments->My, pmModelClassGetName (modelFit->type), modelFit->chisq, modelFit->nPix, modelFit->nIter);
               Nplain ++;
               if (!(modelFit->flags & (PM_MODEL_STATUS_BADARGS | PM_MODEL_STATUS_NONCONVERGE | PM_MODEL_STATUS_OFFIMAGE))) {
@@ -433 +438 @@
               }
           }
+
+          psLogMsg ("psphot", PS_LOG_INFO, "model fit: %7.5f sec, %5d npix, %2d iter, %s type\n", psTimerMark ("psphot.extended.fit"), modelFit->nPix, modelFit->nIter, pmModelClassGetName (modelFit->type));
 
           // save each of the model flux images and store the best