Changeset 10855 for trunk/psastro/src/psastroAstromGuess.c

Timestamp:

Dec 29, 2006, 8:36:33 AM (20 years ago)

Author:

eugene

Message:

extensive work to catch errors and clear leaks; implement mosaics

File:

• trunk/psastro/src/psastroAstromGuess.c (modified) (5 diffs)

trunk/psastro/src/psastroAstromGuess.c

@@ -11 +11 @@
     bool newFPA = true;
     bool status = false;
-    double RAmin = NAN;
-    double RAmax = NAN;
+    double RAmin = FLT_MAX;
+    double RAmax = FLT_MIN;
+    double DECmin = FLT_MAX;
+    double DECmax = FLT_MIN;
+
     double RAminSky = NAN;
     double RAmaxSky = NAN;
-    double DECmin = NAN;
-    double DECmax = NAN;
 
     pmChip *chip = NULL;
     pmCell *cell = NULL;
     pmReadout *readout = NULL;
-
-    DECmin = DECmax = -90;
-    RAmin = RAmax = RAminSky = RAmaxSky = 0;
 
     // select the current recipe
@@ -63 +61 @@
             RAminSky = fpa->toSky->R - M_PI;
             RAmaxSky = fpa->toSky->R + M_PI;
-            RAmin = RAmax = fpa->toSky->R;
-            DECmin = DECmax = fpa->toSky->D;
         }
 
@@ -79 +75 @@
                 if (rawstars == NULL) { continue; }
 
+                FILE *f = fopen ("rawstars.dat", "w");
                 for (int i = 0; i < rawstars->n; i++) {
                     pmAstromObj *raw = rawstars->data[i];
@@ -85 +82 @@
                     psPlaneTransformApply (raw->TP, fpa->toTPA, raw->FP);
                     psDeproject (raw->sky, raw->TP, fpa->toSky);
+
+                    fprintf (f, "%d  %f %f  %f %f  %f %f  %f %f\n", i,
+                             raw->sky->r, raw->sky->d, 
+                             raw->TP->x, raw->TP->y, 
+                             raw->FP->x, raw->FP->y, 
+                             raw->chip->x, raw->chip->y);
 
                     if (i < 10) {
@@ -122 +125 @@
                     DECmax = PS_MAX (raw->sky->d, DECmax);
                 }
+                fclose (f);
             }
         }