X-Git-Url: https://www.chiark.greenend.org.uk/ucgi/~mdw/git/catacomb-python/blobdiff_plain/bc243788edd0e564df8b67d0750a238d978c954b..75e81374e45c342f15e4eaa658770548c472b15e:/mp.c?ds=sidebyside diff --git a/mp.c b/mp.c index 79c6cf2..62d2dfe 100644 --- a/mp.c +++ b/mp.c @@ -544,13 +544,15 @@ end: return ((PyObject *)zz); } -static long mp_pyhash(PyObject *me) +long mphash(mp *x) { - long h; - PyObject *l = mp_topylong(MP_X(me)); h = PyObject_Hash(l); + PyObject *l = mp_topylong(x); + long h = PyObject_Hash(l); Py_DECREF(l); return (h); } +static long mp_pyhash(PyObject *me) { return (mphash(MP_X(me))); } + static PyObject *mpmeth_jacobi(PyObject *me, PyObject *arg) { mp *y = 0; @@ -682,6 +684,19 @@ end: return (z); } +static PyObject *mpmeth_leastcongruent(PyObject *me, PyObject *arg) +{ + mp *z, *b, *m; + PyObject *rc = 0; + + if (!PyArg_ParseTuple(arg, "O&O&:leastcongruent", convmp, &b, convmp, &m)) + goto end; + z = mp_leastcongruent(MP_NEW, b, MP_X(me), m); + rc = mp_pywrap(z); +end: + return (rc); +} + #define STOREOP(name, c) \ static PyObject *mpmeth_##name(PyObject *me, \ PyObject *arg, PyObject *kw) \ @@ -792,6 +807,8 @@ static PyMethodDef mp_pymethods[] = { "X.gcdx(Y) -> (gcd(X, Y), U, V) with X U + Y V = gcd(X, Y)") METH (modinv, "X.modinv(Y) -> multiplicative inverse of Y mod X") METH (modsqrt, "X.modsqrt(Y) -> square root of Y mod X, if X prime") + METH (leastcongruent, + "X.leastcongruent(B, M) -> smallest Z >= B with Z == X (mod M)") KWMETH(primep, "X.primep(rng = rand) -> true/false if X is prime") KWMETH(tostring, "X.tostring(radix = 10) -> STR") KWMETH(storel, "X.storel(len = -1) -> little-endian bytes") @@ -875,11 +892,15 @@ static PyTypeObject mp_pytype_skel = { /* @tp_doc@ */ "Multiprecision integers, similar to `long' but more efficient and\n\ -versatile. Support all the standard arithmetic operations.\n\ +versatile. Support all the standard arithmetic operations, with\n\ +implicit conversions from `PrimeFilter', and other objects which\n\ +convert to `long'.\n\ \n\ -Constructor mp(X, radix = R) attempts to convert X to an `mp'. If\n\ +Constructor MP(X, radix = R) attempts to convert X to an `MP'. If\n\ X is a string, it's read in radix-R form, or we look for a prefix\n\ -if R = 0. Other acceptable things are ints and longs.\n\ +if R = 0. Other acceptable things are field elements, elliptic curve\n\ +points, group elements, Python `int' and `long' objects, and anything\n\ +with an integer conversion.\n\ \n\ Notes:\n\ \n\ @@ -1870,15 +1891,6 @@ end: return ((PyObject *)zz); } -static long gf_pyhash(PyObject *me) -{ - long i = mp_tolong(MP_X(me)); - i ^= 0xc7ecd67c; /* random perturbance */ - if (i == -1) - i = -2; - return (i); -} - static PyObject *gf_pyexp(PyObject *x, PyObject *y, PyObject *z) { mp *xx = 0, *yy = 0, *zz = 0; @@ -2061,7 +2073,7 @@ static PyTypeObject gf_pytype_skel = { &gf_pynumber, /* @tp_as_number@ */ 0, /* @tp_as_sequence@ */ 0, /* @tp_as_mapping@ */ - gf_pyhash, /* @tp_hash@ */ + mp_pyhash, /* @tp_hash@ */ 0, /* @tp_call@ */ mp_pyhex, /* @tp_str@ */ 0, /* @tp_getattro@ */ @@ -2075,9 +2087,11 @@ static PyTypeObject gf_pytype_skel = { "Binary polynomials. Support almost all the standard arithmetic\n\ operations.\n\ \n\ -Constructor gf(X, radix = R) attempts to convert X to a `gf'. If\n\ +Constructor GF(X, radix = R) attempts to convert X to a `GF'. If\n\ X is a string, it's read in radix-R form, or we look for a prefix\n\ -if R = 0. Other acceptable things are ints and longs.\n\ +if R = 0. Other acceptable things are field elements, elliptic curve\n\ +points, group elements, Python `int' and `long' objects, and anything\n\ +with an integer conversion.\n\ \n\ The name is hopelessly wrong from a technical point of view, but\n\ but it's much easier to type than `p2' or `c2' or whatever.\n\