Changeset 39287 for trunk/Ohana/src/relastro/src/bcatalog.c

Timestamp:

Dec 17, 2015, 8:53:07 PM (11 years ago)

Author:

eugene

Message:

adding FitPosPMfixed_IRLS

File:

• trunk/Ohana/src/relastro/src/bcatalog.c (modified) (8 diffs)

trunk/Ohana/src/relastro/src/bcatalog.c

@@ -61 +61 @@
   }
 
+  int myNskip1 = 0, myNskip2 = 0, myNskip3 = 0, myNskip4 = 0, myNskip5 = 0, myNskip6 = 0;
+
   /* exclude stars not in range or with too few measurements */
   for (i = 0; i < catalog[0].Naverage; i++) {
     if (catalog[0].average[i].Nmeasure <= SRC_MEAS_TOOFEW) {
+      myNskip1 ++;
       continue;
     }
@@ -108 +111 @@
         REALLOCATE (subcatalog[0].secfilt, SecFilt, NAVERAGE*Nsecfilt);
       }
+      myNskip2 ++;
       continue;
     }
@@ -146 +150 @@
           Nskip1 ++;
         }
+        myNskip3 ++;
         continue;
       }
@@ -159 +164 @@
           Nskip2 ++;
         }
+        myNskip4 ++;
         continue;
       }
@@ -180 +186 @@
             free (date);
           }
+          myNskip5 ++;
           continue;
         }
@@ -220 +227 @@
     if (Nm <= SRC_MEAS_TOOFEW) {
       Nmeasure -= Nm;
+      myNskip6 ++;
       continue;
     }
@@ -230 +238 @@
     }
   }
+  fprintf (stderr, "skips: %d %d %d %d %d %d\n", myNskip1, myNskip2, myNskip3, myNskip4, myNskip5, myNskip6);
   REALLOCATE (subcatalog[0].average,  Average,     MAX (Naverage, 1));
   REALLOCATE (subcatalog[0].measureT, MeasureTiny, MAX (Nmeasure, 1));
@@ -244 +253 @@
     LimitDensityCatalog_ByNmeasureGrid (subcatalog, catalog);
   } else {
-    if (VERBOSE2) {
+    if (VERBOSE) {
       char *basename = filebasename (catalog[0].filename);
       fprintf (stderr, "subset of "OFF_T_FMT" ("OFF_T_FMT" total) stars, "OFF_T_FMT" ("OFF_T_FMT" total) measures for catalog %s\n",