summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--ChangeLog7
-rw-r--r--complex.c19
-rw-r--r--rational.c3
3 files changed, 26 insertions, 3 deletions
diff --git a/ChangeLog b/ChangeLog
index 3af2e62..194dff7 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,3 +1,10 @@
+Sun Jun 28 22:25:07 2009 Tadayoshi Funaba <tadf@dotrb.org>
+
+ * complex.c (nucomp_expt): convert to a float when the given power
+ is a bignum.
+
+ * rational.c (nurat_expt): ditto.
+
Sun Jun 28 21:16:48 2009 Tadayoshi Funaba <tadf@dotrb.org>
* lib/cmath.rb (sqrt): fixed an issue [ruby-list:45852].
diff --git a/complex.c b/complex.c
index 2a7e602..3863e68 100644
--- a/complex.c
+++ b/complex.c
@@ -227,6 +227,18 @@ k_integer_p(VALUE x)
}
inline static VALUE
+k_fixnum_p(VALUE x)
+{
+ return f_kind_of_p(x, rb_cFixnum);
+}
+
+inline static VALUE
+k_bignum_p(VALUE x)
+{
+ return f_kind_of_p(x, rb_cBignum);
+}
+
+inline static VALUE
k_float_p(VALUE x)
{
return f_kind_of_p(x, rb_cFloat);
@@ -831,7 +843,7 @@ nucomp_expt(VALUE self, VALUE other)
f_mul(dat->imag, m_log_bang(r)));
return f_complex_polar(CLASS_OF(self), nr, ntheta);
}
- if (k_integer_p(other)) {
+ if (k_fixnum_p(other)) {
if (f_gt_p(other, ZERO)) {
VALUE x, z, n;
@@ -862,9 +874,12 @@ nucomp_expt(VALUE self, VALUE other)
if (k_numeric_p(other) && f_real_p(other)) {
VALUE r, theta;
+ if (k_bignum_p(other))
+ rb_warn("in a**b, b may be too big");
+
r = f_abs(self);
theta = f_arg(self);
- return f_complex_polar(CLASS_OF(self), f_expt(r, other),
+ return f_complex_polar(CLASS_OF(self), rb_fexpt(r, other),
f_mul(theta, other));
}
return rb_num_coerce_bin(self, other, id_expt);
diff --git a/rational.c b/rational.c
index 90025fa..8b28ac5 100644
--- a/rational.c
+++ b/rational.c
@@ -902,7 +902,6 @@ nurat_expt(VALUE self, VALUE other)
switch (TYPE(other)) {
case T_FIXNUM:
- case T_BIGNUM:
{
VALUE num, den;
@@ -924,6 +923,8 @@ nurat_expt(VALUE self, VALUE other)
}
return f_rational_new2(CLASS_OF(self), num, den);
}
+ case T_BIGNUM:
+ rb_warn("in a**b, b may be too big");
case T_FLOAT:
case T_RATIONAL:
return rb_fexpt(f_to_f(self), other);