Changeset 2154

Timestamp:

07/04/08 02:38:03 (5 months ago)

Author:

gary

Message:

bugfix, force multiplier - use explicit test for div by zero

Files:

• trunk/source/radius_increment.cpp (modified) (1 diff)

trunk/source/radius_increment.cpp

@@ -503 +503 @@
         wind.AccelMax = (realnum)MAX2(wind.AccelMax,wind.AccelTot);
 
-        /* force multiplier; RELEC*T can be zero for very low densities */
-        if( relec*t > SMALLFLOAT )
-        {
-                wind.fmul = (realnum)((wind.AccelLine + wind.AccelCont)/(relec*t));
-        }
+        /* force multiplier; relec can be zero for very low densities - so use double
+         * form of safe_div - fmul is of order unity - wind.AccelLine and wind.AccelCont 
+         * are both floats to will underflow long before relec will - fmul is only used
+         * in output, not in any physics */
+        relec *= (EN1RYD/SPEEDLIGHT/dense.xMassDensity);
+        if( relec > SMALLFLOAT )
+                wind.fmul = (realnum)( (wind.AccelLine + wind.AccelCont) / relec);
         else
-        {
                 wind.fmul = 0.;
-        }
 
         /* keep track of average acceleration */