Damien George | 438c88d | 2014-02-22 19:25:23 +0000 | [diff] [blame] | 1 | typedef uint16_t mpz_dig_t; |
| 2 | typedef uint32_t mpz_dbl_dig_t; |
| 3 | typedef int32_t mpz_dbl_dig_signed_t; |
| 4 | |
| 5 | typedef struct _mpz_t { |
Damien George | 5104775 | 2014-02-26 17:40:52 +0000 | [diff] [blame] | 6 | machine_uint_t neg : 1; |
Damien George | 06201ff | 2014-03-01 19:50:50 +0000 | [diff] [blame] | 7 | machine_uint_t fixed_dig : 1; |
| 8 | machine_uint_t alloc : 30; |
Damien George | 438c88d | 2014-02-22 19:25:23 +0000 | [diff] [blame] | 9 | machine_uint_t len; |
| 10 | mpz_dig_t *dig; |
| 11 | } mpz_t; |
| 12 | |
Damien George | 06201ff | 2014-03-01 19:50:50 +0000 | [diff] [blame] | 13 | #define MPZ_DIG_SIZE (15) // see mpn_div for why this needs to be at most 15 |
| 14 | #define MPZ_NUM_DIG_FOR_INT (sizeof(machine_int_t) * 8 / MPZ_DIG_SIZE + 1) |
Damien George | bb4a43f | 2014-03-12 15:36:06 +0000 | [diff] [blame] | 15 | #define MPZ_NUM_DIG_FOR_LL (sizeof(long long) * 8 / MPZ_DIG_SIZE + 1) |
Damien George | 06201ff | 2014-03-01 19:50:50 +0000 | [diff] [blame] | 16 | |
| 17 | // convenience macro to declare an mpz with a digit array from the stack, initialised by an integer |
| 18 | #define MPZ_CONST_INT(z, val) mpz_t z; mpz_dig_t z ## _digits[MPZ_NUM_DIG_FOR_INT]; mpz_init_fixed_from_int(&z, z_digits, MPZ_NUM_DIG_FOR_INT, val); |
Damien George | 438c88d | 2014-02-22 19:25:23 +0000 | [diff] [blame] | 19 | |
| 20 | void mpz_init_zero(mpz_t *z); |
| 21 | void mpz_init_from_int(mpz_t *z, machine_int_t val); |
Damien George | 06201ff | 2014-03-01 19:50:50 +0000 | [diff] [blame] | 22 | void mpz_init_fixed_from_int(mpz_t *z, mpz_dig_t *dig, uint dig_alloc, machine_int_t val); |
Damien George | 438c88d | 2014-02-22 19:25:23 +0000 | [diff] [blame] | 23 | void mpz_deinit(mpz_t *z); |
| 24 | |
| 25 | mpz_t *mpz_zero(); |
| 26 | mpz_t *mpz_from_int(machine_int_t i); |
Damien George | bb4a43f | 2014-03-12 15:36:06 +0000 | [diff] [blame] | 27 | mpz_t *mpz_from_ll(long long i); |
Damien George | 438c88d | 2014-02-22 19:25:23 +0000 | [diff] [blame] | 28 | mpz_t *mpz_from_str(const char *str, uint len, bool neg, uint base); |
| 29 | void mpz_free(mpz_t *z); |
| 30 | |
| 31 | mpz_t *mpz_clone(const mpz_t *src); |
| 32 | |
| 33 | void mpz_set(mpz_t *dest, const mpz_t *src); |
| 34 | void mpz_set_from_int(mpz_t *z, machine_int_t src); |
Damien George | bb4a43f | 2014-03-12 15:36:06 +0000 | [diff] [blame] | 35 | void mpz_set_from_ll(mpz_t *z, long long i); |
Damien George | 438c88d | 2014-02-22 19:25:23 +0000 | [diff] [blame] | 36 | uint mpz_set_from_str(mpz_t *z, const char *str, uint len, bool neg, uint base); |
| 37 | |
| 38 | bool mpz_is_zero(const mpz_t *z); |
| 39 | bool mpz_is_pos(const mpz_t *z); |
| 40 | bool mpz_is_neg(const mpz_t *z); |
| 41 | bool mpz_is_odd(const mpz_t *z); |
| 42 | bool mpz_is_even(const mpz_t *z); |
| 43 | |
| 44 | int mpz_cmp(const mpz_t *lhs, const mpz_t *rhs); |
Damien George | 438c88d | 2014-02-22 19:25:23 +0000 | [diff] [blame] | 45 | |
| 46 | mpz_t *mpz_abs(const mpz_t *z); |
| 47 | mpz_t *mpz_neg(const mpz_t *z); |
| 48 | mpz_t *mpz_add(const mpz_t *lhs, const mpz_t *rhs); |
| 49 | mpz_t *mpz_sub(const mpz_t *lhs, const mpz_t *rhs); |
| 50 | mpz_t *mpz_mul(const mpz_t *lhs, const mpz_t *rhs); |
| 51 | mpz_t *mpz_pow(const mpz_t *lhs, const mpz_t *rhs); |
| 52 | |
| 53 | void mpz_abs_inpl(mpz_t *dest, const mpz_t *z); |
| 54 | void mpz_neg_inpl(mpz_t *dest, const mpz_t *z); |
Damien George | 06201ff | 2014-03-01 19:50:50 +0000 | [diff] [blame] | 55 | void mpz_not_inpl(mpz_t *dest, const mpz_t *z); |
| 56 | void mpz_shl_inpl(mpz_t *dest, const mpz_t *lhs, machine_int_t rhs); |
| 57 | void mpz_shr_inpl(mpz_t *dest, const mpz_t *lhs, machine_int_t rhs); |
Damien George | 438c88d | 2014-02-22 19:25:23 +0000 | [diff] [blame] | 58 | void mpz_add_inpl(mpz_t *dest, const mpz_t *lhs, const mpz_t *rhs); |
| 59 | void mpz_sub_inpl(mpz_t *dest, const mpz_t *lhs, const mpz_t *rhs); |
| 60 | void mpz_mul_inpl(mpz_t *dest, const mpz_t *lhs, const mpz_t *rhs); |
| 61 | void mpz_pow_inpl(mpz_t *dest, const mpz_t *lhs, const mpz_t *rhs); |
Paul Sokolovsky | 57207b8 | 2014-03-23 01:52:36 +0200 | [diff] [blame] | 62 | void mpz_and_inpl(mpz_t *dest, const mpz_t *lhs, const mpz_t *rhs); |
| 63 | void mpz_or_inpl(mpz_t *dest, const mpz_t *lhs, const mpz_t *rhs); |
| 64 | void mpz_xor_inpl(mpz_t *dest, const mpz_t *lhs, const mpz_t *rhs); |
Damien George | 438c88d | 2014-02-22 19:25:23 +0000 | [diff] [blame] | 65 | |
| 66 | mpz_t *mpz_gcd(const mpz_t *z1, const mpz_t *z2); |
| 67 | mpz_t *mpz_lcm(const mpz_t *z1, const mpz_t *z2); |
| 68 | void mpz_divmod(const mpz_t *lhs, const mpz_t *rhs, mpz_t **quo, mpz_t **rem); |
| 69 | void mpz_divmod_inpl(mpz_t *dest_quo, mpz_t *dest_rem, const mpz_t *lhs, const mpz_t *rhs); |
| 70 | mpz_t *mpz_div(const mpz_t *lhs, const mpz_t *rhs); |
| 71 | mpz_t *mpz_mod(const mpz_t *lhs, const mpz_t *rhs); |
| 72 | |
Damien George | aca1412 | 2014-02-24 21:32:52 +0000 | [diff] [blame] | 73 | machine_int_t mpz_as_int(const mpz_t *z); |
Damien George | 8270e38 | 2014-04-03 11:00:54 +0000 | [diff] [blame] | 74 | bool mpz_as_int_checked(const mpz_t *z, machine_int_t *value); |
Damien George | 5260810 | 2014-03-08 15:04:54 +0000 | [diff] [blame] | 75 | #if MICROPY_ENABLE_FLOAT |
| 76 | mp_float_t mpz_as_float(const mpz_t *z); |
| 77 | #endif |
Damien George | 438c88d | 2014-02-22 19:25:23 +0000 | [diff] [blame] | 78 | uint mpz_as_str_size(const mpz_t *z, uint base); |
| 79 | char *mpz_as_str(const mpz_t *z, uint base); |
| 80 | uint mpz_as_str_inpl(const mpz_t *z, uint base, char *str); |