summaryrefslogtreecommitdiffstats
path: root/numeric.c
diff options
context:
space:
mode:
Diffstat (limited to 'numeric.c')
-rw-r--r--numeric.c34
1 files changed, 33 insertions, 1 deletions
diff --git a/numeric.c b/numeric.c
index 756f510ca..86728f206 100644
--- a/numeric.c
+++ b/numeric.c
@@ -2259,6 +2259,38 @@ fix_divmod(VALUE x, VALUE y)
}
}
+static VALUE
+int_pow(long x, unsigned long y)
+{
+ int neg = x < 0;
+ long z = 1;
+
+ if (neg) x = -x;
+ if (y & 1) z = x;
+ y &= ~1;
+ do {
+ while (y % 2 == 0) {
+ long x2 = x * x;
+ if (x2 < x || !POSFIXABLE(x2)) {
+ bignum:
+ return rb_big_mul(rb_big_pow(rb_int2big(x), LONG2NUM(y)),
+ rb_int2big(neg ? -z : z));
+ }
+ x = x2;
+ y >>= 1;
+ }
+ {
+ long xz = x * z;
+ if (xz < z || xz < x || !POSFIXABLE(xz)) {
+ goto bignum;
+ }
+ z = xz;
+ }
+ } while (--y);
+ if (neg) z = -z;
+ return LONG2NUM(z);
+}
+
/*
* call-seq:
* fix ** other => Numeric
@@ -2282,7 +2314,7 @@ fix_pow(VALUE x, VALUE y)
if (b == 1) return x;
a = FIX2LONG(x);
if (b > 0) {
- return rb_big_pow(rb_int2big(a), y);
+ return int_pow(a, b);
}
return rb_float_new(pow((double)a, (double)b));
}