Changeset 35639 for branches/eam_branches/ipp-20130509/psModules/src/objects/pmSourceFitModel.c

Timestamp:

Jun 9, 2013, 3:30:23 PM (13 years ago)

Author:

eugene

Message:

add threading for pmPSFtryFitEXT and pmPSFtryFitPSF; pass chisq conv options to pmSourceFitModel

File:

• branches/eam_branches/ipp-20130509/psModules/src/objects/pmSourceFitModel.c (modified) (3 diffs)

branches/eam_branches/ipp-20130509/psModules/src/objects/pmSourceFitModel.c

@@ -62 +62 @@
     opt->poissonErrors = true;
     opt->saveCovariance = false;
+
+    // we default to the old algorithm
+    opt->gainFactorMode = 0;
+    opt->chisqConvergence = true;
 
     return opt;
@@ -241 +245 @@
 
     psMinimization *myMin = psMinimizationAlloc (options->nIter, options->minTol, options->maxTol);
+    myMin->gainFactorMode = options->gainFactorMode;
+    myMin->chisqConvergence = options->chisqConvergence;
 
     psImage *covar = psImageAlloc (params->n, params->n, PS_TYPE_F32);
@@ -274 +280 @@
     model->flags |= PM_MODEL_STATUS_FITTED;
     if (!fitStatus) model->flags |= PM_MODEL_STATUS_NONCONVERGE;
-    if (myMin->lastDelta > myMin->minTol) model->flags |= PM_MODEL_STATUS_WEAK_FIT;
+
+    if (myMin->chisqConvergence) {
+      if (myMin->lastDelta > myMin->minTol) model->flags |= PM_MODEL_STATUS_WEAK_FIT;
+    } else {
+      if (myMin->rParSigma > myMin->minTol*nParams) model->flags |= PM_MODEL_STATUS_WEAK_FIT;
+    }
 
     // get the Gauss-Newton distance for fixed model parameters