Changeset 25678

Timestamp:

Sep 29, 2009, 9:53:32 PM (17 years ago)

Author:

eugene

Message:

do not fail is PAR_7 is 0.0 or less

File:

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

branches/eam_branches/20090715/psModules/src/objects/models/pmModel_PS1_V1.c

@@ -323 +323 @@
 
     // choose a z value guaranteed to be beyond our limit
-    float z0 = pow((1.0 / limit), (1.0 / ALPHA));
-    float z1 = (1.0 / limit) / PAR[PM_PAR_7];
-    z1 = PS_MAX (z0, z1);
-    z0 = 0.0;
+    float z0 = 0.0;
+    float z1 = pow((1.0 / limit), (1.0 / ALPHA));
+    if (PAR[PM_PAR_7] < 0.0) z1 *= 2.0;
 
     // perform a type of bisection to find the value