Changeset 27416

Timestamp:

Mar 23, 2010, 3:18:02 PM (16 years ago)

Author:

Paul Price

Message:

Log getstar command.

File:

• branches/pap_stack/psastro/src/psastroLoadRefstars.c (modified) (6 diffs)

branches/pap_stack/psastro/src/psastroLoadRefstars.c

@@ -1 +1 @@
 /** @file psastroLoadRefstars.c
  *
- *  @brief 
+ *  @brief
  *
  *  @ingroup libpsastro
@@ -59 +59 @@
     psMetadata *catdir_folder = psMetadataLookupMetadata(&status, catdirs, catdir_recipe);
     if (catdir_folder) {
-        // randomly choose one of the entries
-        psLogMsg ("psastro", 3, "choosing catdir_folder\n");
-        int nEntry = catdir_folder->list->n;
-
-        psRandom *rng = psRandomAlloc(PS_RANDOM_TAUS);  
-        double frnd = psRandomUniform(rng);
-        int entry = PS_MIN(nEntry - 1, PS_MAX(0, nEntry * frnd));
-        psFree(rng);
-
-        psMetadataItem *item = psListGet(catdir_folder->list, entry);
-        if (item->type != PS_DATA_STRING) {
-            psError(PSASTRO_ERR_CONFIG, true, "Invalid entry in PSASTRO.CATDIR folder: %s\n", item->name);
-            return false;
-        }
-        catdir_virtual = item->data.str;
+        // randomly choose one of the entries
+        psLogMsg ("psastro", 3, "choosing catdir_folder\n");
+        int nEntry = catdir_folder->list->n;
+
+        psRandom *rng = psRandomAlloc(PS_RANDOM_TAUS);
+        double frnd = psRandomUniform(rng);
+        int entry = PS_MIN(nEntry - 1, PS_MAX(0, nEntry * frnd));
+        psFree(rng);
+
+        psMetadataItem *item = psListGet(catdir_folder->list, entry);
+        if (item->type != PS_DATA_STRING) {
+            psError(PSASTRO_ERR_CONFIG, true, "Invalid entry in PSASTRO.CATDIR folder: %s\n", item->name);
+            return false;
+        }
+        catdir_virtual = item->data.str;
     }
 
@@ -133 +133 @@
     psTrace ("psastro", 3, "%s\n", getstarCommand);
 
+    psLogMsg("psastro", PS_LOG_INFO, "getstar command: %s", getstarCommand);
+
     // XXX use psPipe: catch stderr, stdout, allow for Nsec timeout...
     // use fork to add timeout capability
@@ -160 +162 @@
       // assert (outfile);
       // for (int nn = 0; nn < refstars->n; nn++) {
-      //          pmAstromObj *ref = refstars->data[nn];
-      //          fprintf (outfile, "%lf %lf\n", ref->sky->r*PS_DEG_RAD, ref->sky->d*PS_DEG_RAD);
+      //          pmAstromObj *ref = refstars->data[nn];
+      //          fprintf (outfile, "%lf %lf\n", ref->sky->r*PS_DEG_RAD, ref->sky->d*PS_DEG_RAD);
       // }
       // fclose (outfile);
@@ -217 +219 @@
             ref->sky->d   = RAD_DEG*psMetadataLookupF32 (&status, row, "DEC");
             ref->Mag      = 0.001*psMetadataLookupS32 (&status, row, "MAG");  // ELIXIR uses millimags
-            ref->Color    = 0.0;
+            ref->Color    = 0.0;
         } else {
             ref->sky->r   = RAD_DEG*psMetadataLookupF64 (&status, row, "RA");
             ref->sky->d   = RAD_DEG*psMetadataLookupF64 (&status, row, "DEC");
             ref->Mag      = psMetadataLookupF32 (&status, row, "MAG"); // PANSTARRS uses mags
-            ref->Color    = 0.0;
+            ref->Color    = 0.0;
         }
 
@@ -259 +261 @@
         float MagC1   = psMetadataLookupF32 (&status, row, "MAG_C1");
         float MagC2   = psMetadataLookupF32 (&status, row, "MAG_C2");
-        if (!isnan(MagC1) && !isnan(MagC2)) {
-            ref->Color = MagC1 - MagC2;
-        } else {
-            // XXX save the color and the slope in the table header?
-            ref->Color = 0.0;
-        }
+        if (!isnan(MagC1) && !isnan(MagC2)) {
+            ref->Color = MagC1 - MagC2;
+        } else {
+            // XXX save the color and the slope in the table header?
+            ref->Color = 0.0;
+        }
 
         // XXX VERY temporary hack to avoid M31 bulge