summaryrefslogtreecommitdiff
path: root/bignum.c
diff options
context:
space:
mode:
authorYukihiro Matsumoto <matz@ruby-lang.org>1994-08-10 15:54:46 +0900
committerTakashi Kokubun <takashikkbn@gmail.com>2019-08-17 22:09:30 +0900
commit6e3090413652b6592346556149fed1e9aec5495d (patch)
treebac97139bbeedc8cb67cb2e451a22ed4ddb2b2d4 /bignum.c
parent200e0ee2fd3c1c006c528874a88f684447215524 (diff)
version 0.50v0_50
http://cache.ruby-lang.org/pub/ruby/1.0/ruby-0.50.tar.gz Wed Aug 10 15:54:46 1994 Yukihiro Matsumoto (matz@ix-02) * variable.c: -vオプションが指定されている時は初期化されていない, 大域変数, インスタンス変数, ローカル変数を参照した時点でwarning を出すようにした. Tue Aug 9 11:50:48 1994 Yukihiro Matsumoto (matz@ix-02) * bignum.c: 冪乗に関しても多倍長演算を行なうように. 特に浮動小数点 数の範囲を越えた時の処理を的確に行なうように. * eval.c: メソッド定義後は構文木から, メソッド定義部分を外す. 無駄 な再定義が起こらないようにするためと2重にfree()されないため. * array.c(Fary_aref): 引数が1つでFixnumの時, Range checkを行なわな いように修正. * eval.c: 引数の数をコンパイル時に計算して若干の高速化. Mon Aug 8 13:06:24 1994 Yukihiro Matsumoto (matz@ix-02) * object.c: nilによる比較連鎖をなくした. * parse.y: bit演算子の優先順位を比較演算子よりも強くした. Cとは異 なることになるが, 直観には合致する. * gc.c: クラスを解放する時, 個々のメソッド毎にキャッシュをクリアす るのではなく, クラス単位でクリアするように. Thu Aug 4 18:45:09 1994 Yukihiro Matsumoto (matz@ix-02) * methods.c(method_free): 解放されたメソッドに関してキャッシュをク リアしておく必要があった. * gc.c: Dataクラスのデータ部分をfree()し忘れていた. Wed Aug 3 09:58:14 1994 Yukihiro Matsumoto (matz@ix-02) * parse.y: def func .. end形式による関数メソッドの定義はなくなった. * methods.c: func形式のメソッドをなくした. あっても, あまり意味が ないので. * eval.c: $0への代入でps(1)の出力が変化するように. * io.c(Fsyscall): syscall()を実現. Mon Aug 1 13:41:11 1994 Yukihiro Matsumoto (matz@ix-02) * parse.y: ダブルクォートで囲まれた文字列や正規表現中で"#{変数名}" または"#変数名"という形式で変数の内容を埋め込むことができるよう になった. * io.c: 関数メソッドsystem2()はなくなった. 今はバッククォートがあ るからね. * parse.y: `cmd`によってコマンドを文字列に展開することができるよう になった. * parse.y: __FILE__, __LINE__を追加. それぞれファイル名(文字列), 行番号(整数)を値とする疑似変数. Fri Jul 29 13:16:07 1994 Yukihiro Matsumoto (matz@ix-02) * methods.h: メソッドをオブジェクトとして扱うのをやめる. メソッド のメモリ管理にはリファレンスカウントを使うことにした. これでオブ ジェクトの数が減ってほんの少しだけGCが速くなる(かな). * purifyによってメモリ関係のバグを検査した(見つかる,見つかる…). * gc.c: GCをプログラマが変数をマークする形式から, スタックとレジス タからマークする方法に変更. 移植性が下がるような気もするが, siod やscmでも採用されているから多分大丈夫だろう. Linux on i486でも動 作を確認した. Wed Jul 27 16:13:13 1994 Yukihiro Matsumoto (matz@ix-02) * eval.c(Eval): トップレベルでは構造木をfreeしないように. どうせ解 放されるから時間の無駄である. * array.c, dict.c: "=="を構造一致に変更. Fri Jul 22 10:14:09 1994 Yukihiro Matsumoto (matz@ix-02) * error.c: 組み込みタイプの名前を登録し忘れていた. Thu Jul 21 14:06:48 1994 Yukihiro Matsumoto (matz@ix-02) * parse.y(freenode),eval.c(Eval): 解析木を解放し忘れていた. Mon Jul 18 10:19:15 1994 Yukihiro Matsumoto (matz@ix-02) * parse.y: 多重代入を処理するルールにバグがあって, 3要素以上の多重 代入に失敗していた. * eval.c(rb_eval): 多重代入で, 右辺が配列でない時には`to_a'メソッ ドで配列に変換して代入するようにした. 今までの仕様だと右辺値が第 1要素にそのまま代入されていたが, structなど配列に変換できるもの は変換した方が嬉しい気がする. * dbm.c,dict.c(delete_if): メソッド追加. * process.c(wait,waitpid): システムコールwaitpidまたはwait4がある 時はそちらを使うように. configureもそれらをチェックするように変更. * dbm.c, dict.c(clear): メソッド追加.
Diffstat (limited to 'bignum.c')
-rw-r--r--bignum.c107
1 files changed, 44 insertions, 63 deletions
diff --git a/bignum.c b/bignum.c
index ce124a6..c149242 100644
--- a/bignum.c
+++ b/bignum.c
@@ -10,6 +10,7 @@
#include "ruby.h"
#include <ctype.h>
+#include <math.h>
extern VALUE C_Integer;
VALUE C_Bignum;
@@ -87,7 +88,10 @@ bignorm(x)
while (len-- && !ds[len]) ;
x->len = ++len;
- if (len*sizeof(USHORT) <= sizeof(VALUE)) {
+
+ if (len*sizeof(USHORT) < sizeof(VALUE) ||
+ (len*sizeof(USHORT) == sizeof(VALUE) &&
+ ds[sizeof(VALUE)/sizeof(USHORT)-1] <= 0x3fff)) {
long num = 0;
while (len--) {
num = BIGUP(num) + ds[len];
@@ -282,10 +286,9 @@ big2str(x, base)
Fail("bignum cannot treat base %d", base);
}
- GC_LINK;
- GC_PRO3(t, Fbig_clone(x));
+ t = Fbig_clone(x);
ds = BDIGITS(t);
- GC_PRO3(ss, str_new(0, j));
+ ss = str_new(0, j);
s = RSTRING(ss)->ptr;
s[0] = x->sign ? '+' : '-';
@@ -310,7 +313,7 @@ big2str(x, base)
RSTRING(ss)->len -= x->sign?j:j-1;
memmove(x->sign?s:s+1, s+j, RSTRING(ss)->len);
s[RSTRING(ss)->len] = '\0';
- GC_UNLINK;
+
return ss;
}
@@ -363,20 +366,19 @@ dbl2big(d)
VALUE z;
double u = (d < 0)?-d:d;
- while (0 != floor(u)) {
+ while (0 != (long)u) {
u /= BIGRAD;
i++;
}
- GC_LINK;
- GC_PRO3(z, bignew(i, d>=0));
+ z = bignew(i, d>=0);
digits = BDIGITS(z);
while (i--) {
u *= BIGRAD;
- c = floor(u);
+ c = (long)u;
u -= c;
digits[i] = c;
}
- GC_UNLINK;
+
return bignorm(z);
}
@@ -430,8 +432,6 @@ bigadd(x, y, sign)
i = y->len;
while (i--) zds[i] = BDIGITS(y)[i];
- GC_LINK;
- GC_PRO(z);
i = 0; num = 0;
if (x->sign == z->sign) {
do {
@@ -481,7 +481,7 @@ bigadd(x, y, sign)
}
}
}
- GC_UNLINK;
+
return bignorm(z);
}
@@ -491,14 +491,12 @@ Fbig_plus(x, y)
{
VALUE z;
- GC_LINK;
- GC_PRO(y);
if (FIXNUM_P(y)) y = int2big(FIX2INT(y));
else {
Check_Type(x, T_BIGNUM);
}
z = bigadd(x, y, 1);
- GC_UNLINK;
+
return z;
}
@@ -506,14 +504,11 @@ VALUE
Fbig_minus(x, y)
VALUE x, y;
{
- GC_LINK;
- GC_PRO(y);
if (FIXNUM_P(y)) y = int2big(FIX2INT(y));
else {
Check_Type(y, T_BIGNUM);
}
x = bigadd(x, y, 0);
- GC_UNLINK;
return x;
}
@@ -527,8 +522,7 @@ Fbig_mul(x, y)
VALUE z;
USHORT *zds;
- GC_LINK;
- GC_PRO(y);
+ if (FIXNUM_P(x)) x = (struct RBignum*)int2big(FIX2INT(x));
if (FIXNUM_P(y)) y = (struct RBignum*)int2big(FIX2INT(y));
else {
Check_Type(y, T_BIGNUM);
@@ -552,7 +546,6 @@ Fbig_mul(x, y)
}
}
} while (++i < x->len);
- GC_UNLINK;
return bignorm(z);
}
@@ -579,8 +572,7 @@ bigdivmod(x, y, div, mod)
xds = BDIGITS(x);
if (ny == 1) {
dd = yds[0];
- GC_LINK;
- GC_PRO3(z, Fbig_clone(x));
+ z = Fbig_clone(x);
zds = BDIGITS(z);
t2 = 0; i = nx;
while(i--) {
@@ -593,16 +585,14 @@ bigdivmod(x, y, div, mod)
if (!y->sign) t2 = -t2;
*mod = FIX2INT(t2);
}
- GC_UNLINK;
return;
}
- GC_LINK;
- GC_PRO3(z, bignew(nx==ny?nx+2:nx+1, x->sign==y->sign));
+ z = bignew(nx==ny?nx+2:nx+1, x->sign==y->sign);
zds = BDIGITS(z);
if (nx==ny) zds[nx+1] = 0;
while (!yds[ny-1]) ny--;
if ((dd = BIGRAD/(yds[ny-1]+1)) != 1) {
- GC_PRO3(y, (struct RBignum*)Fbig_clone(y));
+ y = (struct RBignum*)Fbig_clone(y);
tds = BDIGITS(y);
j = 0;
num = 0;
@@ -681,7 +671,6 @@ bigdivmod(x, y, div, mod)
RBIGNUM(*mod)->sign = y->sign;
*mod = bignorm(*mod);
}
- GC_UNLINK;
}
static VALUE
@@ -690,14 +679,12 @@ Fbig_div(x, y)
{
VALUE z;
- GC_LINK;
- GC_PRO(y);
if (FIXNUM_P(y)) y = int2big(FIX2INT(y));
else {
Check_Type(y, T_BIGNUM);
}
bigdivmod(x, y, &z, Qnil);
- GC_UNLINK;
+
return z;
}
@@ -707,14 +694,12 @@ Fbig_mod(x, y)
{
VALUE z;
- GC_LINK;
- GC_PRO(y);
if (FIXNUM_P(y)) y = int2big(FIX2INT(y));
else {
Check_Type(y, T_BIGNUM);
}
bigdivmod(x, y, Qnil, &z);
- GC_UNLINK;
+
return z;
}
@@ -724,38 +709,41 @@ Fbig_divmod(x, y)
{
VALUE div, mod;
- GC_LINK;
- GC_PRO(y);
if (FIXNUM_P(y)) y = int2big(FIX2INT(y));
else {
Check_Type(y, T_BIGNUM);
}
bigdivmod(x, y, &div, &mod);
- GC_UNLINK;
return assoc_new(div, mod);;
}
-static VALUE
+VALUE
Fbig_pow(x, y)
VALUE x, y;
{
- extern double pow();
double d1, d2;
+ VALUE z;
+ int n;
- GC_LINK;
- GC_PRO(y);
- if (FIXNUM_P(y)) y = int2big(FIX2INT(y));
- else {
- Check_Type(y, T_BIGNUM);
+ if (TYPE(y) == T_FLOAT) {
+ return float_new(pow(big2dbl(x), RFLOAT(y)->value));
+ }
+ n = NUM2INT(y);
+ if (n == 0) return INT2FIX(1);
+ if (n < 0) {
+ return float_new(pow(big2dbl(x), (double)n));
}
- d1 = big2dbl(x);
- d2 = big2dbl(y);
- d1 = pow(d1, d2);
- GC_UNLINK;
-
- return dbl2big(d1);
+ z = x;
+ while (--n) {
+ while (!(n % 2)) {
+ n = n /2;
+ x = Fbig_mul(x, x);
+ }
+ z = Fbig_mul(z, x);
+ }
+ return z;
}
VALUE
@@ -774,14 +762,12 @@ Fbig_and(x, y)
Check_Type(y, T_BIGNUM);
}
- GC_LINK;
- GC_PRO(y);
if (!y->sign) {
y = (struct RBignum*)Fbig_clone(y);
big_2comp(y);
}
if (!x->sign) {
- GC_PRO3(x, (struct RBignum*)Fbig_clone(x));
+ x = (struct RBignum*)Fbig_clone(x);
big_2comp(x);
}
if (x->len > y->len) {
@@ -808,7 +794,6 @@ Fbig_and(x, y)
zds[i] = sign?0:ds2[i];
}
if (!RBIGNUM(z)->sign) big_2comp(z);
- GC_UNLINK;
return bignorm(z);
}
@@ -828,14 +813,12 @@ Fbig_or(x, y)
Check_Type(y, T_BIGNUM);
}
- GC_LINK;
- GC_PRO(y);
if (!y->sign) {
y = (struct RBignum*)Fbig_clone(y);
big_2comp(y);
}
if (!x->sign) {
- GC_PRO3(x, (struct RBignum*)Fbig_clone(x));
+ x = (struct RBignum*)Fbig_clone(x);
big_2comp(x);
}
if (x->len > y->len) {
@@ -862,7 +845,7 @@ Fbig_or(x, y)
zds[i] = sign?ds2[i]:(BIGRAD-1);
}
if (!RBIGNUM(z)->sign) big_2comp(z);
- GC_UNLINK;
+
return bignorm(z);
}
@@ -882,14 +865,12 @@ Fbig_xor(x, y)
Check_Type(y, T_BIGNUM);
}
- GC_LINK;
- GC_PRO(y);
if (!y->sign) {
y = (struct RBignum*)Fbig_clone(y);
big_2comp(y);
}
if (!x->sign) {
- GC_PRO3(x, (struct RBignum*)Fbig_clone(x));
+ x = (struct RBignum*)Fbig_clone(x);
big_2comp(x);
}
if (x->len > y->len) {
@@ -918,7 +899,7 @@ Fbig_xor(x, y)
zds[i] = sign?ds2[i]:~ds2[i];
}
if (!RBIGNUM(z)->sign) big_2comp(z);
- GC_UNLINK;
+
return bignorm(z);
}