Changeset 33411 for trunk/psphot

Timestamp:

Mar 5, 2012, 4:56:39 PM (14 years ago)

Author:

bills

Message:

don't extrapolate limitRadius if the result goes beyond MAX_RADIUS

File:

• trunk/psphot/src/psphotRadialProfileWings.c (modified) (1 diff)

trunk/psphot/src/psphotRadialProfileWings.c

@@ -307 +307 @@
             // linearly interpolate to the radius at which we hit the sky, using the last flux and the limiting slope
             if (isfinite(lastFlux)) {
-                limitRadius = lastRadius + lastFlux / 3.0;
+                float interpolatedRadius = lastRadius + lastFlux / 3.0;
+                if (interpolatedRadius < MAX_RADIUS) {
+                    limitRadius = interpolatedRadius;
+                } else {
+                    // XXX should we keep going in this case?
+                    limitRadius = meanRadius;
+                }
             } else {
                 limitRadius = meanRadius;