IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 20190


Ignore:
Timestamp:
Oct 16, 2008, 9:18:27 AM (18 years ago)
Author:
eugene
Message:

add some verbosity to the loading

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Ohana/src/relphot/src/load_catalogs.c

    r15743 r20190  
    33Catalog *load_catalogs (SkyList *skylist, int *Ncatalog) {
    44
    5   int i, Nstar;
     5  int i, Nmeas, Nstar, Nmeas_total, Nstar_total;
    66  Catalog *catalog, tcatalog;
    77
     
    99
    1010  ALLOCATE (catalog, Catalog, skylist[0].Nregions);
     11
     12  Nmeas_total = Nstar_total = 0;
    1113
    1214  // load data from each region file, only use bright stars
     
    2527    }
    2628    if (VERBOSE && !tcatalog.Naves_disk) fprintf (stderr, "no data in %s, skipping\n", tcatalog.filename);
     29    Nstar_total += tcatalog.Naverage;
     30    Nmeas_total += tcatalog.Nmeasure;
    2731
    2832    // select only the brighter stars
     
    3337
    3438  // XXX keep this test?
    35   Nstar = 0;
     39  Nstar = Nmeas = 0;
    3640  for (i = 0; i < skylist[0].Nregions; i++) {
    3741    Nstar += catalog[i].Naverage;
     42    Nmeas += catalog[i].Nmeasure;
    3843  }
    3944  if (Nstar < 2) {
    4045    fprintf (stderr, "warning: insufficient stars %d\n", Nstar);
    4146  }
     47
     48  fprintf (stderr, "using %d of %d stars (%d of %d measurements)\n", Nstar, Nstar_total, Nmeas, Nmeas_total);
    4249
    4350  // XXX consider only returning the populated catalogs
Note: See TracChangeset for help on using the changeset viewer.