Bignum support, here we go!

Bignums, based on Michael Fromberger's MPI library, are integrated
into the input syntax, stream output, equality testing, the garbage
collector, and hashing.

The plus operation handles transitions between fixnums and bignums.
Other operations are still fixnum only.

* Makefile (CFLAGS): Add mpi directory to include file search.
(OBJS): Include new arith.o module and all of MPI_OBJS.
(MPI_OBJS, MPI_OBJS_BASE): New variables.

* configure (mpi_version, have_quilt, have_patch): New variables.
Script detects whether patch and quilt are available. Unpacks
mpi library, applies patches. Detects 128 bit integer type.
Records more information in config.h about the sizes of types.

* dep.mk: Updated.

* depend.txr: Make work with paths that have directory components.

* eval.c (eval_init): Rename of nump to fixnump.

* gc.c (finalize, mark_obj): Handle BGNUM case.

* hash.c: (hash_c_str): Changed to return unsigned long
instead of long.
(equal_hash): Handle BGNUM case.
(eql_hash): Handle bignums with equal-hash, but other
objects as eq.

* lib.c (num_s): Variable renamed to fixnum_s.
(bignum_s): New symbol variable.
(code2type): Follow rename of num_s. Handle BGNUM case.
(typeof): Follow rename of num_s.
(eql): Handle bignums using equal, and other types using eq.
(equal): Handle BGNUM case.
(chk_calloc): New function.
(c_num): Wording change in error message: is not a fixnum.
(nump): Renamed to fixnump.
(bignump): New function.
(plus): Function removed, reimplemented in arith.c.
(int_str): Handle integers which are too large for wcstol
using bignum conversion. Base 0 is no longer passed to
wcstol but converted to 10 because the special semantics
for 0 would be inconsistent for bignums.
(obj_init): Follow rename of num_s. Initialize bignum_s.

1 files changed, 48 insertions, 0 deletions

diff --git a/mpi-patches/add-mp-hash b/mpi-patches/add-mp-hash
new file mode 100644
index 00000000..8a2cd585
--- /dev/null
+++ b/mpi-patches/add-mp-hash
@@ -0,0 +1,48 @@
+Index: mpi-1.8.6/mpi.c
+===================================================================
+--- mpi-1.8.6.orig/mpi.c	2011-12-09 14:10:41.000000000 -0800
++++ mpi-1.8.6/mpi.c	2011-12-09 14:26:02.000000000 -0800
+@@ -1960,6 +1960,30 @@
+ 
+ /* }}} */
+ 
++unsigned long mp_hash(mp_int *a)
++{
++  unsigned long hash = 0;
++  int ix;
++  for (ix = 0; ix < USED(a); ix++) {
++    mp_digit d = DIGIT(a, ix);
++#if SIZEOF_LONG < MP_DIGIT_SIZE
++    int j;
++    for (j = 0; j < MP_DIGIT_SIZE / SIZEOF_LONG; j++) {
++	hash ^= d;
++	d >> (SIZEOF_LONG * CHAR_BIT);
++    }
++#elif SIZEOF_LONG == MP_DIGIT_SIZE
++    hash ^= d;
++#else
++    hash <<= MP_DIGIT_BITS;
++    hash ^= d;
++#endif
++  }
++  if (SIGN(a) == MP_NEG)
++    hash = (hash << 16 | hash >> (SIZEOF_LONG * CHAR_BIT - 16));
++  return hash;
++}
++
+ /*------------------------------------------------------------------------*/
+ /* {{{ Number theoretic functions */
+ 
+Index: mpi-1.8.6/mpi.h
+===================================================================
+--- mpi-1.8.6.orig/mpi.h	2011-12-09 14:10:41.000000000 -0800
++++ mpi-1.8.6/mpi.h	2011-12-09 14:10:41.000000000 -0800
+@@ -165,6 +165,8 @@
+ int    mp_isodd(mp_int *a);
+ int    mp_iseven(mp_int *a);
+ 
++unsigned long mp_hash(mp_int *a);
++
+ /*------------------------------------------------------------------------*/
+ /* Number theoretic                                                       */
+