[ruby-list:19891] Re: build on alpha (1/2)

えぐち＠エスアンドイー です。

＃手違いで [ruby-list:19856] なくしてしまい
＃body は復元したものの、 MID がおかしいかも知れません
＃スレッドが正しくつながらなかったら御容赦ください _o_

>>> In message [ruby-list:19856] Re: build on alpha (1/2)
    On Thu, 23 Dec 1999 14:00:15 +0900, 中川 said:

中川> 中川です。
中川> 
中川> In message "[ruby-list:19825] Re: build on alpha (1/2)"
中川>     on 99/12/22, EGUCHI Osamu <[e-mail]> wrote:
中川> 
中川>   > ＃引用元のマルチパートのエンコードに失敗してしまったので、
中川>   > ＃御見苦しいかも知れませんが、御勘弁を。
中川> 
中川> ごめんなさい。encoding を x-gzip64 で split されてしまったので元に戻らな
中川> かったみたいです。
中川> 
中川>   >   if (n1**2-1) / (n1+1) != (n1-1)
中川> 
中川>   > の / が事故現場だと思います。
中川>   > とすると、 n1 + 1 が 0 になっているのかなぁ？
中川> 
中川> ログを見ていただければ分りますが、n1+1 が 0 になっているわけではなさそう
中川> です。それより << の処理に問題がありそうです。

32bit マシンで 64bitマシンの状況を再現させる方法を、週末考えてみました。

考え付いた、再現と対策の方法は

  USHORT を unsigned int (32bit) にし、桁溢れが発生させる

  コンパイラの警告や make test でエラーが出た部分を
  long long にかえたりキャストしたりする事で黙らせる

というものです。(我ながら力業)

この方法で、 make test の marshal 以外の項目をパスする事が出来ました。
＃面白い事に、加減算などは数％から十数％速度が向上しました
＃が乗除算は数十％遅くなったです ^^;;

このメールのパッチは、 cvs の HEAD へのパッチです。
＃このパッチでは、long long を元の long に直してあります。

	えぐち

Index: bignum.c
===================================================================
RCS file: /usr/home/eguchi/cvs.netlab.co.jp/ruby/bignum.c,v
retrieving revision 1.4
diff -d -u -r1.4 bignum.c
--- bignum.c	1999/12/06 09:03:49	1.4
+++ bignum.c	1999/12/26 13:53:40
@@ -13,15 +13,14 @@
 #include <ctype.h>
 
 VALUE rb_cBignum;
-typedef unsigned short USHORT;
 
 #define BDIGITS(x) RBIGNUM(x)->digits
-#define BITSPERDIG (sizeof(short)*CHAR_BIT)
-#define BIGRAD (1L << BITSPERDIG)
-#define DIGSPERINT ((unsigned int)(sizeof(long)/sizeof(short)))
+#define BITSPERDIG (sizeof(USHORT)*CHAR_BIT)
+#define BIGRAD (((unsigned long)1) << BITSPERDIG)
+#define DIGSPERINT ((unsigned int)(sizeof(long)/sizeof(USHORT)))
 #define BIGUP(x) ((unsigned long)(x) << BITSPERDIG)
 #define BIGDN(x) RSHIFT(x,BITSPERDIG)
-#define BIGLO(x) ((USHORT)((x) & (BIGRAD-1)))
+#define BIGLO(x) ((USHORT)((x)))
 
 static VALUE
 bignew_1(klass, len, sign)
@@ -263,7 +262,7 @@
 	num = c;
 	for (;;) {
 	    while (i<blen) {
-		num += zds[i]*base;
+		num += (unsigned long)zds[i] * base;
 		zds[i++] = BIGLO(num);
 		num = BIGDN(num);
 	    }
@@ -589,7 +588,7 @@
 
     len = RBIGNUM(x)->len;
     for (i = 0, num = 0; i < len; i++) {
-	num += BDIGITS(x)[i] + BDIGITS(y)[i];
+	num += (unsigned long)BDIGITS(x)[i] + BDIGITS(y)[i];
 	BDIGITS(z)[i] = BIGLO(num);
 	num = BIGDN(num);
     }
@@ -680,9 +679,11 @@
 	if (dd == 0) continue;
 	n = 0;
 	for (j = 0; j < RBIGNUM(y)->len; j++) {
-	    int ee = n + dd * BDIGITS(y)[j];
-	    n = zds[i + j] + ee;
-	    if (ee) zds[i + j] = BIGLO(n);
+	    long ee = n + dd * BDIGITS(y)[j];
+	    if (ee) {
+		n = zds[i + j] + ee;
+		zds[i + j] = BIGLO(n);
+	    }
 	    n = BIGDN(n);
 	}
 	if (n) {
@@ -737,13 +738,13 @@
     zds = BDIGITS(z);
     if (nx==ny) zds[nx+1] = 0;
     while (!yds[ny-1]) ny--;
-    if ((dd = BIGRAD/(int)(yds[ny-1]+1)) != 1) {
+    if ((dd = BIGRAD / (yds[ny - 1] + 1)) != 1) {
 	yy = rb_big_clone(y);
 	tds = BDIGITS(yy);
 	j = 0;
 	num = 0;
 	while (j<ny) {
-	    num += (long)yds[j]*dd;
+	    num += (long)yds[j] * dd;
 	    tds[j++] = BIGLO(num);
 	    num = BIGDN(num);
 	}
@@ -751,7 +752,7 @@
 	j = 0;
 	num = 0;
 	while (j<nx) {
-	    num += (long)xds[j]*dd;
+	    num += (long)xds[j] * dd;
 	    zds[j++] = BIGLO(num);
 	    num = BIGDN(num);
 	}
@@ -769,11 +770,13 @@
 	if (q) {
 	    i = 0; num = 0; t2 = 0;
 	    do {			/* multiply and subtract */
-		int ee;
+		long ee;
 		t2 += (long)yds[i] * q;
 		ee = num - BIGLO(t2);
-		num = zds[j - ny + i] + ee;
-		if (ee) zds[j - ny + i] = BIGLO(num);
+		if (ee) {
+		    num = zds[j - ny + i] + ee;
+		    zds[j - ny + i] = BIGLO(num);
+		}
 		num = BIGDN(num);
 		t2 = BIGDN(t2);
 	    } while (++i < ny);
@@ -781,9 +784,11 @@
 	    while (num) {		/* "add back" required */
 		i = 0; num = 0; q--;
 		do {
-		    int ee = num + yds[i];
-		    num = (long) zds[j - ny + i] + ee;
-		    if (ee) zds[j - ny + i] = BIGLO(num);
+		    long ee = num + yds[i];
+		    if (ee) {
+			num = (long)zds[j - ny + i] + ee;
+			zds[j - ny + i] = BIGLO(num);
+		    }
 		    num = BIGDN(num);
 		} while (++i < ny);
 		num--;
@@ -1143,7 +1148,7 @@
 	*zds++ = 0;
     }
     for (i=0; i<len; i++) {
-	num = num | *xds++<<s2;
+	num = num | (unsigned long)*xds++ << s2;
 	*zds++ = BIGLO(num);
 	num = BIGDN(num);
     }