From e9268149834f0d371d44bc8a612cd6ba51c00e3d Mon Sep 17 00:00:00 2001
From: tadf <tadf@b2dd03c8-39d4-4d8f-98ff-823fe69b080e>
Date: Sun, 16 Aug 2009 23:03:45 +0000
Subject: 	* numeric.c (flo_pow,fix_pow): may return complex number.

	* bignum.c (rb_big_pow): ditto.



git-svn-id: http://svn.ruby-lang.org/repos/ruby/trunk@24563 b2dd03c8-39d4-4d8f-98ff-823fe69b080e
---
 numeric.c | 17 +++++++++++++----
 1 file changed, 13 insertions(+), 4 deletions(-)

(limited to 'numeric.c')

diff --git a/numeric.c b/numeric.c
index 4ba69fd5d..dedd8e2ee 100644
--- a/numeric.c
+++ b/numeric.c
@@ -814,7 +814,6 @@ flo_divmod(VALUE x, VALUE y)
  * Raises <code>float</code> the <code>other</code> power.
  *
  *    2.0**3      #=> 8.0
- *    (-8.0)**0.5 #=> NaN  # try (-8.0)**Complex(0.5)
  */
 
 static VALUE
@@ -826,7 +825,13 @@ flo_pow(VALUE x, VALUE y)
       case T_BIGNUM:
 	return DBL2NUM(pow(RFLOAT_VALUE(x), rb_big2dbl(y)));
       case T_FLOAT:
-	return DBL2NUM(pow(RFLOAT_VALUE(x), RFLOAT_VALUE(y)));
+	{
+	    double dx = RFLOAT_VALUE(x);
+	    double dy = RFLOAT_VALUE(y);
+	    if (dx < 0 && dy != round(dy))
+		return rb_funcall(rb_complex_raw1(x), rb_intern("**"), 1, y);
+	    return DBL2NUM(pow(dx, dy));
+	}
       default:
 	return rb_num_coerce_bin(x, y, rb_intern("**"));
     }
@@ -2468,7 +2473,6 @@ int_pow(long x, unsigned long y)
  *    2 ** 3      #=> 8
  *    2 ** -1     #=> 0.5
  *    2 ** 0.5    #=> 1.4142135623731
- *    (-8)**0.5   #=> NaN  # try (-8)**Complex(0.5)
  */
 
 static VALUE
@@ -2517,7 +2521,12 @@ fix_pow(VALUE x, VALUE y)
 	    return DBL2NUM(RFLOAT_VALUE(y) < 0 ? infinite_value() : 0.0);
 	}
 	if (a == 1) return DBL2NUM(1.0);
-	return DBL2NUM(pow((double)a, RFLOAT_VALUE(y)));
+	{
+	    double dy = RFLOAT_VALUE(y);
+	    if (a < 0 && dy != round(dy))
+		return rb_funcall(rb_complex_raw1(x), rb_intern("**"), 1, y);
+	    return DBL2NUM(pow((double)a, dy));
+	}
       default:
 	return rb_num_coerce_bin(x, y, rb_intern("**"));
     }
-- 
cgit