diff options
Diffstat (limited to 'py/objint_mpz.c')
| -rw-r--r-- | py/objint_mpz.c | 181 |
1 files changed, 181 insertions, 0 deletions
diff --git a/py/objint_mpz.c b/py/objint_mpz.c new file mode 100644 index 0000000000..a0889da9e4 --- /dev/null +++ b/py/objint_mpz.c @@ -0,0 +1,181 @@ +#include <stdlib.h> +#include <stdint.h> +#include <string.h> +#include <assert.h> + +#include "nlr.h" +#include "misc.h" +#include "mpconfig.h" +#include "qstr.h" +#include "obj.h" +#include "mpz.h" +#include "objint.h" +#include "runtime0.h" + +#if MICROPY_LONGINT_IMPL == MICROPY_LONGINT_IMPL_MPZ + +STATIC mp_obj_int_t *mp_obj_int_new_mpz(void) { + mp_obj_int_t *o = m_new_obj(mp_obj_int_t); + o->base.type = &int_type; + mpz_init_zero(&o->mpz); + return o; +} + +void int_print(void (*print)(void *env, const char *fmt, ...), void *env, mp_obj_t self_in, mp_print_kind_t kind) { + if (MP_OBJ_IS_SMALL_INT(self_in)) { + print(env, INT_FMT, MP_OBJ_SMALL_INT_VALUE(self_in)); + } else { + // TODO would rather not allocate memory to print... + mp_obj_int_t *self = self_in; + char *str = mpz_as_str(&self->mpz, 10); + print(env, "%s", str); + m_free(str, 0); + } +} + +mp_obj_t int_unary_op(int op, mp_obj_t o_in) { + mp_obj_int_t *o = o_in; + switch (op) { + case RT_UNARY_OP_BOOL: return MP_BOOL(!mpz_is_zero(&o->mpz)); + case RT_UNARY_OP_POSITIVE: return o_in; + case RT_UNARY_OP_NEGATIVE: { mp_obj_int_t *o2 = mp_obj_int_new_mpz(); mpz_neg_inpl(&o2->mpz, &o->mpz); return o2; } + //case RT_UNARY_OP_INVERT: ~ not implemented for mpz + default: return NULL; // op not supported + } +} + +mp_obj_t int_binary_op(int op, mp_obj_t lhs_in, mp_obj_t rhs_in) { + mpz_t *zlhs = &((mp_obj_int_t*)lhs_in)->mpz; + mpz_t *zrhs; + + if (MP_OBJ_IS_SMALL_INT(rhs_in)) { + zrhs = mpz_from_int(MP_OBJ_SMALL_INT_VALUE(rhs_in)); + } else if (MP_OBJ_IS_TYPE(rhs_in, &int_type)) { + zrhs = &((mp_obj_int_t*)rhs_in)->mpz; + } else { + return MP_OBJ_NULL; + } + + if (op == RT_BINARY_OP_TRUE_DIVIDE || op == RT_BINARY_OP_INPLACE_TRUE_DIVIDE) { + machine_float_t flhs = mpz_as_float(zlhs); + machine_float_t frhs = mpz_as_float(zrhs); + return mp_obj_new_float(flhs / frhs); + + } else if (op <= RT_BINARY_OP_POWER) { + mp_obj_int_t *res = mp_obj_int_new_mpz(); + + switch (op) { + case RT_BINARY_OP_ADD: + case RT_BINARY_OP_INPLACE_ADD: + mpz_add_inpl(&res->mpz, zlhs, zrhs); + break; + case RT_BINARY_OP_SUBTRACT: + case RT_BINARY_OP_INPLACE_SUBTRACT: + mpz_sub_inpl(&res->mpz, zlhs, zrhs); + break; + case RT_BINARY_OP_MULTIPLY: + case RT_BINARY_OP_INPLACE_MULTIPLY: + mpz_mul_inpl(&res->mpz, zlhs, zrhs); + break; + case RT_BINARY_OP_FLOOR_DIVIDE: + case RT_BINARY_OP_INPLACE_FLOOR_DIVIDE: { + mpz_t rem; mpz_init_zero(&rem); + mpz_divmod_inpl(&res->mpz, &rem, zlhs, zrhs); + mpz_deinit(&rem); + break; + } + + //case RT_BINARY_OP_MODULO: + //case RT_BINARY_OP_INPLACE_MODULO: + + //case RT_BINARY_OP_AND: + //case RT_BINARY_OP_INPLACE_AND: + //case RT_BINARY_OP_OR: + //case RT_BINARY_OP_INPLACE_OR: + //case RT_BINARY_OP_XOR: + //case RT_BINARY_OP_INPLACE_XOR: + + //case RT_BINARY_OP_LSHIFT: + //case RT_BINARY_OP_INPLACE_LSHIFT: + //case RT_BINARY_OP_RSHIFT: + //case RT_BINARY_OP_INPLACE_RSHIFT: + + case RT_BINARY_OP_POWER: + case RT_BINARY_OP_INPLACE_POWER: + mpz_pow_inpl(&res->mpz, zlhs, zrhs); + break; + + default: + return MP_OBJ_NULL; + } + + return res; + + } else { + int cmp = mpz_cmp(zlhs, zrhs); + switch (op) { + case RT_BINARY_OP_LESS: + return MP_BOOL(cmp < 0); + case RT_BINARY_OP_MORE: + return MP_BOOL(cmp > 0); + case RT_BINARY_OP_LESS_EQUAL: + return MP_BOOL(cmp <= 0); + case RT_BINARY_OP_MORE_EQUAL: + return MP_BOOL(cmp >= 0); + case RT_BINARY_OP_EQUAL: + return MP_BOOL(cmp == 0); + case RT_BINARY_OP_NOT_EQUAL: + return MP_BOOL(cmp != 0); + + default: + return MP_OBJ_NULL; + } + } +} + +mp_obj_t mp_obj_new_int(machine_int_t value) { + if (MP_OBJ_FITS_SMALL_INT(value)) { + return MP_OBJ_NEW_SMALL_INT(value); + } + return mp_obj_new_int_from_ll(value); +} + +mp_obj_t mp_obj_new_int_from_ll(long long val) { + mp_obj_int_t *o = mp_obj_int_new_mpz(); + mpz_set_from_int(&o->mpz, val); + return o; +} + +mp_obj_t mp_obj_new_int_from_uint(machine_uint_t value) { + // SMALL_INT accepts only signed numbers, of one bit less size + // than word size, which totals 2 bits less for unsigned numbers. + if ((value & (WORD_MSBIT_HIGH | (WORD_MSBIT_HIGH >> 1))) == 0) { + return MP_OBJ_NEW_SMALL_INT(value); + } + return mp_obj_new_int_from_ll(value); +} + +mp_obj_t mp_obj_new_int_from_long_str(const char *str) { + mp_obj_int_t *o = mp_obj_int_new_mpz(); + uint len = strlen(str); + uint n = mpz_set_from_str(&o->mpz, str, len, false, 10); + if (n != len) { + nlr_jump(mp_obj_new_exception_msg(&mp_type_SyntaxError, "invalid syntax for number")); + } + return o; +} + +machine_int_t mp_obj_int_get(mp_obj_t self_in) { + if (MP_OBJ_IS_SMALL_INT(self_in)) { + return MP_OBJ_SMALL_INT_VALUE(self_in); + } + mp_obj_int_t *self = self_in; + return mpz_as_int(&self->mpz); +} + +machine_int_t mp_obj_int_get_checked(mp_obj_t self_in) { + // TODO: Check overflow + return mp_obj_int_get(self_in); +} + +#endif |