Changeset 30784 for branches/eam_branches/ipp-20110213/psphot/src/psphotFitSourcesLinearStack.c

Timestamp:

Mar 3, 2011, 3:11:09 PM (15 years ago)

Author:

eugene

Message:

fix the calculation of chisq values for constant errors (re-calculated, do not attempt to correct raw value)

File:

• branches/eam_branches/ipp-20110213/psphot/src/psphotFitSourcesLinearStack.c (modified) (1 diff)

branches/eam_branches/ipp-20110213/psphot/src/psphotFitSourcesLinearStack.c

@@ -163 +163 @@
     for (int i = 0; i < fitSources->n; i++) {
         pmSource *source = fitSources->data[i];
-        if (source->mode & PM_SOURCE_MODE_NONLINEAR_FIT) continue;
         pmModel *model = pmSourceGetModel (NULL, source);
-        pmSourceChisq (model, source->pixels, source->maskObj, source->variance, maskVal, COVAR_FACTOR, 1);
+        if (!(source->mode & PM_SOURCE_MODE_NONLINEAR_FIT)) {
+            model->nPar = 1; // LINEAR-only sources have 1 parameter; NONLINEAR sources have their original value
+        }
+        pmSourceChisq (model, source->pixels, source->maskObj, source->variance, maskVal);
     }
     psLogMsg ("psphot.ensemble", PS_LOG_MINUTIA, "get chisqs: %f sec (%d elements)\n", psTimerMark ("psphot.linear"), sparse->Nelem);