Changeset 25718 for branches/pap/psModules/src/objects/models/pmModel_PS1_V1.c

Timestamp:

Oct 1, 2009, 11:14:09 AM (17 years ago)

Author:

Paul Price

Message:

Explicit handling of PAR_7 == 0, otherwise radius = Inf.

File:

• branches/pap/psModules/src/objects/models/pmModel_PS1_V1.c (modified) (1 diff)

branches/pap/psModules/src/objects/models/pmModel_PS1_V1.c

@@ -265 +265 @@
     psF32 *PAR = params->data.F32;
 
-    if (flux <= 0)
-        return (1.0);
-    if (PAR[PM_PAR_I0] <= 0)
-        return (1.0);
-    if (flux >= PAR[PM_PAR_I0])
-        return (1.0);
+    if (flux <= 0) {
+        return 1.0;
+    }
+    if (PAR[PM_PAR_I0] <= 0) {
+        return 1.0;
+    }
+    if (flux >= PAR[PM_PAR_I0]) {
+        return 1.0;
+    }
+    if (PAR[PM_PAR_7] == 0.0) {
+        return powf(PAR[PM_PAR_I0] / flux - 1.0, 1.0 / ALPHA);
+    }
 
     shape.sx  = PAR[PM_PAR_SXX] / M_SQRT2;