Mercurial > hg > xemacs-beta
annotate src/number-gmp.h @ 5781:0853e1ec8529
Use alloca_{rawbytes,ibytes} in #'copy-file, #'insert-file-contents-internal
src/ChangeLog addition:
2014-01-20 Aidan Kehoe <kehoea@parhasard.net>
* fileio.c (Fcopy_file, Finsert_file_contents_internal):
Use alloca_{rawbytes,ibytes} here instead of the implicit alloca
on the stack; doesn't change where the buffers are allocated for
these two functions, but does mean that decisions about alloca
vs. malloc based on buffer size are made in the same place
(ultimately, the ALLOCA() macro).
| author | Aidan Kehoe <kehoea@parhasard.net> |
|---|---|
| date | Mon, 20 Jan 2014 17:53:07 +0000 |
| parents | a2912073be85 |
| 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_ */ |