diff options
Diffstat (limited to 'usr/rsc/gmp/gcc.c')
| -rw-r--r-- | usr/rsc/gmp/gcc.c | 295 |
1 files changed, 295 insertions, 0 deletions
diff --git a/usr/rsc/gmp/gcc.c b/usr/rsc/gmp/gcc.c index 2e1c884e2..a28266df6 100644 --- a/usr/rsc/gmp/gcc.c +++ b/usr/rsc/gmp/gcc.c @@ -7,8 +7,25 @@ #include <string.h> typedef int32_t int32; +typedef uint32_t uint32; +typedef int64_t int64; typedef uint64_t uint64; +typedef struct Slice Slice; +struct Slice +{ + void *data; + uint32 len; + uint32 cap; +}; + +typedef struct String String; +struct String +{ + void *data; + uint32 len; +}; + typedef struct Int Int; struct Int { @@ -42,6 +59,197 @@ gmp_addInt(void *v) } void +gmp_subInt(void *v) +{ + struct { + Int *z; + Int *x; + Int *y; + Int *ret; + } *a = v; + + a->ret = a->z; + mpz_sub(*a->z->mp, *a->x->mp, *a->y->mp); +} + +void +gmp_mulInt(void *v) +{ + struct { + Int *z; + Int *x; + Int *y; + Int *ret; + } *a = v; + + a->ret = a->z; + mpz_mul(*a->z->mp, *a->x->mp, *a->y->mp); +} + +void +gmp_setInt64Int(void *v) +{ + struct { + Int *z; + int64 x; + Int *ret; + } *a = v; + + a->ret = a->z; + mpz_set_si(*a->z->mp, a->x); +} + +void +gmp_int64Int(void *v) +{ + struct { + Int *z; + int64 ret; + } *a = v; + + a->ret = mpz_get_si(*a->z->mp); +} + +void +gmp_divInt(void *v) +{ + struct { + Int *z; + Int *x; + Int *y; + Int *ret; + } *a = v; + + a->ret = a->z; + mpz_div(*a->z->mp, *a->x->mp, *a->y->mp); +} + +void +gmp_modInt(void *v) +{ + struct { + Int *z; + Int *x; + Int *y; + Int *ret; + } *a = v; + + a->ret = a->z; + mpz_mod(*a->z->mp, *a->x->mp, *a->y->mp); +} + +void +gmp_divModInt(void *v) +{ + struct { + Int *d; + Int *m; + Int *x; + Int *y; + } *a = v; + + mpz_tdiv_qr(*a->d->mp, *a->m->mp, *a->x->mp, *a->y->mp); +} + +void +gmp_lshInt(void *v) +{ + struct { + Int *z; + Int *x; + uint32 s; + Int *ret; + } *a = v; + + a->ret = a->z; + mpz_mul_2exp(*a->z->mp, *a->x->mp, a->s); +} + +void +gmp_rshInt(void *v) +{ + struct { + Int *z; + Int *x; + int32 s; + Int *ret; + } *a = v; + + a->ret = a->z; + mpz_div_2exp(*a->z->mp, *a->x->mp, a->s); +} + + +void +gmp_expInt(void *v) +{ + struct { + Int *z; + Int *x; + Int *y; + Int *w; + Int *ret; + } *a = v; + + a->ret = a->z; + mpz_powm(*a->z->mp, *a->x->mp, *a->y->mp, *a->w->mp); +} + +void +gmp_gcdInt(void *v) +{ + struct { + Int *z; + Int *x; + Int *y; + Int *a; + Int *b; + Int *ret; + } *a = v; + + a->ret = a->z; + mpz_gcdext(*a->z->mp, *a->x->mp, *a->y->mp, *a->a->mp, *a->b->mp); +} + +void +gmp_negInt(void *v) +{ + struct { + Int *z; + Int *x; + Int *ret; + } *a = v; + + a->ret = a->z; + mpz_neg(*a->z->mp, *a->x->mp); +} + +void +gmp_absInt(void *v) +{ + struct { + Int *z; + Int *x; + Int *ret; + } *a = v; + + a->ret = a->z; + mpz_abs(*a->z->mp, *a->x->mp); +} + +void +gmp_cmpInt(void *v) +{ + struct { + Int *x; + Int *y; + int32 ret; + } *a = v; + + a->ret = mpz_cmp(*a->x->mp, *a->y->mp); +} + +void gmp_stringInt(void *v) { struct { @@ -52,3 +260,90 @@ gmp_stringInt(void *v) a->p = mpz_get_str(NULL, 10, *a->z->mp); } +void +gmp_setInt(void *v) +{ + struct { + Int *z; + Int *x; + Int *ret; + } *a = v; + + a->ret = a->z; + mpz_set(*a->z->mp, *a->x->mp); +} + +void +gmp_setBytesInt(void *v) +{ + struct { + Int *z; + Slice b; + Int *ret; + } *a = v; + + a->ret = a->z; + mpz_import(*a->z->mp, a->b.len, 1, 1, 1, 0, a->b.data); +} + +void +gmp_lenInt(void *v) +{ + struct { + Int *z; + int32 ret; + } *a = v; + + a->ret = mpz_sizeinbase(*a->z->mp, 2); +} + +void +gmp_bytesInt(void *v) +{ + struct { + Int *z; + Slice b; + } *a = v; + size_t n; + char *p; + + n = (mpz_sizeinbase(*a->z->mp, 2) + 7) >> 3; + p = malloc(n); // TODO: mallocgc + mpz_export(p, &n, 1, 1, 1, 0, *a->z->mp); + a->b.data = p; + a->b.len = n; + a->b.cap = n; +} + +void +gmp_setStringInt(void *v) +{ + struct { + Int *z; + String s; + int32 base; + int32 pad; + int32 ret; + } *a = v; + char *p; + + p = malloc(a->s.len+1); + memmove(p, a->s.data, a->s.len); + p[a->s.len] = 0; + a->ret = mpz_set_str(*a->z->mp, p, a->base); + free(p); +} + +void +gmp_probablyPrimeInt(void *v) +{ + struct { + Int *z; + int32 nreps; + int32 pad; + int32 ret; + } *a = v; + + a->ret = mpz_probab_prime_p(*a->z->mp, a->nreps); +} + |