Mercurial > hg > xemacs-beta
annotate src/number-gmp.h @ 5739:a2912073be85
Support bignums with MPIR. Add documentation on the bignum, ratio,
and bigfloat implementations. See xemacs-patches message with ID
<CAHCOHQkytZao7Uk9ggeo1HKKJtN1bqO054X2mPsGYyQFjbHrZA@mail.gmail.com>
and following messages.
| author | Jerry James <james@xemacs.org> |
|---|---|
| date | Wed, 19 Jun 2013 09:30:30 -0600 |
| parents | 3192994c49ca |
| children | 574f0cded429 |
| rev | line source |
|---|---|
|
5739
a2912073be85
Support bignums with MPIR. Add documentation on the bignum, ratio,
Jerry James <james@xemacs.org>
parents:
5736
diff
changeset
|
1 /* Definitions of numeric types for XEmacs using the GMP or MPIR library. |
|
a2912073be85
Support bignums with MPIR. Add documentation on the bignum, ratio,
Jerry James <james@xemacs.org>
parents:
5736
diff
changeset
|
2 Copyright (C) 2004,2013 Jerry James. |
| 1983 | 3 |
| 4 This file is part of XEmacs. | |
| 5 | |
|
5405
2aa9cd456ae7
Move src/ to GPLv3.
Mike Sperber <sperber@deinprogramm.de>
parents:
5231
diff
changeset
|
6 XEmacs is free software: you can redistribute it and/or modify it |
| 1983 | 7 under the terms of the GNU General Public License as published by the |
|
5405
2aa9cd456ae7
Move src/ to GPLv3.
Mike Sperber <sperber@deinprogramm.de>
parents:
5231
diff
changeset
|
8 Free Software Foundation, either version 3 of the License, or (at your |
|
2aa9cd456ae7
Move src/ to GPLv3.
Mike Sperber <sperber@deinprogramm.de>
parents:
5231
diff
changeset
|
9 option) any later version. |
| 1983 | 10 |
| 11 XEmacs is distributed in the hope that it will be useful, but WITHOUT | |
| 12 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
| 13 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |
| 14 for more details. | |
| 15 | |
| 16 You should have received a copy of the GNU General Public License | |
|
5405
2aa9cd456ae7
Move src/ to GPLv3.
Mike Sperber <sperber@deinprogramm.de>
parents:
5231
diff
changeset
|
17 along with XEmacs. If not, see <http://www.gnu.org/licenses/>. */ |
| 1983 | 18 |
| 19 /* Synched up with: Not in FSF. */ | |
| 20 | |
| 21 /* This library defines the following types: | |
| 22 bignum = mpz_t | |
| 23 ratio = mpq_t | |
| 24 bigfloat = mpf_t | |
| 25 */ | |
| 26 | |
| 27 #ifndef INCLUDED_number_gmp_h_ | |
| 28 #define INCLUDED_number_gmp_h_ | |
| 29 | |
| 2286 | 30 #ifdef UNO |
| 31 /* Uno complains about several inline functions that include conditions with | |
| 32 assignments and side effects if we don't do this */ | |
| 33 #undef __GNUC__ | |
| 34 #endif | |
| 35 | |
| 2956 | 36 #ifdef _MSC_VER |
| 37 /* "unary minus operator applied to unsigned type, result still unsigned": | |
| 38 Occurs on line 1596 of gmp.h in version 4.1.4. */ | |
| 39 #pragma warning ( disable : 4146 ) | |
| 40 #endif | |
|
5739
a2912073be85
Support bignums with MPIR. Add documentation on the bignum, ratio,
Jerry James <james@xemacs.org>
parents:
5736
diff
changeset
|
41 #ifdef WITH_GMP |
| 1983 | 42 #include <gmp.h> |
|
5739
a2912073be85
Support bignums with MPIR. Add documentation on the bignum, ratio,
Jerry James <james@xemacs.org>
parents:
5736
diff
changeset
|
43 #endif |
|
a2912073be85
Support bignums with MPIR. Add documentation on the bignum, ratio,
Jerry James <james@xemacs.org>
parents:
5736
diff
changeset
|
44 #ifdef WITH_MPIR |
|
a2912073be85
Support bignums with MPIR. Add documentation on the bignum, ratio,
Jerry James <james@xemacs.org>
parents:
5736
diff
changeset
|
45 #include <mpir.h> |
|
a2912073be85
Support bignums with MPIR. Add documentation on the bignum, ratio,
Jerry James <james@xemacs.org>
parents:
5736
diff
changeset
|
46 #endif |
| 2956 | 47 #ifdef _MSC_VER |
| 48 #pragma warning ( default : 4146 ) | |
| 49 #endif | |
| 1983 | 50 |
| 51 typedef mpz_t bignum; | |
| 52 typedef mpq_t ratio; | |
| 53 typedef mpf_t bigfloat; | |
| 54 | |
| 55 extern void init_number_gmp(void); | |
| 56 | |
| 57 | |
| 58 /********************************* Bignums **********************************/ | |
| 59 | |
| 60 #define HAVE_BIGNUM 1 | |
| 61 | |
| 62 extern gmp_randstate_t random_state; | |
| 63 | |
| 64 /***** Bignum: basic functions *****/ | |
| 65 #define bignum_init(b) mpz_init (b) | |
| 66 #define bignum_fini(b) mpz_clear (b) | |
| 67 #define bignum_hashcode(b) mpz_get_ui (b) | |
| 68 #define bignum_sign(b) mpz_sgn (b) | |
| 69 #define bignum_evenp(b) mpz_even_p (b) | |
| 70 #define bignum_oddp(b) mpz_odd_p (b) | |
| 71 | |
| 72 /***** Bignum: size *****/ | |
| 73 #define bignum_fits_int_p(b) mpz_fits_sint_p (b) | |
| 74 #define bignum_fits_uint_p(b) mpz_fits_uint_p (b) | |
| 75 #define bignum_fits_long_p(b) mpz_fits_slong_p (b) | |
| 76 #define bignum_fits_ulong_p(b) mpz_fits_ulong_p (b) | |
|
5736
3192994c49ca
Convert C (un)signed long long values to bignums properly.
Jerry James <james@xemacs.org>
parents:
5405
diff
changeset
|
77 #define bignum_fits_llong_p(b) \ |
|
3192994c49ca
Convert C (un)signed long long values to bignums properly.
Jerry James <james@xemacs.org>
parents:
5405
diff
changeset
|
78 (mpz_sizeinbase (b, 2) <= (sizeof(long long) << 3) - 1U) |
|
3192994c49ca
Convert C (un)signed long long values to bignums properly.
Jerry James <james@xemacs.org>
parents:
5405
diff
changeset
|
79 #define bignum_fits_ullong_p(b) \ |
|
3192994c49ca
Convert C (un)signed long long values to bignums properly.
Jerry James <james@xemacs.org>
parents:
5405
diff
changeset
|
80 (mpz_sgn (b) >= 0 && \ |
|
3192994c49ca
Convert C (un)signed long long values to bignums properly.
Jerry James <james@xemacs.org>
parents:
5405
diff
changeset
|
81 mpz_sizeinbase (b, 2) <= (sizeof(unsigned long long) << 3)) |
| 1983 | 82 |
| 83 /***** Bignum: conversions *****/ | |
| 84 #define bignum_to_string(b,base) mpz_get_str (NULL, base, b) | |
| 85 #define bignum_to_int(b) ((int) mpz_get_si (b)) | |
| 86 #define bignum_to_uint(b) ((unsigned int) mpz_get_ui (b)) | |
| 87 #define bignum_to_long(b) mpz_get_si (b) | |
| 88 #define bignum_to_ulong(b) mpz_get_ui (b) | |
|
5736
3192994c49ca
Convert C (un)signed long long values to bignums properly.
Jerry James <james@xemacs.org>
parents:
5405
diff
changeset
|
89 extern long long bignum_to_llong(const bignum b); |
|
3192994c49ca
Convert C (un)signed long long values to bignums properly.
Jerry James <james@xemacs.org>
parents:
5405
diff
changeset
|
90 extern unsigned long long bignum_to_ullong(const bignum b); |
| 1983 | 91 #define bignum_to_double(b) mpz_get_d (b) |
| 92 | |
| 93 /***** Bignum: converting assignments *****/ | |
| 94 #define bignum_set(b1,b2) mpz_set (b1, b2) | |
| 95 #define bignum_set_string(b,s,base) mpz_set_str (b, s, base) | |
| 96 #define bignum_set_long(b,l) mpz_set_si (b, l) | |
| 97 #define bignum_set_ulong(b,l) mpz_set_ui (b, l) | |
|
5736
3192994c49ca
Convert C (un)signed long long values to bignums properly.
Jerry James <james@xemacs.org>
parents:
5405
diff
changeset
|
98 extern void bignum_set_llong(bignum b, long long l); |
|
3192994c49ca
Convert C (un)signed long long values to bignums properly.
Jerry James <james@xemacs.org>
parents:
5405
diff
changeset
|
99 #define bignum_set_ullong(b,l) mpz_import (b,1U,1,sizeof (l),0,0U,&l) |
| 1983 | 100 #define bignum_set_double(b,f) mpz_set_d (b, f) |
| 101 #define bignum_set_ratio(b,r) mpz_set_q (b, r) | |
| 102 #define bignum_set_bigfloat(b,f) mpz_set_f (b, f) | |
| 103 | |
| 104 /***** Bignum: comparisons *****/ | |
| 105 #define bignum_cmp(b1,b2) mpz_cmp (b1, b2) | |
| 106 #define bignum_lt(b1,b2) (mpz_cmp (b1, b2) < 0) | |
| 107 #define bignum_le(b1,b2) (mpz_cmp (b1, b2) <= 0) | |
| 108 #define bignum_eql(b1,b2) (mpz_cmp (b1, b2) == 0) | |
| 109 #define bignum_ge(b1,b2) (mpz_cmp (b1, b2) >= 0) | |
| 110 #define bignum_gt(b1,b2) (mpz_cmp (b1, b2) > 0) | |
| 111 | |
| 112 /***** Bignum: arithmetic *****/ | |
| 113 #define bignum_neg(b,b2) mpz_neg (b, b2) | |
| 114 #define bignum_abs(b,b2) mpz_abs (b, b2) | |
| 115 #define bignum_add(b,b1,b2) mpz_add (b, b1, b2) | |
| 116 #define bignum_sub(b,b1,b2) mpz_sub (b, b1, b2) | |
| 117 #define bignum_mul(b,b1,b2) mpz_mul (b, b1, b2) | |
| 118 #define bignum_divisible_p(b1,b2) mpz_divisible_p (b1, b2) | |
| 119 #define bignum_div(b,b1,b2) mpz_tdiv_q (b, b1, b2) | |
| 120 #define bignum_ceil(b,b1,b2) mpz_cdiv_q (b, b1, b2) | |
| 121 #define bignum_floor(b,b1,b2) mpz_fdiv_q (b, b1, b2) | |
| 122 #define bignum_mod(b,b1,b2) mpz_mod (b, b1, b2) | |
| 123 #define bignum_pow(res,b,pow) mpz_pow_ui (res, b, pow) | |
| 124 #define bignum_gcd(res,b1,b2) mpz_gcd (res, b1, b2) | |
| 125 #define bignum_lcm(res,b1,b2) mpz_lcm (res, b1, b2) | |
| 126 | |
| 127 /***** Bignum: bit manipulations *****/ | |
| 128 #define bignum_and(res,b1,b2) mpz_and (res, b1, b2) | |
| 129 #define bignum_ior(res,b1,b2) mpz_ior (res, b1, b2) | |
| 130 #define bignum_xor(res,b1,b2) mpz_xor (res, b1, b2) | |
| 131 #define bignum_not(res,b) mpz_com (res, b) | |
| 132 #define bignum_setbit(b,bit) mpz_setbit (b, bit) | |
| 133 #define bignum_clrbit(b,bit) mpz_clrbit (b, bit) | |
| 134 #define bignum_testbit(b,bit) mpz_tstbit (b, bit) | |
| 135 #define bignum_lshift(res,b,bits) mpz_mul_2exp (res, b, bits) | |
| 136 #define bignum_rshift(res,b,bits) mpz_fdiv_q_2exp (res, b, bits) | |
| 137 | |
| 138 /***** Bignum: random numbers *****/ | |
| 139 #define bignum_random_seed(seed) gmp_randseed_ui (random_state, seed) | |
| 140 #define bignum_random(res,limit) mpz_urandomm (res, random_state, limit) | |
| 141 | |
| 142 | |
| 143 /********************************** Ratios **********************************/ | |
| 144 | |
| 145 #define HAVE_RATIO 1 | |
| 146 | |
| 147 /***** Ratio: basic functions *****/ | |
| 148 #define ratio_init(r) mpq_init (r) | |
| 149 #define ratio_fini(r) mpq_clear (r) | |
| 150 #define ratio_hashcode(r) \ | |
| 151 (mpz_get_ui (mpq_denref (r)) * mpz_get_ui (mpq_numref (r))) | |
| 152 #define ratio_sign(r) mpq_sgn (r) | |
| 153 #define ratio_numerator(r) mpq_numref (r) | |
| 154 #define ratio_denominator(r) mpq_denref (r) | |
| 155 #define ratio_canonicalize(r) mpq_canonicalize (r) | |
| 156 | |
| 157 /***** Ratio: conversions *****/ | |
| 158 #define ratio_to_string(r,base) mpq_get_str (NULL, base, r) | |
| 159 #define ratio_to_int(r) ((int) (mpq_get_d (r))) | |
| 160 #define ratio_to_uint(r) ((unsigned int) (mpq_get_d (r))) | |
| 161 #define ratio_to_long(r) ((long) (mpq_get_d (r))) | |
| 162 #define ratio_to_ulong(r) ((unsigned long) (mpq_get_d (r))) | |
| 163 #define ratio_to_double(r) mpq_get_d (r) | |
| 164 | |
| 165 /***** Ratio: converting assignments *****/ | |
| 166 #define ratio_set(r1,r2) mpq_set (r1, r2) | |
| 167 #define ratio_set_string(r,s,base) mpq_set_str (r, s, base) | |
| 168 #define ratio_set_long(r,l) mpq_set_si (r, l, 1UL) | |
| 169 #define ratio_set_ulong(r,l) mpq_set_ui (r, l, 1UL) | |
| 170 #define ratio_set_double(r,f) mpq_set_d (r, f) | |
| 171 #define ratio_set_bignum(r,b) mpq_set_z (r, b) | |
| 172 #define ratio_set_bigfloat(r,f) mpq_set_f (r, f) | |
|
5739
a2912073be85
Support bignums with MPIR. Add documentation on the bignum, ratio,
Jerry James <james@xemacs.org>
parents:
5736
diff
changeset
|
173 #define ratio_set_long_ulong(r,num,den) do { \ |
|
a2912073be85
Support bignums with MPIR. Add documentation on the bignum, ratio,
Jerry James <james@xemacs.org>
parents:
5736
diff
changeset
|
174 mpq_set_si (r, num, den); \ |
|
a2912073be85
Support bignums with MPIR. Add documentation on the bignum, ratio,
Jerry James <james@xemacs.org>
parents:
5736
diff
changeset
|
175 mpq_canonicalize (r); \ |
|
a2912073be85
Support bignums with MPIR. Add documentation on the bignum, ratio,
Jerry James <james@xemacs.org>
parents:
5736
diff
changeset
|
176 } while (0) |
|
a2912073be85
Support bignums with MPIR. Add documentation on the bignum, ratio,
Jerry James <james@xemacs.org>
parents:
5736
diff
changeset
|
177 #define ratio_set_ulong_ulong(r,num,den) do { \ |
|
a2912073be85
Support bignums with MPIR. Add documentation on the bignum, ratio,
Jerry James <james@xemacs.org>
parents:
5736
diff
changeset
|
178 mpq_set_ui (r, num, den); \ |
|
a2912073be85
Support bignums with MPIR. Add documentation on the bignum, ratio,
Jerry James <james@xemacs.org>
parents:
5736
diff
changeset
|
179 mpq_canonicalize (r); \ |
|
a2912073be85
Support bignums with MPIR. Add documentation on the bignum, ratio,
Jerry James <james@xemacs.org>
parents:
5736
diff
changeset
|
180 } while (0) |
| 1983 | 181 #define ratio_set_bignum_bignum(r,num,den) do { \ |
| 182 mpz_set (mpq_numref (r), num); \ | |
| 183 mpz_set (mpq_denref (r), den); \ | |
| 184 mpq_canonicalize (r); \ | |
| 185 } while (0) | |
| 186 | |
| 187 /***** Ratio: comparisons *****/ | |
| 188 #define ratio_cmp(r1,r2) mpq_cmp (r1, r2) | |
| 189 #define ratio_lt(r1,r2) (mpq_cmp (r1, r2) < 0) | |
| 190 #define ratio_le(r1,r2) (mpq_cmp (r1, r2) <= 0) | |
| 191 #define ratio_eql(r1,r2) mpq_equal (r1, r2) | |
| 192 #define ratio_ge(r1,r2) (mpq_cmp (r1, r2) >= 0) | |
| 193 #define ratio_gt(r1,r2) (mpq_cmp (r1, r2) > 0) | |
| 194 | |
| 195 /***** Ratio: arithmetic *****/ | |
| 196 #define ratio_neg(q,q2) mpq_neg (q, q2) | |
| 197 #define ratio_abs(q,q2) mpq_abs (q, q2) | |
| 198 #define ratio_inv(q,q2) mpq_inv (q, q2) | |
| 199 #define ratio_add(res,q1,q2) mpq_add (res, q1, q2) | |
| 200 #define ratio_sub(res,q1,q2) mpq_sub (res, q1, q2) | |
| 201 #define ratio_mul(res,q1,q2) mpq_mul (res, q1, q2) | |
| 202 #define ratio_div(res,q1,q2) mpq_div (res, q1, q2) | |
| 203 | |
| 204 | |
| 205 /******************************** Bigfloats *********************************/ | |
| 206 | |
| 207 #define HAVE_BIGFLOAT 1 | |
| 208 | |
| 209 /***** Bigfloat: basic functions *****/ | |
| 210 #define bigfloat_init(f) mpf_init (f) | |
| 211 #define bigfloat_init_prec(f,prec) mpf_init2 (f, prec) | |
| 212 #define bigfloat_fini(f) mpf_clear (f) | |
| 213 #define bigfloat_hashcode(f) mpf_get_ui (f) | |
| 214 #define bigfloat_sign(f) mpf_sgn (f) | |
| 215 #define bigfloat_get_prec(f) mpf_get_prec (f) | |
| 216 #define bigfloat_set_prec(f, prec) mpf_set_prec (f, prec) | |
| 217 #define bigfloat_set_default_prec(prec) mpf_set_default_prec(prec) | |
| 218 #define bigfloat_get_default_prec() mpf_get_default_prec () | |
| 219 | |
| 220 /***** Bigfloat: conversions *****/ | |
| 221 extern CIbyte *bigfloat_to_string (bigfloat f, int base); | |
| 222 #define bigfloat_to_int(f) ((int) mpf_get_si (f)) | |
| 223 #define bigfloat_to_uint(f) ((unsigned int) mpf_get_ui (f)) | |
| 224 #define bigfloat_to_long(f) mpf_get_si (f) | |
| 225 #define bigfloat_to_ulong(f) mpf_get_ui (f) | |
| 226 #define bigfloat_to_double(f) mpf_get_d (f) | |
| 227 | |
| 228 /***** Bigfloat: converting assignments *****/ | |
| 229 #define bigfloat_set(f1,f2) mpf_set (f1, f2) | |
| 230 #define bigfloat_set_string(f,str,base) mpf_set_str (f, str, base) | |
| 231 #define bigfloat_set_long(f,l) mpf_set_si (f, l) | |
| 232 #define bigfloat_set_ulong(f,l) mpf_set_ui (f, l) | |
| 233 #define bigfloat_set_double(d,f) mpf_set_d (d, f) | |
| 234 #define bigfloat_set_bignum(f,b) mpf_set_z (f, b) | |
| 235 #define bigfloat_set_ratio(f,r) mpf_set_q (f, r) | |
| 236 | |
| 237 /***** Bigfloat: comparisons *****/ | |
| 238 #define bigfloat_cmp(f1,f2) mpf_cmp (f1, f2) | |
| 239 #define bigfloat_lt(f1,f2) (mpf_cmp (f1, f2) < 0) | |
| 240 #define bigfloat_le(f1,f2) (mpf_cmp (f1, f2) <= 0) | |
| 241 #define bigfloat_eql(f1,f2) (mpf_cmp (f1, f2) == 0) | |
| 242 #define bigfloat_eql_bits(f1,f2,bits) mpf_eq (f1, f2, bits) | |
| 243 #define bigfloat_ge(f1,f2) (mpf_cmp (f1, f2) >= 0) | |
| 244 #define bigfloat_gt(f1,f2) (mpf_cmp (f1, f2) > 0) | |
| 245 | |
| 246 /***** Bigfloat: arithmetic *****/ | |
| 247 #define bigfloat_neg(f,f2) mpf_neg (f, f2) | |
| 248 #define bigfloat_abs(f,f2) mpf_abs (f, f2) | |
| 249 #define bigfloat_add(res,f1,f2) mpf_add (res, f1, f2) | |
| 250 #define bigfloat_sub(res,f1,f2) mpf_sub (res, f1, f2) | |
| 251 #define bigfloat_mul(res,f1,f2) mpf_mul (res, f1, f2) | |
| 252 #define bigfloat_div(res,f1,f2) mpf_div (res, f1, f2) | |
| 253 #define bigfloat_ceil(res,f) mpf_ceil (res, f) | |
| 254 #define bigfloat_floor(res,f) mpf_floor (res, f) | |
| 255 #define bigfloat_trunc(res,f) mpf_trunc (res, f) | |
| 256 #define bigfloat_sqrt(res,f) mpf_sqrt (res, f) | |
| 257 #define bigfloat_pow(res,f,exp) mpf_pow_ui (res, f, exp) | |
| 258 | |
| 259 #endif /* INCLUDED_number_gmp_h_ */ |