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