| 1 | #include <isl_int.h> |
| 2 | |
| 3 | uint32_t isl_imath_hash(mp_int v, uint32_t hash) |
| 4 | { |
| 5 | unsigned const char *data = (unsigned char *)v->digits; |
| 6 | unsigned const char *end = data + v->used * sizeof(v->digits[0]); |
| 7 | |
| 8 | if (v->sign == 1) |
| 9 | isl_hash_byte(hash, 0xFF); |
| 10 | for (; data < end; ++data) |
| 11 | isl_hash_byte(hash, *data); |
| 12 | return hash; |
| 13 | } |
| 14 | |
| 15 | /* Try a standard conversion that fits into a long. |
| 16 | */ |
| 17 | int isl_imath_fits_slong_p(mp_int op) |
| 18 | { |
| 19 | long out; |
| 20 | mp_result res = mp_int_to_int(z: op, out: &out); |
| 21 | return res == MP_OK; |
| 22 | } |
| 23 | |
| 24 | /* Try a standard conversion that fits into an unsigned long. |
| 25 | */ |
| 26 | int isl_imath_fits_ulong_p(mp_int op) |
| 27 | { |
| 28 | unsigned long out; |
| 29 | mp_result res = mp_int_to_uint(z: op, out: &out); |
| 30 | return res == MP_OK; |
| 31 | } |
| 32 | |
| 33 | void isl_imath_addmul_ui(mp_int rop, mp_int op1, unsigned long op2) |
| 34 | { |
| 35 | mpz_t temp; |
| 36 | mp_int_init(z: &temp); |
| 37 | |
| 38 | mp_int_set_uvalue(z: &temp, uvalue: op2); |
| 39 | mp_int_mul(a: op1, b: &temp, c: &temp); |
| 40 | mp_int_add(a: rop, b: &temp, c: rop); |
| 41 | |
| 42 | mp_int_clear(z: &temp); |
| 43 | } |
| 44 | |
| 45 | void isl_imath_submul_ui(mp_int rop, mp_int op1, unsigned long op2) |
| 46 | { |
| 47 | mpz_t temp; |
| 48 | mp_int_init(z: &temp); |
| 49 | |
| 50 | mp_int_set_uvalue(z: &temp, uvalue: op2); |
| 51 | mp_int_mul(a: op1, b: &temp, c: &temp); |
| 52 | mp_int_sub(a: rop, b: &temp, c: rop); |
| 53 | |
| 54 | mp_int_clear(z: &temp); |
| 55 | } |
| 56 | |
| 57 | /* Compute the division of lhs by a rhs of type unsigned long, rounding towards |
| 58 | * positive infinity (Ceil). |
| 59 | */ |
| 60 | void isl_imath_cdiv_q_ui(mp_int rop, mp_int lhs, unsigned long rhs) |
| 61 | { |
| 62 | mpz_t temp; |
| 63 | mp_int_init(z: &temp); |
| 64 | |
| 65 | mp_int_set_uvalue(z: &temp, uvalue: rhs); |
| 66 | impz_cdiv_q(q: rop, n: lhs, d: &temp); |
| 67 | |
| 68 | mp_int_clear(z: &temp); |
| 69 | } |
| 70 | |
| 71 | /* Compute the division of lhs by a rhs of type unsigned long, rounding towards |
| 72 | * negative infinity (Floor). |
| 73 | */ |
| 74 | void isl_imath_fdiv_q_ui(mp_int rop, mp_int lhs, unsigned long rhs) |
| 75 | { |
| 76 | mpz_t temp; |
| 77 | mp_int_init(z: &temp); |
| 78 | |
| 79 | mp_int_set_uvalue(z: &temp, uvalue: rhs); |
| 80 | impz_fdiv_q(q: rop, n: lhs, d: &temp); |
| 81 | |
| 82 | mp_int_clear(z: &temp); |
| 83 | } |
| 84 | |