Changeset 27610 for trunk/Ohana/src/opihi/dvo/avmatch.c

Timestamp:

Apr 5, 2010, 3:47:10 PM (16 years ago)

Author:

eugene

Message:

fixing avmatch

File:

• trunk/Ohana/src/opihi/dvo/avmatch.c (modified) (5 diffs)

trunk/Ohana/src/opihi/dvo/avmatch.c

@@ -83 +83 @@
     if ((vec[i] = SelectVector (name, ANYVECTOR, TRUE)) == NULL) goto escape;
     ResetVector (vec[i], fields[i].type, NPTS);
+    for (n = 0; n < vec[i][0].Nelements; n++) {
+      if (vec[i][0].type == OPIHI_FLT) {
+        vec[i][0].elements.Flt[n] = NAN;
+      } else {
+        vec[i][0].elements.Int[n] = 0; // or NAN_INT?
+      }
+    }
   }
   ALLOCATE (index, off_t, NPTS);
@@ -103 +110 @@
     }
     dvo_catalog_unlock (&catalog);  /// we can unlock here since this is read-only (do not block other access)
+    if (catalog.Naverage == 0) {
+      dvo_catalog_free (&catalog);
+      continue;
+    }
 
     find_matches_by_vectors (skylist[0].regions[i], &catalog, RAvec, DECvec, RADIUS, index);
@@ -111 +122 @@
 
       if (Ncat == -1) continue;
-      if (Ncat == -2) {
-        for (n = 0; n < Nfields; n++) {
-          if (vec[n][0].type == OPIHI_FLT) {
-            vec[n][0].elements.Flt[Npts] = NAN;
-          } else {
-            vec[n][0].elements.Int[Npts] = 0; // or NAN_INT?
-          }
-        }
-        continue;
-      }
+      if (Ncat == -2) continue;
+
       // XXX set a 'found' vector to double check we catch everything?
       m = catalog.average[Ncat].measureOffset;
@@ -148 +151 @@
   dbFreeFields (fields, Nfields);
   SkyListFree (skylist);
-  dvo_catalog_free (&catalog);
   return (TRUE);
 
@@ -156 +158 @@
   dbFreeFields (fields, Nfields);
   SkyListFree (skylist);
-  dvo_catalog_free (&catalog);
   return (FALSE);