github
/
ruby
зеркало из https://github.com/github/ruby.git
1
0
Форкнуть 0
Код Задачи Пакеты Проекты Релизы Вики Активность

matz 

git-svn-id: svn+ssh://ci.ruby-lang.org/ruby/trunk@1026 b2dd03c8-39d4-4d8f-98ff-823fe69b080e
This commit is contained in:
matz 2000-11-01 08:49:40 +00:00
Родитель 3685e2bf2c
Коммит d1bcd482a7
2 изменённых файлов: 37 добавлений и 34 удалений
Показать статистику Скачать Patch файл Скачать Diff файл Показать все файлы Свернуть все файлы

68
bignum.c
Просмотреть файл

@ -22,23 +22,23 @@ VALUE rb_cBignum;
#if SIZEOF_LONG*2 <= SIZEOF_LONG_LONG #if SIZEOF_LONG*2 <= SIZEOF_LONG_LONG
typedef unsigned long BDIGIT; typedef unsigned long BDIGIT;
typedef unsigned long long BDIGIT2; typedef unsigned long long BDIGIT_DBL;
typedef long long BDIGIT2_SIGNED; typedef long long BDIGIT_DBL_SIGNED;
#elif SIZEOF_INT*2 <= SIZEOF_LONG_LONG #elif SIZEOF_INT*2 <= SIZEOF_LONG_LONG
typedef unsigned int BDIGIT; typedef unsigned int BDIGIT;
typedef unsigned long long BDIGIT2; typedef unsigned long long BDIGIT_DBL;
typedef long long BDIGIT2_SIGNED; typedef long long BDIGIT_DBL_SIGNED;
#else #else
typedef unsigned short BDIGIT; typedef unsigned short BDIGIT;
typedef unsigned long BDIGIT2; typedef unsigned long BDIGIT_DBL;
typedef long BDIGIT2_SIGNED; typedef long BDIGIT_DBL_SIGNED;
#endif #endif
#define BDIGITS(x) ((BDIGIT*)RBIGNUM(x)->digits) #define BDIGITS(x) ((BDIGIT*)RBIGNUM(x)->digits)
#define BITSPERDIG (sizeof(BDIGIT)*CHAR_BIT) #define BITSPERDIG (sizeof(BDIGIT)*CHAR_BIT)
#define BIGRAD ((BDIGIT2)1 << BITSPERDIG) #define BIGRAD ((BDIGIT_DBL)1 << BITSPERDIG)
#define DIGSPERLONG ((unsigned int)(sizeof(long)/sizeof(BDIGIT))) #define DIGSPERLONG ((unsigned int)(sizeof(long)/sizeof(BDIGIT)))
#define BIGUP(x) ((BDIGIT2)(x) << BITSPERDIG) #define BIGUP(x) ((BDIGIT_DBL)(x) << BITSPERDIG)
#define BIGDN(x) RSHIFT(x,BITSPERDIG) #define BIGDN(x) RSHIFT(x,BITSPERDIG)
#define BIGLO(x) ((BDIGIT)((x) & (BIGRAD-1))) #define BIGLO(x) ((BDIGIT)((x) & (BIGRAD-1)))
@ -75,7 +75,7 @@ rb_big_2comp(x) /* get 2's complement */
{ {
long i = RBIGNUM(x)->len; long i = RBIGNUM(x)->len;
BDIGIT *ds = BDIGITS(x); BDIGIT *ds = BDIGITS(x);
BDIGIT2 num; BDIGIT_DBL num;
while (i--) ds[i] = ~ds[i]; while (i--) ds[i] = ~ds[i];
i = 0; num = 1; i = 0; num = 1;
@ -132,7 +132,7 @@ VALUE
rb_uint2big(n) rb_uint2big(n)
unsigned long n; unsigned long n;
{ {
BDIGIT2 num = n; BDIGIT_DBL num = n;
long i = 0; long i = 0;
BDIGIT *digits; BDIGIT *digits;
VALUE big; VALUE big;
@ -194,7 +194,7 @@ rb_cstr2inum(str, base)
char *end; char *end;
int badcheck = (base==0)?1:0; int badcheck = (base==0)?1:0;
char sign = 1, c; char sign = 1, c;
BDIGIT2 num; BDIGIT_DBL num;
long len, blen = 1; long len, blen = 1;
long i; long i;
VALUE z; VALUE z;
@ -306,7 +306,7 @@ rb_cstr2inum(str, base)
num = c; num = c;
for (;;) { for (;;) {
while (i<blen) { while (i<blen) {
num += (BDIGIT2)zds[i]*base; num += (BDIGIT_DBL)zds[i]*base;
zds[i++] = BIGLO(num); zds[i++] = BIGLO(num);
num = BIGDN(num); num = BIGDN(num);
} }
@ -390,7 +390,7 @@ rb_big2str(x, base)
s[0] = RBIGNUM(x)->sign ? '+' : '-'; s[0] = RBIGNUM(x)->sign ? '+' : '-';
while (i && j) { while (i && j) {
long k = i; long k = i;
BDIGIT2 num = 0; BDIGIT_DBL num = 0;
while (k--) { while (k--) {
num = BIGUP(num) + ds[k]; num = BIGUP(num) + ds[k];
@ -427,7 +427,7 @@ big2ulong(x, type)
char *type; char *type;
{ {
long len = RBIGNUM(x)->len; long len = RBIGNUM(x)->len;
BDIGIT2 num; BDIGIT_DBL num;
BDIGIT *ds; BDIGIT *ds;
if (len > sizeof(long)/sizeof(BDIGIT)) if (len > sizeof(long)/sizeof(BDIGIT))
@ -617,7 +617,7 @@ bigsub(x, y)
{ {
VALUE z = 0; VALUE z = 0;
BDIGIT *zds; BDIGIT *zds;
BDIGIT2_SIGNED num; BDIGIT_DBL_SIGNED num;
long i; long i;
i = RBIGNUM(x)->len; i = RBIGNUM(x)->len;
@ -642,7 +642,7 @@ bigsub(x, y)
zds = BDIGITS(z); zds = BDIGITS(z);
for (i = 0, num = 0; i < RBIGNUM(y)->len; i++) { for (i = 0, num = 0; i < RBIGNUM(y)->len; i++) {
num += (BDIGIT2_SIGNED)BDIGITS(x)[i] - BDIGITS(y)[i]; num += (BDIGIT_DBL_SIGNED)BDIGITS(x)[i] - BDIGITS(y)[i];
zds[i] = BIGLO(num); zds[i] = BIGLO(num);
num = BIGDN(num); num = BIGDN(num);
} }
@ -665,7 +665,7 @@ bigadd(x, y, sign)
char sign; char sign;
{ {
VALUE z; VALUE z;
BDIGIT2 num; BDIGIT_DBL num;
long i, len; long i, len;
sign = (sign == RBIGNUM(y)->sign); sign = (sign == RBIGNUM(y)->sign);
@ -685,7 +685,7 @@ bigadd(x, y, sign)
len = RBIGNUM(x)->len; len = RBIGNUM(x)->len;
for (i = 0, num = 0; i < len; i++) { for (i = 0, num = 0; i < len; i++) {
num += BDIGITS(x)[i] + BDIGITS(y)[i]; num += (BDIGIT_DBL)BDIGITS(x)[i] + BDIGITS(y)[i];
BDIGITS(z)[i] = BIGLO(num); BDIGITS(z)[i] = BIGLO(num);
num = BIGDN(num); num = BIGDN(num);
} }
@ -747,7 +747,7 @@ rb_big_mul(x, y)
VALUE x, y; VALUE x, y;
{ {
long i, j; long i, j;
BDIGIT2 n = 0; BDIGIT_DBL n = 0;
VALUE z; VALUE z;
BDIGIT *zds; BDIGIT *zds;
@ -772,11 +772,11 @@ rb_big_mul(x, y)
zds = BDIGITS(z); zds = BDIGITS(z);
while (j--) zds[j] = 0; while (j--) zds[j] = 0;
for (i = 0; i < RBIGNUM(x)->len; i++) { for (i = 0; i < RBIGNUM(x)->len; i++) {
BDIGIT2 dd = BDIGITS(x)[i]; BDIGIT_DBL dd = BDIGITS(x)[i];
if (dd == 0) continue; if (dd == 0) continue;
n = 0; n = 0;
for (j = 0; j < RBIGNUM(y)->len; j++) { for (j = 0; j < RBIGNUM(y)->len; j++) {
BDIGIT2 ee = n + (BDIGIT2)dd * BDIGITS(y)[j]; BDIGIT_DBL ee = n + (BDIGIT_DBL)dd * BDIGITS(y)[j];
n = zds[i + j] + ee; n = zds[i + j] + ee;
if (ee) zds[i + j] = BIGLO(n); if (ee) zds[i + j] = BIGLO(n);
n = BIGDN(n); n = BIGDN(n);
@ -798,8 +798,8 @@ bigdivrem(x, y, divp, modp)
long i, j; long i, j;
VALUE yy, z; VALUE yy, z;
BDIGIT *xds, *yds, *zds, *tds; BDIGIT *xds, *yds, *zds, *tds;
BDIGIT2 t2; BDIGIT_DBL t2;
BDIGIT2_SIGNED num; BDIGIT_DBL_SIGNED num;
BDIGIT dd, q; BDIGIT dd, q;
yds = BDIGITS(y); yds = BDIGITS(y);
@ -830,13 +830,13 @@ bigdivrem(x, y, divp, modp)
zds = BDIGITS(z); zds = BDIGITS(z);
if (nx==ny) zds[nx+1] = 0; if (nx==ny) zds[nx+1] = 0;
while (!yds[ny-1]) ny--; while (!yds[ny-1]) ny--;
if ((dd = BIGRAD/(BDIGIT2_SIGNED)(yds[ny-1]+1)) != 1) { if ((dd = BIGRAD/(BDIGIT_DBL_SIGNED)(yds[ny-1]+1)) != 1) {
yy = rb_big_clone(y); yy = rb_big_clone(y);
tds = BDIGITS(yy); tds = BDIGITS(yy);
j = 0; j = 0;
num = 0; num = 0;
while (j<ny) { while (j<ny) {
num += (BDIGIT2)yds[j]*dd; num += (BDIGIT_DBL)yds[j]*dd;
tds[j++] = BIGLO(num); tds[j++] = BIGLO(num);
num = BIGDN(num); num = BIGDN(num);
} }
@ -844,7 +844,7 @@ bigdivrem(x, y, divp, modp)
j = 0; j = 0;
num = 0; num = 0;
while (j<nx) { while (j<nx) {
num += (BDIGIT2)xds[j]*dd; num += (BDIGIT_DBL)xds[j]*dd;
zds[j++] = BIGLO(num); zds[j++] = BIGLO(num);
num = BIGDN(num); num = BIGDN(num);
} }
@ -862,20 +862,20 @@ bigdivrem(x, y, divp, modp)
if (q) { if (q) {
i = 0; num = 0; t2 = 0; i = 0; num = 0; t2 = 0;
do { /* multiply and subtract */ do { /* multiply and subtract */
BDIGIT2 ee; BDIGIT_DBL ee;
t2 += (BDIGIT2)yds[i] * q; t2 += (BDIGIT_DBL)yds[i] * q;
ee = num - BIGLO(t2); ee = num - BIGLO(t2);
num = zds[j - ny + i] + ee; num = (BDIGIT_DBL)zds[j - ny + i] + ee;
if (ee) zds[j - ny + i] = BIGLO(num); if (ee) zds[j - ny + i] = BIGLO(num);
num = BIGDN(num); num = BIGDN(num);
t2 = BIGDN(t2); t2 = BIGDN(t2);
} while (++i < ny); } while (++i < ny);
num += zds[j - ny + i] - t2; /* borrow from high digit; don't update */ num += zds[j - ny + i] - t2;/* borrow from high digit; don't update */
while (num) { /* "add back" required */ while (num) { /* "add back" required */
i = 0; num = 0; q--; i = 0; num = 0; q--;
do { do {
BDIGIT2 ee = num + yds[i]; BDIGIT_DBL ee = num + yds[i];
num = (BDIGIT2)zds[j - ny + i] + ee; num = (BDIGIT_DBL)zds[j - ny + i] + ee;
if (ee) zds[j - ny + i] = BIGLO(num); if (ee) zds[j - ny + i] = BIGLO(num);
num = BIGDN(num); num = BIGDN(num);
} while (++i < ny); } while (++i < ny);
@ -1230,7 +1230,7 @@ rb_big_lshift(x, y)
int s1 = shift/BITSPERDIG; int s1 = shift/BITSPERDIG;
int s2 = shift%BITSPERDIG; int s2 = shift%BITSPERDIG;
VALUE z; VALUE z;
BDIGIT2 num = 0; BDIGIT_DBL num = 0;
long len, i; long len, i;
if (shift < 0) return rb_big_rshift(x, INT2FIX(-shift)); if (shift < 0) return rb_big_rshift(x, INT2FIX(-shift));
@ -1259,7 +1259,7 @@ rb_big_rshift(x, y)
int s1 = shift/BITSPERDIG; int s1 = shift/BITSPERDIG;
int s2 = shift%BITSPERDIG; int s2 = shift%BITSPERDIG;
VALUE z; VALUE z;
BDIGIT2 num = 0; BDIGIT_DBL num = 0;
long i = RBIGNUM(x)->len; long i = RBIGNUM(x)->len;
long j; long j;

3
sample/test.rb
Просмотреть файл

@ -641,6 +641,9 @@ test_ok(7 == a.remainder(-b))
test_ok(-7 == (-a).remainder(b)) test_ok(-7 == (-a).remainder(b))
test_ok(-7 == (-a).remainder(-b)) test_ok(-7 == (-a).remainder(-b))
test_ok(10**40+10**20 == 10000000000000000000100000000000000000000)
test_ok(10**40/10**20 == 100000000000000000000)
test_check "string & char" test_check "string & char"
test_ok("abcd" == "abcd") test_ok("abcd" == "abcd")