2014-05-03 16:27:38 -06:00
|
|
|
/*
|
|
|
|
* This file is part of the Micro Python project, http://micropython.org/
|
|
|
|
*
|
|
|
|
* The MIT License (MIT)
|
|
|
|
*
|
|
|
|
* Copyright (c) 2013, 2014 Damien P. George
|
|
|
|
*
|
|
|
|
* Permission is hereby granted, free of charge, to any person obtaining a copy
|
|
|
|
* of this software and associated documentation files (the "Software"), to deal
|
|
|
|
* in the Software without restriction, including without limitation the rights
|
|
|
|
* to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
|
|
|
|
* copies of the Software, and to permit persons to whom the Software is
|
|
|
|
* furnished to do so, subject to the following conditions:
|
|
|
|
*
|
|
|
|
* The above copyright notice and this permission notice shall be included in
|
|
|
|
* all copies or substantial portions of the Software.
|
|
|
|
*
|
|
|
|
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
|
|
|
|
* IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
|
|
|
|
* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
|
|
|
|
* AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
|
|
|
|
* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
|
|
|
|
* OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
|
|
|
|
* THE SOFTWARE.
|
|
|
|
*/
|
2015-01-01 13:27:54 -07:00
|
|
|
#ifndef __MICROPY_INCLUDED_PY_MPZ_H__
|
|
|
|
#define __MICROPY_INCLUDED_PY_MPZ_H__
|
|
|
|
|
|
|
|
#include <stdint.h>
|
|
|
|
|
|
|
|
#include "py/mpconfig.h"
|
|
|
|
#include "py/misc.h"
|
2014-05-03 16:27:38 -06:00
|
|
|
|
2014-09-06 10:15:34 -06:00
|
|
|
// This mpz module implements arbitrary precision integers.
|
|
|
|
//
|
|
|
|
// The storage for each digit is defined by mpz_dig_t. The actual number of
|
|
|
|
// bits in mpz_dig_t that are used is defined by MPZ_DIG_SIZE. The machine must
|
|
|
|
// also provide a type that is twice as wide as mpz_dig_t, in both signed and
|
|
|
|
// unsigned versions.
|
|
|
|
//
|
|
|
|
// MPZ_DIG_SIZE can be between 4 and 8*sizeof(mpz_dig_t), but it makes most
|
2015-03-28 15:45:57 -06:00
|
|
|
// sense to have it as large as possible. If MPZ_DIG_SIZE is not already
|
|
|
|
// defined then it is auto-detected below, depending on the machine. The types
|
|
|
|
// are then set based on the value of MPZ_DIG_SIZE (although they can be freely
|
|
|
|
// changed so long as the constraints mentioned above are met).
|
2014-09-06 10:15:34 -06:00
|
|
|
|
2015-03-28 15:45:57 -06:00
|
|
|
#ifndef MPZ_DIG_SIZE
|
|
|
|
#if defined(__x86_64__) || defined(_WIN64)
|
|
|
|
// 64-bit machine, using 32-bit storage for digits
|
|
|
|
#define MPZ_DIG_SIZE (32)
|
|
|
|
#else
|
|
|
|
// default: 32-bit machine, using 16-bit storage for digits
|
|
|
|
#define MPZ_DIG_SIZE (16)
|
|
|
|
#endif
|
|
|
|
#endif
|
|
|
|
|
|
|
|
#if MPZ_DIG_SIZE > 16
|
2014-09-06 10:15:34 -06:00
|
|
|
typedef uint32_t mpz_dig_t;
|
|
|
|
typedef uint64_t mpz_dbl_dig_t;
|
|
|
|
typedef int64_t mpz_dbl_dig_signed_t;
|
2015-03-28 15:45:57 -06:00
|
|
|
#elif MPZ_DIG_SIZE > 8
|
2014-02-22 12:25:23 -07:00
|
|
|
typedef uint16_t mpz_dig_t;
|
|
|
|
typedef uint32_t mpz_dbl_dig_t;
|
|
|
|
typedef int32_t mpz_dbl_dig_signed_t;
|
2015-03-28 15:45:57 -06:00
|
|
|
#elif MPZ_DIG_SIZE > 4
|
|
|
|
typedef uint8_t mpz_dig_t;
|
|
|
|
typedef uint16_t mpz_dbl_dig_t;
|
|
|
|
typedef int16_t mpz_dbl_dig_signed_t;
|
|
|
|
#else
|
|
|
|
typedef uint8_t mpz_dig_t;
|
|
|
|
typedef uint8_t mpz_dbl_dig_t;
|
|
|
|
typedef int8_t mpz_dbl_dig_signed_t;
|
2014-09-06 10:15:34 -06:00
|
|
|
#endif
|
|
|
|
|
2014-10-30 07:39:22 -06:00
|
|
|
#ifdef _WIN64
|
2015-11-20 07:59:06 -07:00
|
|
|
#ifdef __MINGW32__
|
|
|
|
#define MPZ_LONG_1 1LL
|
|
|
|
#else
|
|
|
|
#define MPZ_LONG_1 1i64
|
|
|
|
#endif
|
2014-10-30 07:39:22 -06:00
|
|
|
#else
|
|
|
|
#define MPZ_LONG_1 1L
|
|
|
|
#endif
|
|
|
|
|
2015-09-15 09:15:57 -06:00
|
|
|
// these define the maximum storage needed to hold an int or long long
|
|
|
|
#define MPZ_NUM_DIG_FOR_INT ((sizeof(mp_int_t) * 8 + MPZ_DIG_SIZE - 1) / MPZ_DIG_SIZE)
|
|
|
|
#define MPZ_NUM_DIG_FOR_LL ((sizeof(long long) * 8 + MPZ_DIG_SIZE - 1) / MPZ_DIG_SIZE)
|
2014-02-22 12:25:23 -07:00
|
|
|
|
|
|
|
typedef struct _mpz_t {
|
2017-02-15 21:50:28 -07:00
|
|
|
size_t neg : 1;
|
|
|
|
size_t fixed_dig : 1;
|
|
|
|
size_t alloc : 8 * sizeof(size_t) - 2;
|
|
|
|
size_t len;
|
2014-02-22 12:25:23 -07:00
|
|
|
mpz_dig_t *dig;
|
|
|
|
} mpz_t;
|
|
|
|
|
2014-03-01 12:50:50 -07:00
|
|
|
// convenience macro to declare an mpz with a digit array from the stack, initialised by an integer
|
|
|
|
#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);
|
2014-02-22 12:25:23 -07:00
|
|
|
|
|
|
|
void mpz_init_zero(mpz_t *z);
|
2014-07-03 06:25:24 -06:00
|
|
|
void mpz_init_from_int(mpz_t *z, mp_int_t val);
|
2017-02-15 21:50:28 -07:00
|
|
|
void mpz_init_fixed_from_int(mpz_t *z, mpz_dig_t *dig, size_t dig_alloc, mp_int_t val);
|
2014-02-22 12:25:23 -07:00
|
|
|
void mpz_deinit(mpz_t *z);
|
|
|
|
|
|
|
|
void mpz_set(mpz_t *dest, const mpz_t *src);
|
2014-07-03 06:25:24 -06:00
|
|
|
void mpz_set_from_int(mpz_t *z, mp_int_t src);
|
2014-09-10 15:10:33 -06:00
|
|
|
void mpz_set_from_ll(mpz_t *z, long long i, bool is_signed);
|
2015-01-02 05:39:22 -07:00
|
|
|
#if MICROPY_PY_BUILTINS_FLOAT
|
|
|
|
void mpz_set_from_float(mpz_t *z, mp_float_t src);
|
|
|
|
#endif
|
2017-02-15 21:50:28 -07:00
|
|
|
size_t mpz_set_from_str(mpz_t *z, const char *str, size_t len, bool neg, mp_uint_t base);
|
|
|
|
void mpz_set_from_bytes(mpz_t *z, bool big_endian, size_t len, const byte *buf);
|
2014-02-22 12:25:23 -07:00
|
|
|
|
|
|
|
bool mpz_is_zero(const mpz_t *z);
|
2014-10-03 11:44:14 -06:00
|
|
|
int mpz_cmp(const mpz_t *lhs, const mpz_t *rhs);
|
2014-02-22 12:25:23 -07:00
|
|
|
|
|
|
|
void mpz_abs_inpl(mpz_t *dest, const mpz_t *z);
|
|
|
|
void mpz_neg_inpl(mpz_t *dest, const mpz_t *z);
|
2014-03-01 12:50:50 -07:00
|
|
|
void mpz_not_inpl(mpz_t *dest, const mpz_t *z);
|
2015-10-01 11:01:37 -06:00
|
|
|
void mpz_shl_inpl(mpz_t *dest, const mpz_t *lhs, mp_uint_t rhs);
|
|
|
|
void mpz_shr_inpl(mpz_t *dest, const mpz_t *lhs, mp_uint_t rhs);
|
2014-02-22 12:25:23 -07:00
|
|
|
void mpz_add_inpl(mpz_t *dest, const mpz_t *lhs, const mpz_t *rhs);
|
|
|
|
void mpz_sub_inpl(mpz_t *dest, const mpz_t *lhs, const mpz_t *rhs);
|
|
|
|
void mpz_mul_inpl(mpz_t *dest, const mpz_t *lhs, const mpz_t *rhs);
|
|
|
|
void mpz_pow_inpl(mpz_t *dest, const mpz_t *lhs, const mpz_t *rhs);
|
2017-02-01 16:41:22 -07:00
|
|
|
void mpz_pow3_inpl(mpz_t *dest, const mpz_t *lhs, const mpz_t *rhs, const mpz_t *mod);
|
2014-03-22 17:52:36 -06:00
|
|
|
void mpz_and_inpl(mpz_t *dest, const mpz_t *lhs, const mpz_t *rhs);
|
|
|
|
void mpz_or_inpl(mpz_t *dest, const mpz_t *lhs, const mpz_t *rhs);
|
|
|
|
void mpz_xor_inpl(mpz_t *dest, const mpz_t *lhs, const mpz_t *rhs);
|
2014-02-22 12:25:23 -07:00
|
|
|
void mpz_divmod_inpl(mpz_t *dest_quo, mpz_t *dest_rem, const mpz_t *lhs, const mpz_t *rhs);
|
|
|
|
|
2016-10-10 20:11:32 -06:00
|
|
|
static inline size_t mpz_max_num_bits(const mpz_t *z) { return z->len * MPZ_DIG_SIZE; }
|
2014-07-24 07:21:37 -06:00
|
|
|
mp_int_t mpz_hash(const mpz_t *z);
|
2014-07-03 06:25:24 -06:00
|
|
|
bool mpz_as_int_checked(const mpz_t *z, mp_int_t *value);
|
2014-07-31 07:41:43 -06:00
|
|
|
bool mpz_as_uint_checked(const mpz_t *z, mp_uint_t *value);
|
2017-02-15 21:50:28 -07:00
|
|
|
void mpz_as_bytes(const mpz_t *z, bool big_endian, size_t len, byte *buf);
|
2014-06-01 06:32:54 -06:00
|
|
|
#if MICROPY_PY_BUILTINS_FLOAT
|
2014-03-08 08:04:54 -07:00
|
|
|
mp_float_t mpz_as_float(const mpz_t *z);
|
|
|
|
#endif
|
2017-02-15 21:50:28 -07:00
|
|
|
size_t mpz_as_str_inpl(const mpz_t *z, mp_uint_t base, const char *prefix, char base_char, char comma, char *str);
|
2015-01-01 13:27:54 -07:00
|
|
|
|
|
|
|
#endif // __MICROPY_INCLUDED_PY_MPZ_H__
|