Changeset 37575 for branches/eam_branches/ipp-20140904/Ohana/src/fakeastro

Timestamp:

Nov 10, 2014, 11:14:57 AM (12 years ago)

Author:

eugene

Message:

fix AstromOffsetTableFree, use 1st order only for fakeastro fitting

File:

• branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/fit_fake_stars.c (modified) (4 diffs)

branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/fit_fake_stars.c

@@ -33 +33 @@
     return FALSE;
   }
+  order_use = 1; // hard-wire linear solutions (nothing higher is needed for now)
 
   // I have Nstars with "reference" positions stars[].Rref,Dref and 
@@ -40 +41 @@
 
   int CHIPMAP = (image[0].coords.Npolyterms == -1);
-  image[0].coords.Npolyterms = CHIPMAP ? 1 : order_use;
 
-  CoordFit *fit = fit_init (image[0].coords.Npolyterms);
+  // XX1 image[0].coords.Npolyterms = CHIPMAP ? 1 : order_use;
+  CoordFit *fit = CHIPMAP ? fit_init (1) : fit_init (order_use);
 
   double *L, *M;
   ALLOCATE (L, double, Nstars);
   ALLOCATE (M, double, Nstars);
+
+  // static int Ntest = 0;
+  // char name[64];
+  // snprintf (name, 64, "test.%03d.dat", Ntest); Ntest ++;
+  // FILE *f = fopen (name, "w");
 
   // generate the fit matches
@@ -53 +59 @@
     RD_to_XY (&L[i], &M[i], stars[i].Rref, stars[i].Dref, image[0].coords.mosaic);
     fit_add (fit, stars[i].measure.Xccd, stars[i].measure.Yccd, L[i], M[i], 0.05);
+
+    // double Lo, Mo, Xo, Yo;
+    // RD_to_XY (&Lo, &Mo, stars[i].average.R, stars[i].average.D, image[0].coords.mosaic);
+    // RD_to_XY (&Xo, &Yo, stars[i].average.R, stars[i].average.D, &image[0].coords);
+    // fprintf (f, "%d : %f %f : %f %f : %f %f : %f %f : %f %f : %f %f\n", i, stars[i].Rref, stars[i].Dref, stars[i].average.R, stars[i].average.D, stars[i].measure.Xccd, stars[i].measure.Yccd, L[i], M[i], Lo, Mo, Xo, Yo);
   }
+  // fclose (f);
 
   // measure the fit, update the coords & object coordinates
@@ -66 +78 @@
   }
 
+  image[0].coords.Npolyterms = CHIPMAP ? 1 : order_use;
   if (!fit_apply_coords (fit, &image[0].coords, FALSE)) {
     fprintf (stderr, "failed to fit new model\n");