Changeset 17561 for trunk/psphot/src/psphotIsophotal.c

Timestamp:

May 7, 2008, 11:14:09 AM (18 years ago)

Author:

Paul Price

Message:

Reworking due to compiler error which seems to be due to initialising static value to NAN: 'initializer element is not constant'. Got rid of the static --- no need for this, unless the psMetadataLookup is too heavy (shouldn't be, and there's more processing going on in these functions than is going on in psMetadataLookup).

File:

• trunk/psphot/src/psphotIsophotal.c (modified) (5 diffs)

trunk/psphot/src/psphotIsophotal.c

@@ -1 +1 @@
 # include "psphotInternal.h"
-
-static float ISOPHOT_FLUX = NAN;
 
 bool psphotIsophotal (pmSource *source, psMetadata *recipe, psMaskType maskVal) {
@@ -16 +14 @@
 
   // flux at which to measure isophotal parameters
-  if (!isfinite(ISOPHOT_FLUX)) {
-    // XXX ISOPHOTAL_FLUX should be specified in mags, need the zero point to get counts/sec
-    ISOPHOT_FLUX = psMetadataLookupF32 (&status, recipe, "ISOPHOTAL_FLUX");
-    assert (status);
-  }
+  // XXX ISOPHOTAL_FLUX should be specified in mags, need the zero point to get counts/sec
+  float ISOPHOT_FLUX = psMetadataLookupF32 (&status, recipe, "ISOPHOTAL_FLUX");
+  assert (status);
 
   // find the first bin below the flux level and the last above the level
   // XXX can this be done faster with bisection?
-  // XXX do I need to worry about crazy outliers?  
+  // XXX do I need to worry about crazy outliers?
   // XXX should i be smoothing or fitting the curve?
   int firstBelow = -1;
@@ -44 +40 @@
     return false;
   }
-  
+
   // need to examine pixels in this vicinity
   float isophotalFluxFirst = 0;
@@ -65 +61 @@
     source->extpars->isophot = pmSourceIsophotalValuesAlloc ();
   }
-  
+
   // these are uncalibrated: instrumental mags and pixel units
   source->extpars->isophot->mag    = -2.5*log10(isophotalFlux);
@@ -74 +70 @@
 
   psTrace ("psphot", 5, "Isophot flux:%f +/- %f @ %f +/- %f for %f, %f\n",
-           source->extpars->isophot->mag, source->extpars->isophot->magErr,
-           source->extpars->isophot->rad, source->extpars->isophot->radErr,
-           source->peak->xf, source->peak->yf);
+           source->extpars->isophot->mag, source->extpars->isophot->magErr,
+           source->extpars->isophot->rad, source->extpars->isophot->radErr,
+           source->peak->xf, source->peak->yf);
 
   return true;