Changeset 24034 for trunk/psModules/src/astrom/pmAstrometryObjects.c

Timestamp:

May 1, 2009, 5:11:04 PM (17 years ago)

Author:

eugene

Message:

prevent multiple matches to the same input or reference source

File:

• trunk/psModules/src/astrom/pmAstrometryObjects.c (modified) (5 diffs)

trunk/psModules/src/astrom/pmAstrometryObjects.c

@@ -83 +83 @@
 {
     psArray *matches = psArrayAllocEmpty(x1->n);
+    psVector *found1 = psVectorAlloc(x1->n, PS_TYPE_S8);
+    psVector *found2 = psVectorAlloc(x2->n, PS_TYPE_S8);
 
     const double RADIUS_SQR = PS_SQR(RADIUS);
     double dX, dY, dR;
+
+    psVectorInit (found1, 0);
+    psVectorInit (found2, 0);
 
     int jStart;
@@ -100 +105 @@
         }
 
+        if (found1->data.S8[i]) {
+            i++;
+            continue;
+        }
+        if (found2->data.S8[j]) {
+            j++;
+            continue;
+        }
+
         jStart = j;
-        while (fabs(dX) < RADIUS && j < x2->n) {
+        while ((fabs(dX) < RADIUS) && (j < x2->n)) {
 
             dX = x1->data.F64[i] - x2->data.F64[j];
@@ -111 +125 @@
                 continue;
             }
+            if (found2->data.S8[j]) {
+                j++;
+                continue;
+            }
 
             // got a match; add to output list
@@ -117 +135 @@
             psFree (match);
 
+            found1->data.S8[i] = 1;
+            found2->data.S8[j] = 1;
+
             j++;
         }
@@ -122 +143 @@
         i++;
     }
+    psFree (found1);
+    psFree (found2);
+
     return (matches);
 }