Changeset 30333

Timestamp:

Jan 21, 2011, 10:56:41 AM (15 years ago)

Author:

eugene

Message:

exclude NANs from chisq vs flux fit

File:

• branches/eam_branches/ipp-20101205/psModules/src/imcombine/pmSubtractionEquation.c (modified) (1 diff)

branches/eam_branches/ipp-20101205/psModules/src/imcombine/pmSubtractionEquation.c

@@ -1426 +1426 @@
     psVectorAppend(momentVector, moment);
     psVectorAppend(fluxesVector, flux);
-    psVectorAppend(stampMask, 0);
+
+    // check that the last appended values are ok:
+    int Nelem = fluxesVector->n - 1;
+    bool valid = true;
+    valid &= isfinite(chisqVector->data.F32[Nelem]);
+    valid &= isfinite(fluxesVector->data.F32[Nelem]);
+    valid &= isfinite(momentVector->data.F32[Nelem]);
+    if (valid) {
+      psVectorAppend(stampMask, 0);
+    } else {
+      psVectorAppend(stampMask, 0x02);
+    }
     return true;
 }