2 files changed, 30 insertions, 0 deletions
diff --git a/mpi/mpi.c b/mpi/mpi.c
index 26334854..48c7203f 100644
--- a/mpi/mpi.c
+++ b/mpi/mpi.c
@@ -1958,6 +1958,34 @@ int    mp_iseven(mp_int *a)
 
 /* }}} */
 
+unsigned long mp_hash(mp_int *a)
+{
+#if SIZEOF_LONG > MP_DIGIT_SIZE
+  unsigned long hash;
+  int ix;
+
+  if (USED(a) >= 2 * SIZEOF_LONG / MP_DIGIT_SIZE) {
+    unsigned long omega = 0;
+    unsigned long alpha = 0;
+    for (ix = 0; ix < SIZEOF_LONG / MP_DIGIT_SIZE; ix++)
+	omega = (omega << MP_DIGIT_BIT) | DIGIT(a, ix);
+    for (ix = USED(a) - SIZEOF_LONG / MP_DIGIT_SIZE; ix < USED(a); ix++)
+	alpha = (alpha << MP_DIGIT_BIT) | DIGIT(a, ix);
+    hash = alpha + omega;
+  } else {
+    hash = 0;
+
+    for (ix = 0; ix < USED(a); ix++)
+	hash = (hash << MP_DIGIT_BIT) | DIGIT(a, ix);
+  }
+#else
+  mp_digit omega = DIGIT(a, 0);
+  mp_digit alpha = DIGIT(a, USED(a) - 1);
+  unsigned long hash = alpha + omega;
+#endif
+  return SIGN(a) == MP_NEG ? ~hash : hash;
+}
+
 /*------------------------------------------------------------------------*/
 /* {{{ Number theoretic functions */
 
diff --git a/mpi/mpi.h b/mpi/mpi.h
index c32b0d3a..51ffc00d 100644
--- a/mpi/mpi.h
+++ b/mpi/mpi.h
@@ -165,6 +165,8 @@ int    mp_cmp_int(mp_int *a, long z);
 int    mp_isodd(mp_int *a);
 int    mp_iseven(mp_int *a);
 
+unsigned long mp_hash(mp_int *a);
+
 /*------------------------------------------------------------------------*/
 /* Number theoretic                                                       */