Changeset 27565

Timestamp:

Apr 1, 2010, 3:19:02 PM (16 years ago)

Author:

eugene

Message:

do not allow z to become negative in radius search

Location:

trunk/psModules/src/objects/models

Files:

• pmModel_PGAUSS.c (modified) (1 diff)
• pmModel_PS1_V1.c (modified) (1 diff)
• pmModel_QGAUSS.c (modified) (1 diff)
• pmModel_RGAUSS.c (modified) (1 diff)

trunk/psModules/src/objects/models/pmModel_PGAUSS.c

@@ -287 +287 @@
         // fprintf (stderr, "%f %f %f : %f %f\n", f, z, dz, dfdz, q);
         z -= dz;
+        z = PS_MAX(z, 0.0);
     }
 

trunk/psModules/src/objects/models/pmModel_PS1_V1.c

@@ -313 +313 @@
         // fprintf (stderr, "%f %f %f : %f %f\n", f, z, dz, dfdz, q);
         z -= dz;
+        z = PS_MAX(z, 0.0);
     }
     psF64 radius = sigma * sqrt (2.0 * z);

trunk/psModules/src/objects/models/pmModel_QGAUSS.c

@@ -314 +314 @@
         // fprintf (stderr, "%f %f %f : %f %f\n", f, z, dz, dfdz, q);
         z -= dz;
+        z = PS_MAX(z, 0.0);
     }
     psF64 radius = sigma * sqrt (2.0 * z);

trunk/psModules/src/objects/models/pmModel_RGAUSS.c

@@ -307 +307 @@
         // fprintf (stderr, "%f %f %f : %f %f\n", f, z, dz, dfdz, q);
         z -= dz;
+        z = PS_MAX(z, 0.0);
     }