# include "relastro.h"
// we are sharing meas positions for all objects which (a) I own and (b) which have unowned detections

# define D_NMEASPOS 10000
int share_meas_pos (Catalog *catalog, int Ncatalog, RegionHostTable *regionHosts, int nloop) {

  int i;
  off_t j, k;

  off_t Nmeaspos = 0;
  off_t NMEASPOS = D_NMEASPOS;

  MeasPos *measpos = NULL;
  ALLOCATE (measpos, MeasPos, NMEASPOS);

  int myHost = regionHosts->index[REGION_HOST_ID];

  // XXX skip some catalogs based on UserPatch?
  for (i = 0; i < Ncatalog; i++) {
    for (j = 0; j < catalog[i].Naverage; j++) {

      // does this object have missing detections (does someone else need it?)
      if (catalog[i].nOwn_t[j] == catalog[i].average[j].Nmeasure) continue;

      off_t m = catalog[i].average[j].measureOffset;
      for (k = 0; k < catalog[i].average[j].Nmeasure; k++, m++) {
	if (!catalog[i].measureT[m].myDet) continue;

	set_meas_pos (&measpos[Nmeaspos], &catalog[i].average[j], &catalog[i].measureT[m]);
	Nmeaspos ++;
	CHECK_REALLOCATE (measpos, MeasPos, NMEASPOS, Nmeaspos, D_NMEASPOS);
      }
    }
  }

  // write out the measmag fits table AND write state in some file
  char *hostname = regionHosts->hosts[myHost].hostname;

  char *posfile = make_filename (CATDIR, hostname, REGION_HOST_ID, "measpos.fits");
  MeasPosSave (posfile, measpos, Nmeaspos);
  free (measpos);
  free (posfile);

  char *syncfile = make_filename (CATDIR, hostname, REGION_HOST_ID, "measpos.sync");
  update_sync_file (syncfile, nloop);
  free (syncfile);

  return TRUE;
}

int slurp_meas_pos (Catalog *catalog, int Ncatalog, RegionHostTable *regionHosts, int nloop) {

  off_t i, k;

  int Nmeaspos = 0;
  MeasPos *measpos = NULL;
  ALLOCATE (measpos, MeasPos, 1);

  fprintf (stderr, "grabbing meas object pos from other hosts...\n");

  for (i = 0; i < regionHosts->Nhosts; i++) {
    if (regionHosts->hosts[i].hostID == REGION_HOST_ID) continue;
    if (REGION_HOST_ID && !regionHosts->hosts[i].isNeighbor) continue;

    char *syncfile = make_filename (CATDIR, regionHosts->hosts[i].hostname, regionHosts->hosts[i].hostID, "measpos.sync");
    check_sync_file (syncfile, nloop);
    free (syncfile);
    
    off_t Nsubset = 0;
    char *posfile = make_filename (CATDIR, regionHosts->hosts[i].hostname, regionHosts->hosts[i].hostID, "measpos.fits");
    MeasPos *measposSubset = MeasPosLoad (posfile, &Nsubset);
    free (posfile);

    // merge_meas_pos reallocs measpos and frees the input measposSubset
    measpos = merge_meas_pos (measpos, &Nmeaspos, measposSubset, Nsubset);
  }

  for (i = 0; i < Nmeaspos; i++) {
    int objID = measpos[i].objID;
    int catID = measpos[i].catID;

    // set the meas mag
    int catSeq;
    off_t objSeq;
    if (!catID_and_objID_to_seq (catID, objID, &catSeq, &objSeq)) {
	// XXX what should I do if this does not match?
	continue;
    }
    off_t m = catalog[catSeq].average[objSeq].measureOffset;
    for (k = 0; k < catalog[catSeq].average[objSeq].Nmeasure; k++, m++) {
      if (catalog[catSeq].measureT[m].imageID != measpos[i].imageID) continue;
      catalog[catSeq].measureT[m].R = measpos[i].R;
      catalog[catSeq].measureT[m].D = measpos[i].D;
      if (catalog[catSeq].measure) {
	catalog[catSeq].measure[m].R = measpos[i].R;
	catalog[catSeq].measure[m].D = measpos[i].D;
      }
    }
  }
  free (measpos);

  fprintf (stderr, "DONE grabbing meas object pos from other hosts...\n");

  return TRUE;
}

int set_meas_pos (MeasPos *measpos, Average *average, MeasureTiny *measure) {

  measpos->R       = measure->R;
  measpos->D       = measure->D;
  measpos->objID   = average->objID;
  measpos->catID   = average->catID;
  measpos->imageID = measure->imageID;

  return TRUE;
}

MeasPos *merge_meas_pos (MeasPos *target, int *ntarget, MeasPos *source, int Nsource) {

  off_t i;

  REALLOCATE (target, MeasPos, *ntarget + Nsource);
  for (i = 0; i < Nsource; i++) {
    off_t n = i + *ntarget;
    target[n] = source[i];
  }
  
  free (source);

  *ntarget += Nsource;
  return (target);
}

