Changeset 33479

Timestamp:

Mar 11, 2012, 2:54:09 PM (14 years ago)

Author:

eugene

Message:

add testobj1 and testobj2

File:

• branches/eam_branches/ipp-20111122/Ohana/src/relastro/src/high_speed_objects.c (modified) (5 diffs)

branches/eam_branches/ipp-20111122/Ohana/src/relastro/src/high_speed_objects.c

@@ -10 +10 @@
   off_t l, i1, NAVERAGE, NMEASURE;
   int *slowMoving, *groupA, *groupB, status, Nmatch, Nmatchmeas, Nepoch, nv[2], Nmatchmeasobj;
-  int foundA, foundB;
+  int foundA, foundB, XVERB;
   double *X1, *Y1;
   double dX, dY, dR, RADIUS2;
@@ -76 +76 @@
   // mark (exclude) objects with both sets of target photcodes
   for (i = 0; i < catalog[0].Naverage; i++) {
+
+    XVERB  = (catalog[0].average[i].objID == OBJ_ID_SRC) && (catalog[0].average[i].catID == CAT_ID_SRC);
+    XVERB |= (catalog[0].average[i].objID == OBJ_ID_DST) && (catalog[0].average[i].catID == CAT_ID_DST);
+    if (XVERB) {
+      fprintf (stderr, "test object\n");
+    }
+
     // do any of the measures for this object match group A?
     m = catalog[0].average[i].measureOffset;
@@ -173 +180 @@
     nj = N1[j];
 
+    XVERB  = (catalog[0].average[ni].objID == OBJ_ID_SRC) && (catalog[0].average[ni].catID == CAT_ID_SRC);
+    XVERB |= (catalog[0].average[nj].objID == OBJ_ID_DST) && (catalog[0].average[nj].catID == CAT_ID_DST);
+    if (XVERB) {
+      fprintf (stderr, "test object %d or %d\n", (int) ni, (int) nj);
+    }
+
     if (slowMoving[ni]) NEXT_I;
     if (slowMoving[nj]) NEXT_J;
@@ -185 +198 @@
     // look for pairs that are within the maximum separation
     dX = X1[i] - X1[j];
+    if (XVERB) {
+      fprintf (stderr, "%d %d : %f\n", (int) ni, (int) nj, dX);
+    }
 
     if (dX <= -1.02*RADIUS) NEXT_I; // negative dX: i is too small
@@ -197 +213 @@
       nj = N1[J];
       if (!groupB[nj]) continue;
+
+      XVERB  = (catalog[0].average[ni].objID == OBJ_ID_SRC) && (catalog[0].average[ni].catID == CAT_ID_SRC);
+      XVERB |= (catalog[0].average[nj].objID == OBJ_ID_DST) && (catalog[0].average[nj].catID == CAT_ID_DST);
+      if (XVERB) {
+        fprintf (stderr, "test object pt2 %d or %d : %f\n", (int) ni, (int) nj, dX);
+        fprintf (stderr, ".\n");
+      }
 
       dY = Y1[i] - Y1[J];