[ruby-dev:11414] Re: [PATCH] an analogue of `long long'

まつもと ゆきひろです

In message "[ruby-dev:11410] Re: [PATCH] an analogue of `long long'"
    on 00/11/07, Masahiro Tanaka <[e-mail]> writes:

|> と思ったのですが、この部分、乗算・除算を行っていたのを左右ビッ
|> トシフトに置き換えたら、かなり速くなりました。やはり人の意見
|> には耳を傾けた方が良いようです。^^;;;
|
|あらら、そうですか。私がおかしいのかな。

どうなんでしょう。私のところでは1.6.1に比べて倍近く速いみた
いです。修正はこんな感じです。

--- bignum.c	2000/11/02 09:04:53	1.26
+++ bignum.c	2000/11/07 16:10:01
@@ -832,3 +832,10 @@
     while (!yds[ny-1]) ny--;
-    if ((dd = BIGRAD/(BDIGIT_DBL_SIGNED)(yds[ny-1]+1)) != 1) {
+
+    dd = 0;
+    q = yds[ny-1];
+    while ((q & (1<<(BITSPERDIG-1))) == 0) {
+	q <<= 1;
+	dd++;
+    }
+    if (dd) {
 	yy = rb_big_clone(y);
@@ -836,7 +843,7 @@
 	j = 0;
-	num = 0;
+	t2 = 0;
 	while (j<ny) {
-	    num += (BDIGIT_DBL)yds[j]*dd;
-	    tds[j++] = BIGLO(num);
-	    num = BIGDN(num);
+	    t2 += (BDIGIT_DBL)yds[j]<<dd;
+	    tds[j++] = BIGLO(t2);
+	    t2 = BIGDN(t2);
 	}
@@ -844,9 +851,9 @@
 	j = 0;
-	num = 0;
+	t2 = 0;
 	while (j<nx) {
-	    num += (BDIGIT_DBL)xds[j]*dd;
-	    zds[j++] = BIGLO(num);
-	    num = BIGDN(num);
+	    t2 += (BDIGIT_DBL)xds[j]<<dd;
+	    zds[j++] = BIGLO(t2);
+	    t2 = BIGDN(t2);
 	}
-	zds[j] = (BDIGIT)num;
+	zds[j] = (BDIGIT)t2;
     }
@@ -899,5 +907,6 @@
 	    while(i--) {
-		t2 = BIGUP(t2) + zds[i];
-		zds[i] = (BDIGIT)(t2 / dd);
-		t2 %= dd;
+		t2 = (t2 | zds[i]) >> dd;
+		q = zds[i];
+		zds[i] = BIGLO(t2);
+		t2 = BIGUP(q);
 	    }