Changeset 13567 for trunk/psModules/src/objects/models/pmModel_QGAUSS.c

Timestamp:

May 31, 2007, 2:30:19 AM (19 years ago)

Author:

rhl

Message:

Set PM_SOURCE_MODE_BADPSF when we fail to get a good PSF estimate for an object

File:

• trunk/psModules/src/objects/models/pmModel_QGAUSS.c (modified) (5 diffs)

trunk/psModules/src/objects/models/pmModel_QGAUSS.c

@@ -88 +88 @@
         f2 = 1.0 / PS_SQR(params[PM_PAR_SYY]) - 1.0 / PS_SQR(params[PM_PAR_SXX]);
         q1 = PS_SQR(f1)*AR_RATIO - PS_SQR(f2);
-        assert (q1 > 0);
+        if (q1 <= 0) {
+            psError(PM_ERR_PSF, true, "Q1 < 0 (%g)", q1);
+            return false;
+        }
         q2  = 0.5*sqrt (q1);
     }
@@ -124 +127 @@
         if (fabs(beta[nParam]) > fabs(beta_lim)) {
             beta[nParam] = (beta[nParam] > 0) ? fabs(beta_lim) : -fabs(beta_lim);
+            psError(PM_ERR_PSF, true, "|beta[nParam==%d]| > |beta_lim|; %g v. %g",
+                    nParam, beta[nParam], beta_lim);
             return false;
         }
@@ -158 +163 @@
         if (params[nParam] < params_min) {
             params[nParam] = params_min;
+            psError(PM_ERR_PSF, true, "params[nParam==%d] < params_min; %g v. %g",
+                    nParam, params[nParam], params_min);
             return false;
         }
@@ -192 +199 @@
         if (params[nParam] > params_max) {
             params[nParam] = params_max;
+            psError(PM_ERR_PSF, true, "params[nParam==%d] > params_max; %g v. %g",
+                    nParam, params[nParam], params_max);
             return false;
         }
@@ -362 +371 @@
         if (psf->params_NEW->data[i] == NULL)
             continue;
-        PM_MODEL_LIMITS (PS_MINIMIZE_PARAM_MIN, i, out, NULL);
-        PM_MODEL_LIMITS (PS_MINIMIZE_PARAM_MAX, i, out, NULL);
+
+        if (!PM_MODEL_LIMITS(PS_MINIMIZE_PARAM_MIN, i, out, NULL) ||
+            !PM_MODEL_LIMITS(PS_MINIMIZE_PARAM_MAX, i, out, NULL)) {
+            psError(PM_ERR_PSF, false, "Setting limits at (r,c) = (%.1f, %.1f)",
+                    in[PM_PAR_XPOS], in[PM_PAR_YPOS]);
+            return false;
+        }
     }