# include "getstar.h"

/* add selected catalog objects to the output catalog */
int select_by_region (Catalog *output, Catalog *catalog, SkyRegion *region, int start, int end) {

  int i, j, n, Nm, offset, Nsecfilt;
  int Nave, NAVE, Nmeas, NMEAS;
  double R, D;

  Nsecfilt = output[0].Nsecfilt;

  if (output == NULL) Shutdown ("output not defined");

  /* identify selection criteria */
  if (VERBOSE) fprintf (stderr, "extracting for time range %d to %d\n", start, end);
  if (VERBOSE) fprintf (stderr, "extracting for region %f,%f to %f,%f\n", 
			region[0].Rmin, region[0].Dmin, region[0].Rmax, region[0].Dmax);

  /* select the average objects in this region */
  Nave = output[0].Naverage;
  NAVE = output[0].Naverage + 1000;
  REALLOCATE (output[0].average, Average, NAVE);
  REALLOCATE (output[0].secfilt, SecFilt, NAVE*Nsecfilt);

  Nmeas = output[0].Nmeasure;
  NMEAS = output[0].Nmeasure + 1000;
  REALLOCATE (output[0].measure, Measure, NMEAS);

  for (i = 0; i < catalog[0].Naverage; i++) {
    n = catalog[0].measure[i].averef;

    R = catalog[0].average[i].R;
    D = catalog[0].average[i].D;
    
    if (R < region[0].Rmin) continue;
    if (R > region[0].Rmax) continue;
    if (D < region[0].Dmin) continue;
    if (D > region[0].Dmax) continue;

    if (MagLimitUse && (catalog[0].average[i].M > MagLimitValue)) continue;

    output[0].average[Nave] = catalog[0].average[i];
    output[0].average[Nave].offset = Nmeas;
    for (j = 0; j < Nsecfilt; j++) {
      output[0].secfilt[Nsecfilt*Nave + j] = catalog[0].secfilt[Nsecfilt*i + j];
    }

    Nm = 0;
    offset = catalog[0].average[i].offset;

    for (j = 0; j < catalog[0].average[i].Nm; j++) {
      output[0].measure[Nmeas] = catalog[0].measure[offset + j];
      output[0].measure[Nmeas].averef = Nave;
      Nmeas ++;
      
      CHECK_REALLOCATE (output[0].measure, Measure, NMEAS, Nmeas, 1000);
    }      

    Nave ++;
    if (Nave == NAVE) {
      NAVE += 1000;
      REALLOCATE (output[0].average, Average, NAVE);
      REALLOCATE (output[0].secfilt, SecFilt, NAVE*Nsecfilt);
    }
  }
  output[0].Naverage = Nave;
  output[0].Nmeasure = Nmeas;
  fprintf (stderr, "output catalog has %d stars (%d measures, %d secfilt)\n",
	   output[0].Naverage, output[0].Nmeasure, output[0].Nsecfilt);
  return (TRUE);
}
