| 1 | + | /* mini-gmp, a minimalistic implementation of a GNU GMP subset. |
| 2 | + | |
| 3 | + | Copyright 2011-2015 Free Software Foundation, Inc. |
| 4 | + | |
| 5 | + | This file is part of the GNU MP Library. |
| 6 | + | |
| 7 | + | The GNU MP Library is free software; you can redistribute it and/or modify |
| 8 | + | it under the terms of either: |
| 9 | + | |
| 10 | + | * the GNU Lesser General Public License as published by the Free |
| 11 | + | Software Foundation; either version 3 of the License, or (at your |
| 12 | + | option) any later version. |
| 13 | + | |
| 14 | + | or |
| 15 | + | |
| 16 | + | * the GNU General Public License as published by the Free Software |
| 17 | + | Foundation; either version 2 of the License, or (at your option) any |
| 18 | + | later version. |
| 19 | + | |
| 20 | + | or both in parallel, as here. |
| 21 | + | |
| 22 | + | The GNU MP Library is distributed in the hope that it will be useful, but |
| 23 | + | WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY |
| 24 | + | or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License |
| 25 | + | for more details. |
| 26 | + | |
| 27 | + | You should have received copies of the GNU General Public License and the |
| 28 | + | GNU Lesser General Public License along with the GNU MP Library. If not, |
| 29 | + | see https://www.gnu.org/licenses/. */ |
| 30 | + | |
| 31 | + | /* About mini-gmp: This is a minimal implementation of a subset of the |
| 32 | + | GMP interface. It is intended for inclusion into applications which |
| 33 | + | have modest bignums needs, as a fallback when the real GMP library |
| 34 | + | is not installed. |
| 35 | + | |
| 36 | + | This file defines the public interface. */ |
| 37 | + | |
| 38 | + | #ifndef __MINI_GMP_H__ |
| 39 | + | #define __MINI_GMP_H__ |
| 40 | + | |
| 41 | + | /* For size_t */ |
| 42 | + | #include <stddef.h> |
| 43 | + | /* for CHAR_BIT */ |
| 44 | + | |
| 45 | + | #if defined (__cplusplus) |
| 46 | + | extern "C" { |
| 47 | + | #endif |
| 48 | + | |
| 49 | + | void mp_set_memory_functions(void *(*) (size_t), |
| 50 | + | void *(*) (void *, size_t, size_t), |
| 51 | + | void(*) (void *, size_t)); |
| 52 | + | |
| 53 | + | void mp_get_memory_functions(void *(**) (size_t), |
| 54 | + | void *(**) (void *, size_t, size_t), |
| 55 | + | void(**) (void *, size_t)); |
| 56 | + | |
| 57 | + | typedef unsigned long mp_limb_t; |
| 58 | + | typedef long mp_size_t; |
| 59 | + | typedef unsigned long mp_bitcnt_t; |
| 60 | + | |
| 61 | + | typedef mp_limb_t *mp_ptr; |
| 62 | + | typedef const mp_limb_t *mp_srcptr; |
| 63 | + | |
| 64 | + | typedef struct |
| 65 | + | { |
| 66 | + | int _mp_alloc; /* Number of *limbs* allocated and pointed |
| 67 | + | to by the _mp_d field. */ |
| 68 | + | int _mp_size; /* abs(_mp_size) is the number of limbs the |
| 69 | + | last field points to. If _mp_size is |
| 70 | + | negative this is a negative number. */ |
| 71 | + | mp_limb_t *_mp_d; /* Pointer to the limbs. */ |
| 72 | + | } __mpz_struct; |
| 73 | + | |
| 74 | + | typedef __mpz_struct mpz_t[1]; |
| 75 | + | |
| 76 | + | typedef __mpz_struct *mpz_ptr; |
| 77 | + | typedef const __mpz_struct *mpz_srcptr; |
| 78 | + | |
| 79 | + | extern const int mp_bits_per_limb; |
| 80 | + | |
| 81 | + | void mpn_copyi(mp_ptr, mp_srcptr, mp_size_t); |
| 82 | + | void mpn_copyd(mp_ptr, mp_srcptr, mp_size_t); |
| 83 | + | void mpn_zero(mp_ptr, mp_size_t); |
| 84 | + | |
| 85 | + | int mpn_cmp(mp_srcptr, mp_srcptr, mp_size_t); |
| 86 | + | int mpn_zero_p(mp_srcptr, mp_size_t); |
| 87 | + | |
| 88 | + | mp_limb_t mpn_add_1(mp_ptr, mp_srcptr, mp_size_t, mp_limb_t); |
| 89 | + | mp_limb_t mpn_add_n(mp_ptr, mp_srcptr, mp_srcptr, mp_size_t); |
| 90 | + | mp_limb_t mpn_add(mp_ptr, mp_srcptr, mp_size_t, mp_srcptr, mp_size_t); |
| 91 | + | |
| 92 | + | mp_limb_t mpn_sub_1(mp_ptr, mp_srcptr, mp_size_t, mp_limb_t); |
| 93 | + | mp_limb_t mpn_sub_n(mp_ptr, mp_srcptr, mp_srcptr, mp_size_t); |
| 94 | + | mp_limb_t mpn_sub(mp_ptr, mp_srcptr, mp_size_t, mp_srcptr, mp_size_t); |
| 95 | + | |
| 96 | + | mp_limb_t mpn_mul_1(mp_ptr, mp_srcptr, mp_size_t, mp_limb_t); |
| 97 | + | mp_limb_t mpn_addmul_1(mp_ptr, mp_srcptr, mp_size_t, mp_limb_t); |
| 98 | + | mp_limb_t mpn_submul_1(mp_ptr, mp_srcptr, mp_size_t, mp_limb_t); |
| 99 | + | |
| 100 | + | mp_limb_t mpn_mul(mp_ptr, mp_srcptr, mp_size_t, mp_srcptr, mp_size_t); |
| 101 | + | void mpn_mul_n(mp_ptr, mp_srcptr, mp_srcptr, mp_size_t); |
| 102 | + | void mpn_sqr(mp_ptr, mp_srcptr, mp_size_t); |
| 103 | + | int mpn_perfect_square_p(mp_srcptr, mp_size_t); |
| 104 | + | mp_size_t mpn_sqrtrem(mp_ptr, mp_ptr, mp_srcptr, mp_size_t); |
| 105 | + | |
| 106 | + | mp_limb_t mpn_lshift(mp_ptr, mp_srcptr, mp_size_t, unsigned int); |
| 107 | + | mp_limb_t mpn_rshift(mp_ptr, mp_srcptr, mp_size_t, unsigned int); |
| 108 | + | |
| 109 | + | mp_bitcnt_t mpn_scan0(mp_srcptr, mp_bitcnt_t); |
| 110 | + | mp_bitcnt_t mpn_scan1(mp_srcptr, mp_bitcnt_t); |
| 111 | + | |
| 112 | + | void mpn_com(mp_ptr, mp_srcptr, mp_size_t); |
| 113 | + | mp_limb_t mpn_neg(mp_ptr, mp_srcptr, mp_size_t); |
| 114 | + | |
| 115 | + | mp_bitcnt_t mpn_popcount(mp_srcptr, mp_size_t); |
| 116 | + | |
| 117 | + | mp_limb_t mpn_invert_3by2(mp_limb_t, mp_limb_t); |
| 118 | + | #define mpn_invert_limb(x) mpn_invert_3by2 ((x), 0) |
| 119 | + | |
| 120 | + | size_t mpn_get_str(unsigned char *, int, mp_ptr, mp_size_t); |
| 121 | + | mp_size_t mpn_set_str(mp_ptr, const unsigned char *, size_t, int); |
| 122 | + | |
| 123 | + | void mpz_init(mpz_t); |
| 124 | + | void mpz_init2(mpz_t, mp_bitcnt_t); |
| 125 | + | void mpz_clear(mpz_t); |
| 126 | + | |
| 127 | + | #define mpz_odd_p(z) (((z)->_mp_size != 0) & (int) (z)->_mp_d[0]) |
| 128 | + | #define mpz_even_p(z) (! mpz_odd_p (z)) |
| 129 | + | |
| 130 | + | int mpz_sgn(const mpz_t); |
| 131 | + | int mpz_cmp_si(const mpz_t, long); |
| 132 | + | int mpz_cmp_ui(const mpz_t, unsigned long); |
| 133 | + | int mpz_cmp(const mpz_t, const mpz_t); |
| 134 | + | int mpz_cmpabs_ui(const mpz_t, unsigned long); |
| 135 | + | int mpz_cmpabs(const mpz_t, const mpz_t); |
| 136 | + | int mpz_cmp_d(const mpz_t, double); |
| 137 | + | int mpz_cmpabs_d(const mpz_t, double); |
| 138 | + | |
| 139 | + | void mpz_abs(mpz_t, const mpz_t); |
| 140 | + | void mpz_neg(mpz_t, const mpz_t); |
| 141 | + | void mpz_swap(mpz_t, mpz_t); |
| 142 | + | |
| 143 | + | void mpz_add_ui(mpz_t, const mpz_t, unsigned long); |
| 144 | + | void mpz_add(mpz_t, const mpz_t, const mpz_t); |
| 145 | + | void mpz_sub_ui(mpz_t, const mpz_t, unsigned long); |
| 146 | + | void mpz_ui_sub(mpz_t, unsigned long, const mpz_t); |
| 147 | + | void mpz_sub(mpz_t, const mpz_t, const mpz_t); |
| 148 | + | |
| 149 | + | void mpz_mul_si(mpz_t, const mpz_t, long int); |
| 150 | + | void mpz_mul_ui(mpz_t, const mpz_t, unsigned long int); |
| 151 | + | void mpz_mul(mpz_t, const mpz_t, const mpz_t); |
| 152 | + | void mpz_mul_2exp(mpz_t, const mpz_t, mp_bitcnt_t); |
| 153 | + | void mpz_addmul_ui(mpz_t, const mpz_t, unsigned long int); |
| 154 | + | void mpz_addmul(mpz_t, const mpz_t, const mpz_t); |
| 155 | + | void mpz_submul_ui(mpz_t, const mpz_t, unsigned long int); |
| 156 | + | void mpz_submul(mpz_t, const mpz_t, const mpz_t); |
| 157 | + | |
| 158 | + | void mpz_cdiv_qr(mpz_t, mpz_t, const mpz_t, const mpz_t); |
| 159 | + | void mpz_fdiv_qr(mpz_t, mpz_t, const mpz_t, const mpz_t); |
| 160 | + | void mpz_tdiv_qr(mpz_t, mpz_t, const mpz_t, const mpz_t); |
| 161 | + | void mpz_cdiv_q(mpz_t, const mpz_t, const mpz_t); |
| 162 | + | void mpz_fdiv_q(mpz_t, const mpz_t, const mpz_t); |
| 163 | + | void mpz_tdiv_q(mpz_t, const mpz_t, const mpz_t); |
| 164 | + | void mpz_cdiv_r(mpz_t, const mpz_t, const mpz_t); |
| 165 | + | void mpz_fdiv_r(mpz_t, const mpz_t, const mpz_t); |
| 166 | + | void mpz_tdiv_r(mpz_t, const mpz_t, const mpz_t); |
| 167 | + | |
| 168 | + | void mpz_cdiv_q_2exp(mpz_t, const mpz_t, mp_bitcnt_t); |
| 169 | + | void mpz_fdiv_q_2exp(mpz_t, const mpz_t, mp_bitcnt_t); |
| 170 | + | void mpz_tdiv_q_2exp(mpz_t, const mpz_t, mp_bitcnt_t); |
| 171 | + | void mpz_cdiv_r_2exp(mpz_t, const mpz_t, mp_bitcnt_t); |
| 172 | + | void mpz_fdiv_r_2exp(mpz_t, const mpz_t, mp_bitcnt_t); |
| 173 | + | void mpz_tdiv_r_2exp(mpz_t, const mpz_t, mp_bitcnt_t); |
| 174 | + | |
| 175 | + | void mpz_mod(mpz_t, const mpz_t, const mpz_t); |
| 176 | + | |
| 177 | + | void mpz_divexact(mpz_t, const mpz_t, const mpz_t); |
| 178 | + | |
| 179 | + | int mpz_divisible_p(const mpz_t, const mpz_t); |
| 180 | + | int mpz_congruent_p(const mpz_t, const mpz_t, const mpz_t); |
| 181 | + | |
| 182 | + | unsigned long mpz_cdiv_qr_ui(mpz_t, mpz_t, const mpz_t, unsigned long); |
| 183 | + | unsigned long mpz_fdiv_qr_ui(mpz_t, mpz_t, const mpz_t, unsigned long); |
| 184 | + | unsigned long mpz_tdiv_qr_ui(mpz_t, mpz_t, const mpz_t, unsigned long); |
| 185 | + | unsigned long mpz_cdiv_q_ui(mpz_t, const mpz_t, unsigned long); |
| 186 | + | unsigned long mpz_fdiv_q_ui(mpz_t, const mpz_t, unsigned long); |
| 187 | + | unsigned long mpz_tdiv_q_ui(mpz_t, const mpz_t, unsigned long); |
| 188 | + | unsigned long mpz_cdiv_r_ui(mpz_t, const mpz_t, unsigned long); |
| 189 | + | unsigned long mpz_fdiv_r_ui(mpz_t, const mpz_t, unsigned long); |
| 190 | + | unsigned long mpz_tdiv_r_ui(mpz_t, const mpz_t, unsigned long); |
| 191 | + | unsigned long mpz_cdiv_ui(const mpz_t, unsigned long); |
| 192 | + | unsigned long mpz_fdiv_ui(const mpz_t, unsigned long); |
| 193 | + | unsigned long mpz_tdiv_ui(const mpz_t, unsigned long); |
| 194 | + | |
| 195 | + | unsigned long mpz_mod_ui(mpz_t, const mpz_t, unsigned long); |
| 196 | + | |
| 197 | + | void mpz_divexact_ui(mpz_t, const mpz_t, unsigned long); |
| 198 | + | |
| 199 | + | int mpz_divisible_ui_p(const mpz_t, unsigned long); |
| 200 | + | |
| 201 | + | unsigned long mpz_gcd_ui(mpz_t, const mpz_t, unsigned long); |
| 202 | + | void mpz_gcd(mpz_t, const mpz_t, const mpz_t); |
| 203 | + | void mpz_gcdext(mpz_t, mpz_t, mpz_t, const mpz_t, const mpz_t); |
| 204 | + | void mpz_lcm_ui(mpz_t, const mpz_t, unsigned long); |
| 205 | + | void mpz_lcm(mpz_t, const mpz_t, const mpz_t); |
| 206 | + | int mpz_invert(mpz_t, const mpz_t, const mpz_t); |
| 207 | + | |
| 208 | + | void mpz_sqrtrem(mpz_t, mpz_t, const mpz_t); |
| 209 | + | void mpz_sqrt(mpz_t, const mpz_t); |
| 210 | + | int mpz_perfect_square_p(const mpz_t); |
| 211 | + | |
| 212 | + | void mpz_pow_ui(mpz_t, const mpz_t, unsigned long); |
| 213 | + | void mpz_ui_pow_ui(mpz_t, unsigned long, unsigned long); |
| 214 | + | void mpz_powm(mpz_t, const mpz_t, const mpz_t, const mpz_t); |
| 215 | + | void mpz_powm_ui(mpz_t, const mpz_t, unsigned long, const mpz_t); |
| 216 | + | |
| 217 | + | void mpz_rootrem(mpz_t, mpz_t, const mpz_t, unsigned long); |
| 218 | + | int mpz_root(mpz_t, const mpz_t, unsigned long); |
| 219 | + | |
| 220 | + | void mpz_fac_ui(mpz_t, unsigned long); |
| 221 | + | void mpz_bin_uiui(mpz_t, unsigned long, unsigned long); |
| 222 | + | |
| 223 | + | int mpz_probab_prime_p(const mpz_t, int); |
| 224 | + | |
| 225 | + | int mpz_tstbit(const mpz_t, mp_bitcnt_t); |
| 226 | + | void mpz_setbit(mpz_t, mp_bitcnt_t); |
| 227 | + | void mpz_clrbit(mpz_t, mp_bitcnt_t); |
| 228 | + | void mpz_combit(mpz_t, mp_bitcnt_t); |
| 229 | + | |
| 230 | + | void mpz_com(mpz_t, const mpz_t); |
| 231 | + | void mpz_and(mpz_t, const mpz_t, const mpz_t); |
| 232 | + | void mpz_ior(mpz_t, const mpz_t, const mpz_t); |
| 233 | + | void mpz_xor(mpz_t, const mpz_t, const mpz_t); |
| 234 | + | |
| 235 | + | mp_bitcnt_t mpz_popcount(const mpz_t); |
| 236 | + | mp_bitcnt_t mpz_hamdist(const mpz_t, const mpz_t); |
| 237 | + | mp_bitcnt_t mpz_scan0(const mpz_t, mp_bitcnt_t); |
| 238 | + | mp_bitcnt_t mpz_scan1(const mpz_t, mp_bitcnt_t); |
| 239 | + | |
| 240 | + | int mpz_fits_slong_p(const mpz_t); |
| 241 | + | int mpz_fits_ulong_p(const mpz_t); |
| 242 | + | long int mpz_get_si(const mpz_t); |
| 243 | + | unsigned long int mpz_get_ui(const mpz_t); |
| 244 | + | double mpz_get_d(const mpz_t); |
| 245 | + | size_t mpz_size(const mpz_t); |
| 246 | + | mp_limb_t mpz_getlimbn(const mpz_t, mp_size_t); |
| 247 | + | |
| 248 | + | void mpz_realloc2(mpz_t, mp_bitcnt_t); |
| 249 | + | mp_srcptr mpz_limbs_read(mpz_srcptr); |
| 250 | + | mp_ptr mpz_limbs_modify(mpz_t, mp_size_t); |
| 251 | + | mp_ptr mpz_limbs_write(mpz_t, mp_size_t); |
| 252 | + | void mpz_limbs_finish(mpz_t, mp_size_t); |
| 253 | + | mpz_srcptr mpz_roinit_n(mpz_t, mp_srcptr, mp_size_t); |
| 254 | + | |
| 255 | + | #define MPZ_ROINIT_N(xp, xs) {{0, (xs),(xp) }} |
| 256 | + | |
| 257 | + | void mpz_set_si(mpz_t, signed long int); |
| 258 | + | void mpz_set_ui(mpz_t, unsigned long int); |
| 259 | + | void mpz_set(mpz_t, const mpz_t); |
| 260 | + | void mpz_set_d(mpz_t, double); |
| 261 | + | |
| 262 | + | void mpz_init_set_si(mpz_t, signed long int); |
| 263 | + | void mpz_init_set_ui(mpz_t, unsigned long int); |
| 264 | + | void mpz_init_set(mpz_t, const mpz_t); |
| 265 | + | void mpz_init_set_d(mpz_t, double); |
| 266 | + | |
| 267 | + | size_t mpz_sizeinbase(const mpz_t, int); |
| 268 | + | char *mpz_get_str(char *, int, const mpz_t); |
| 269 | + | int mpz_set_str(mpz_t, const char *, int); |
| 270 | + | int mpz_init_set_str(mpz_t, const char *, int); |
| 271 | + | |
| 272 | + | /* This long list taken from gmp.h. */ |
| 273 | + | /* For reference, "defined(EOF)" cannot be used here. In g++ 2.95.4, |
| 274 | + | <iostream> defines EOF but not FILE. */ |
| 275 | + | #if defined (FILE) \ |
| 276 | + | || defined (H_STDIO) \ |
| 277 | + | || defined (_H_STDIO) /* AIX */ \ |
| 278 | + | || defined (_STDIO_H) /* glibc, Sun, SCO */ \ |
| 279 | + | || defined (_STDIO_H_) /* BSD, OSF */ \ |
| 280 | + | || defined (__STDIO_H) /* Borland */ \ |
| 281 | + | || defined (__STDIO_H__) /* IRIX */ \ |
| 282 | + | || defined (_STDIO_INCLUDED) /* HPUX */ \ |
| 283 | + | || defined (__dj_include_stdio_h_) /* DJGPP */ \ |
| 284 | + | || defined (_FILE_DEFINED) /* Microsoft */ \ |
| 285 | + | || defined (__STDIO__) /* Apple MPW MrC */ \ |
| 286 | + | || defined (_MSL_STDIO_H) /* Metrowerks */ \ |
| 287 | + | || defined (_STDIO_H_INCLUDED) /* QNX4 */ \ |
| 288 | + | || defined (_ISO_STDIO_ISO_H) /* Sun C++ */ \ |
| 289 | + | || defined (__STDIO_LOADED) /* VMS */ |
| 290 | + | size_t mpz_out_str(FILE *, int, const mpz_t); |
| 291 | + | #endif |
| 292 | + | |
| 293 | + | void mpz_import(mpz_t, size_t, int, size_t, int, size_t, const void *); |
| 294 | + | void *mpz_export(void *, size_t *, int, size_t, int, size_t, const mpz_t); |
| 295 | + | |
| 296 | + | #if defined (__cplusplus) |
| 297 | + | } |
| 298 | + | #endif |
| 299 | + | |
| 300 | + | // additional defines that some program uses |
| 301 | + | #ifndef GMP_LIMB_BITS |
| 302 | + | #define GMP_LIMB_BITS (sizeof(mp_limb_t) * CHAR_BIT) |
| 303 | + | #endif |
| 304 | + | #ifndef GMP_NAIL_BITS |
| 305 | + | #define GMP_NAIL_BITS 0 |
| 306 | + | #endif |
| 307 | + | #ifndef GMP_NUMB_BITS |
| 308 | + | #define GMP_NUMB_BITS GMP_LIMB_BITS |
| 309 | + | #endif |
| 310 | + | |
| 311 | + | |
| 312 | + | #endif /* __MINI_GMP_H__ */ |
| 313 | + | |
| 314 | + | |
| 315 | + | |