# HG changeset patch
# User Micha¸«³l Cadilhac <michael.cadilhac@lrde.org>
# Date 1191230234 0
# Node ID 259a9fdd14501545d9d129002c4c3ba13302fe1a
# Parent  0157774b2cf92941c5fc0f67e2cddea821e8cd0a
(Fexpt): Manually check for overflows, so that a power
of a non-zero value can't yield zero.

diff -r 0157774b2cf9 -r 259a9fdd1450 src/floatfns.c
--- a/src/floatfns.c	Mon Oct 01 06:36:58 2007 +0000
+++ b/src/floatfns.c	Mon Oct 01 09:17:14 2007 +0000
@@ -454,7 +454,7 @@
      (arg1, arg2)
      register Lisp_Object arg1, arg2;
 {
-  double f1, f2;
+  double f1, f2, f3;
 
   CHECK_NUMBER_OR_FLOAT (arg1);
   CHECK_NUMBER_OR_FLOAT (arg2);
@@ -500,8 +500,11 @@
   else if ((f1 == 0.0 && f2 < 0.0) || (f1 < 0 && f2 != floor(f2)))
     domain_error2 ("expt", arg1, arg2);
 #endif
-  IN_FLOAT2 (f1 = pow (f1, f2), "expt", arg1, arg2);
-  return make_float (f1);
+  IN_FLOAT2 (f3 = pow (f1, f2), "expt", arg1, arg2);
+  /* Check for overflow in the result.  */
+  if (f1 != 0.0 && f3 == 0.0)
+    range_error ("expt", arg1);
+  return make_float (f3);
 }
 
 DEFUN ("log", Flog, Slog, 1, 2, 0,