Changeset 6949 for trunk/psphot/src/psphotEnsemblePSF.c

Timestamp:

Apr 21, 2006, 11:50:22 AM (20 years ago)

Author:

eugene

Message:

added tests for NaN results

File:

• trunk/psphot/src/psphotEnsemblePSF.c (modified) (3 diffs)

trunk/psphot/src/psphotEnsemblePSF.c

@@ -29 +29 @@
     psArray  *models = psArrayAlloc (sources->n);
     psVector *index  = psVectorAlloc (sources->n, PS_TYPE_U32);
-    models->n = index->n = 0;
+    // DROP models->n = index->n = 0;
 
     // option to limit analysis to a specific region
@@ -170 +170 @@
 
     psSparseConstraint constraint;
-    constraint.paramMin   = 0;
+    constraint.paramMin   = 0.0;
     constraint.paramMax   = 1e8;
     constraint.paramDelta = 1e8;
@@ -191 +191 @@
 
         // assign linearly-fitted normalization
+        if (isnan(norm->data.F32[i])) {
+            psAbort ("psphot", "ensemble source is nan");
+        }
         Fi->modelPSF->params->data.F32[1] = norm->data.F32[i];
         Fi->modelPSF->dparams->data.F32[1] = sqrt(sqrt(2) * norm->data.F32[i] / (sparse->Bfj->data.F32[i] * weight->data.F32[i]));