Changeset 42894 for trunk/psModules/src/objects/pmSourceMatch.c

Timestamp:

Jun 7, 2025, 2:27:08 PM (13 months ago)

Author:

eugene

Message:

filter out images with NAN zero points from psuedo-relphot analysis

File:

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

trunk/psModules/src/objects/pmSourceMatch.c

@@ -429 +429 @@
             numMeasurements++;
             int index = match->image->data.U32[j]; // Image index
+            if (badImage->data.U8[index]) { continue; }
+
             float mag = match->mag->data.F32[j]; // Measured magnitude
             double magErr2 = PS_SQR(match->magErr->data.F32[j]) + sysErr2; // Error in measured magnitude
@@ -459 +461 @@
             }
             int index = match->image->data.U32[j]; // Image index
+            if (badImage->data.U8[index]) { continue; }
+
             float mag = match->mag->data.F32[j]; // Measured magnitude
             double magErr2 = PS_SQR(match->magErr->data.F32[j]) + sysErr2; // Error in measured magnitude
@@ -487 +491 @@
             }
             int index = match->image->data.U32[j]; // Image index
-            if (badImage->data.U8[index]) {
-                continue;
-            }
+            if (badImage->data.U8[index]) { continue; }
+
             float mag = match->mag->data.F32[j]; // Measured magnitude
             float magErr2 = PS_SQR(match->magErr->data.F32[j]) + sysErr2; // Error in measured magnitude
@@ -607 +610 @@
             numMeasurements++;
             int index = match->image->data.U32[j]; // Image index
-            if (badImage->data.U8[index]) {
-                continue;
-            }
+            if (badImage->data.U8[index]) { continue; }
+
             float mag = match->mag->data.F32[j]; // Measured magnitude
             float magErr = match->magErr->data.F32[j]; // Error in measured magnitude
@@ -857 +859 @@
 
             int index = match->image->data.U32[j]; // Image index
+
             float dx = match->x->data.F32[j] - xShift->data.F32[index] - xStar->data.F32[i];
             float dy = match->y->data.F32[j] - yShift->data.F32[index] - yStar->data.F32[i];