Changeset 6751 for branches/rel10_ifa/psModules/src/objects/pmGrowthCurve.c

Timestamp:

Mar 31, 2006, 4:47:20 PM (20 years ago)

Author:

magnier

Message:

fixed up low-level issues with sky levels, sky bias

File:

• branches/rel10_ifa/psModules/src/objects/pmGrowthCurve.c (modified) (4 diffs)

branches/rel10_ifa/psModules/src/objects/pmGrowthCurve.c

@@ -5 +5 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.1.4.2 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-03-09 03:14:23 $
+ *  @version $Revision: 1.1.4.3 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-04-01 02:47:20 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -26 +26 @@
 }
 
-pmGrowthCurve *pmGrowthCurveAlloc (psF32 minRadius, psF32 maxRadius, psF32 dRadius)
+pmGrowthCurve *pmGrowthCurveAlloc (psF32 minRadius, psF32 maxRadius, psF32 refRadius)
 {
 
@@ -32 +32 @@
     psMemSetDeallocator(growth, (psFreeFunc) pmGrowthCurveFree);
 
-    growth->radius = psVectorCreate (NULL, minRadius, maxRadius, dRadius, PS_TYPE_F32);
+    // fractional pixel radii are not well defined; use integer steps
+    growth->radius = psVectorCreate (NULL, minRadius, maxRadius, 1.0, PS_TYPE_F32);
     growth->apMag  = psVectorAlloc (growth->radius->n, PS_TYPE_F32);
 
     // XXX may want to extend this to allow for a different refRadius;
-    growth->refRadius = maxRadius;
+    growth->refRadius = refRadius;
     growth->maxRadius = maxRadius;
     growth->apLoss = 0.0;
@@ -111 +112 @@
 
     if (n1 == 0) {
-        int n0 = psVectorBracket (index, key, FALSE);
+        n0 = psVectorBracket (index, key, FALSE);
         n1 = n0 + 1;
     }