|
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
|
|
|
18 the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
|
|
|
19 Boston, MA 02111-1307, USA. */
|
|
|
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
|
|
1983
|
38 #include <gmp.h>
|
|
|
39
|
|
|
40 typedef mpz_t bignum;
|
|
|
41 typedef mpq_t ratio;
|
|
|
42 typedef mpf_t bigfloat;
|
|
|
43
|
|
|
44 extern void init_number_gmp(void);
|
|
|
45
|
|
|
46
|
|
|
47 /********************************* Bignums **********************************/
|
|
|
48
|
|
|
49 #define HAVE_BIGNUM 1
|
|
|
50
|
|
|
51 extern gmp_randstate_t random_state;
|
|
|
52
|
|
|
53 /***** Bignum: basic functions *****/
|
|
|
54 #define bignum_init(b) mpz_init (b)
|
|
|
55 #define bignum_fini(b) mpz_clear (b)
|
|
|
56 #define bignum_hashcode(b) mpz_get_ui (b)
|
|
|
57 #define bignum_sign(b) mpz_sgn (b)
|
|
|
58 #define bignum_evenp(b) mpz_even_p (b)
|
|
|
59 #define bignum_oddp(b) mpz_odd_p (b)
|
|
|
60
|
|
|
61 /***** Bignum: size *****/
|
|
|
62 #define bignum_fits_int_p(b) mpz_fits_sint_p (b)
|
|
|
63 #define bignum_fits_uint_p(b) mpz_fits_uint_p (b)
|
|
|
64 #define bignum_fits_long_p(b) mpz_fits_slong_p (b)
|
|
|
65 #define bignum_fits_ulong_p(b) mpz_fits_ulong_p (b)
|
|
|
66
|
|
|
67 /***** Bignum: conversions *****/
|
|
|
68 #define bignum_to_string(b,base) mpz_get_str (NULL, base, b)
|
|
|
69 #define bignum_to_int(b) ((int) mpz_get_si (b))
|
|
|
70 #define bignum_to_uint(b) ((unsigned int) mpz_get_ui (b))
|
|
|
71 #define bignum_to_long(b) mpz_get_si (b)
|
|
|
72 #define bignum_to_ulong(b) mpz_get_ui (b)
|
|
|
73 #define bignum_to_double(b) mpz_get_d (b)
|
|
|
74
|
|
|
75 /***** Bignum: converting assignments *****/
|
|
|
76 #define bignum_set(b1,b2) mpz_set (b1, b2)
|
|
|
77 #define bignum_set_string(b,s,base) mpz_set_str (b, s, base)
|
|
|
78 #define bignum_set_long(b,l) mpz_set_si (b, l)
|
|
|
79 #define bignum_set_ulong(b,l) mpz_set_ui (b, l)
|
|
|
80 #define bignum_set_double(b,f) mpz_set_d (b, f)
|
|
|
81 #define bignum_set_ratio(b,r) mpz_set_q (b, r)
|
|
|
82 #define bignum_set_bigfloat(b,f) mpz_set_f (b, f)
|
|
|
83
|
|
|
84 /***** Bignum: comparisons *****/
|
|
|
85 #define bignum_cmp(b1,b2) mpz_cmp (b1, b2)
|
|
|
86 #define bignum_lt(b1,b2) (mpz_cmp (b1, b2) < 0)
|
|
|
87 #define bignum_le(b1,b2) (mpz_cmp (b1, b2) <= 0)
|
|
|
88 #define bignum_eql(b1,b2) (mpz_cmp (b1, b2) == 0)
|
|
|
89 #define bignum_ge(b1,b2) (mpz_cmp (b1, b2) >= 0)
|
|
|
90 #define bignum_gt(b1,b2) (mpz_cmp (b1, b2) > 0)
|
|
|
91
|
|
|
92 /***** Bignum: arithmetic *****/
|
|
|
93 #define bignum_neg(b,b2) mpz_neg (b, b2)
|
|
|
94 #define bignum_abs(b,b2) mpz_abs (b, b2)
|
|
|
95 #define bignum_add(b,b1,b2) mpz_add (b, b1, b2)
|
|
|
96 #define bignum_sub(b,b1,b2) mpz_sub (b, b1, b2)
|
|
|
97 #define bignum_mul(b,b1,b2) mpz_mul (b, b1, b2)
|
|
|
98 #define bignum_divisible_p(b1,b2) mpz_divisible_p (b1, b2)
|
|
|
99 #define bignum_div(b,b1,b2) mpz_tdiv_q (b, b1, b2)
|
|
|
100 #define bignum_ceil(b,b1,b2) mpz_cdiv_q (b, b1, b2)
|
|
|
101 #define bignum_floor(b,b1,b2) mpz_fdiv_q (b, b1, b2)
|
|
|
102 #define bignum_mod(b,b1,b2) mpz_mod (b, b1, b2)
|
|
|
103 #define bignum_pow(res,b,pow) mpz_pow_ui (res, b, pow)
|
|
|
104 #define bignum_gcd(res,b1,b2) mpz_gcd (res, b1, b2)
|
|
|
105 #define bignum_lcm(res,b1,b2) mpz_lcm (res, b1, b2)
|
|
|
106
|
|
|
107 /***** Bignum: bit manipulations *****/
|
|
|
108 #define bignum_and(res,b1,b2) mpz_and (res, b1, b2)
|
|
|
109 #define bignum_ior(res,b1,b2) mpz_ior (res, b1, b2)
|
|
|
110 #define bignum_xor(res,b1,b2) mpz_xor (res, b1, b2)
|
|
|
111 #define bignum_not(res,b) mpz_com (res, b)
|
|
|
112 #define bignum_setbit(b,bit) mpz_setbit (b, bit)
|
|
|
113 #define bignum_clrbit(b,bit) mpz_clrbit (b, bit)
|
|
|
114 #define bignum_testbit(b,bit) mpz_tstbit (b, bit)
|
|
|
115 #define bignum_lshift(res,b,bits) mpz_mul_2exp (res, b, bits)
|
|
|
116 #define bignum_rshift(res,b,bits) mpz_fdiv_q_2exp (res, b, bits)
|
|
|
117
|
|
|
118 /***** Bignum: random numbers *****/
|
|
|
119 #define bignum_random_seed(seed) gmp_randseed_ui (random_state, seed)
|
|
|
120 #define bignum_random(res,limit) mpz_urandomm (res, random_state, limit)
|
|
|
121
|
|
|
122
|
|
|
123 /********************************** Ratios **********************************/
|
|
|
124
|
|
|
125 #define HAVE_RATIO 1
|
|
|
126
|
|
|
127 /***** Ratio: basic functions *****/
|
|
|
128 #define ratio_init(r) mpq_init (r)
|
|
|
129 #define ratio_fini(r) mpq_clear (r)
|
|
|
130 #define ratio_hashcode(r) \
|
|
|
131 (mpz_get_ui (mpq_denref (r)) * mpz_get_ui (mpq_numref (r)))
|
|
|
132 #define ratio_sign(r) mpq_sgn (r)
|
|
|
133 #define ratio_numerator(r) mpq_numref (r)
|
|
|
134 #define ratio_denominator(r) mpq_denref (r)
|
|
|
135 #define ratio_canonicalize(r) mpq_canonicalize (r)
|
|
|
136
|
|
|
137 /***** Ratio: conversions *****/
|
|
|
138 #define ratio_to_string(r,base) mpq_get_str (NULL, base, r)
|
|
|
139 #define ratio_to_int(r) ((int) (mpq_get_d (r)))
|
|
|
140 #define ratio_to_uint(r) ((unsigned int) (mpq_get_d (r)))
|
|
|
141 #define ratio_to_long(r) ((long) (mpq_get_d (r)))
|
|
|
142 #define ratio_to_ulong(r) ((unsigned long) (mpq_get_d (r)))
|
|
|
143 #define ratio_to_double(r) mpq_get_d (r)
|
|
|
144
|
|
|
145 /***** Ratio: converting assignments *****/
|
|
|
146 #define ratio_set(r1,r2) mpq_set (r1, r2)
|
|
|
147 #define ratio_set_string(r,s,base) mpq_set_str (r, s, base)
|
|
|
148 #define ratio_set_long(r,l) mpq_set_si (r, l, 1UL)
|
|
|
149 #define ratio_set_ulong(r,l) mpq_set_ui (r, l, 1UL)
|
|
|
150 #define ratio_set_double(r,f) mpq_set_d (r, f)
|
|
|
151 #define ratio_set_bignum(r,b) mpq_set_z (r, b)
|
|
|
152 #define ratio_set_bigfloat(r,f) mpq_set_f (r, f)
|
|
|
153 #define ratio_set_long_ulong(r,num,den) mpq_set_si (r, num, den)
|
|
|
154 #define ratio_set_ulong_ulong(r,num,den) mpq_set_ui (r, num, den)
|
|
2551
|
155 /* FIXME: Why does this canonicalize, but the previous 2 don't? */
|
|
1983
|
156 #define ratio_set_bignum_bignum(r,num,den) do { \
|
|
|
157 mpz_set (mpq_numref (r), num); \
|
|
|
158 mpz_set (mpq_denref (r), den); \
|
|
|
159 mpq_canonicalize (r); \
|
|
|
160 } while (0)
|
|
|
161
|
|
|
162 /***** Ratio: comparisons *****/
|
|
|
163 #define ratio_cmp(r1,r2) mpq_cmp (r1, r2)
|
|
|
164 #define ratio_lt(r1,r2) (mpq_cmp (r1, r2) < 0)
|
|
|
165 #define ratio_le(r1,r2) (mpq_cmp (r1, r2) <= 0)
|
|
|
166 #define ratio_eql(r1,r2) mpq_equal (r1, r2)
|
|
|
167 #define ratio_ge(r1,r2) (mpq_cmp (r1, r2) >= 0)
|
|
|
168 #define ratio_gt(r1,r2) (mpq_cmp (r1, r2) > 0)
|
|
|
169
|
|
|
170 /***** Ratio: arithmetic *****/
|
|
|
171 #define ratio_neg(q,q2) mpq_neg (q, q2)
|
|
|
172 #define ratio_abs(q,q2) mpq_abs (q, q2)
|
|
|
173 #define ratio_inv(q,q2) mpq_inv (q, q2)
|
|
|
174 #define ratio_add(res,q1,q2) mpq_add (res, q1, q2)
|
|
|
175 #define ratio_sub(res,q1,q2) mpq_sub (res, q1, q2)
|
|
|
176 #define ratio_mul(res,q1,q2) mpq_mul (res, q1, q2)
|
|
|
177 #define ratio_div(res,q1,q2) mpq_div (res, q1, q2)
|
|
|
178
|
|
|
179
|
|
|
180 /******************************** Bigfloats *********************************/
|
|
|
181
|
|
|
182 #define HAVE_BIGFLOAT 1
|
|
|
183
|
|
|
184 /***** Bigfloat: basic functions *****/
|
|
|
185 #define bigfloat_init(f) mpf_init (f)
|
|
|
186 #define bigfloat_init_prec(f,prec) mpf_init2 (f, prec)
|
|
|
187 #define bigfloat_fini(f) mpf_clear (f)
|
|
|
188 #define bigfloat_hashcode(f) mpf_get_ui (f)
|
|
|
189 #define bigfloat_sign(f) mpf_sgn (f)
|
|
|
190 #define bigfloat_get_prec(f) mpf_get_prec (f)
|
|
|
191 #define bigfloat_set_prec(f, prec) mpf_set_prec (f, prec)
|
|
|
192 #define bigfloat_set_default_prec(prec) mpf_set_default_prec(prec)
|
|
|
193 #define bigfloat_get_default_prec() mpf_get_default_prec ()
|
|
|
194
|
|
|
195 /***** Bigfloat: conversions *****/
|
|
|
196 extern CIbyte *bigfloat_to_string (bigfloat f, int base);
|
|
|
197 #define bigfloat_to_int(f) ((int) mpf_get_si (f))
|
|
|
198 #define bigfloat_to_uint(f) ((unsigned int) mpf_get_ui (f))
|
|
|
199 #define bigfloat_to_long(f) mpf_get_si (f)
|
|
|
200 #define bigfloat_to_ulong(f) mpf_get_ui (f)
|
|
|
201 #define bigfloat_to_double(f) mpf_get_d (f)
|
|
|
202
|
|
|
203 /***** Bigfloat: converting assignments *****/
|
|
|
204 #define bigfloat_set(f1,f2) mpf_set (f1, f2)
|
|
|
205 #define bigfloat_set_string(f,str,base) mpf_set_str (f, str, base)
|
|
|
206 #define bigfloat_set_long(f,l) mpf_set_si (f, l)
|
|
|
207 #define bigfloat_set_ulong(f,l) mpf_set_ui (f, l)
|
|
|
208 #define bigfloat_set_double(d,f) mpf_set_d (d, f)
|
|
|
209 #define bigfloat_set_bignum(f,b) mpf_set_z (f, b)
|
|
|
210 #define bigfloat_set_ratio(f,r) mpf_set_q (f, r)
|
|
|
211
|
|
|
212 /***** Bigfloat: comparisons *****/
|
|
|
213 #define bigfloat_cmp(f1,f2) mpf_cmp (f1, f2)
|
|
|
214 #define bigfloat_lt(f1,f2) (mpf_cmp (f1, f2) < 0)
|
|
|
215 #define bigfloat_le(f1,f2) (mpf_cmp (f1, f2) <= 0)
|
|
|
216 #define bigfloat_eql(f1,f2) (mpf_cmp (f1, f2) == 0)
|
|
|
217 #define bigfloat_eql_bits(f1,f2,bits) mpf_eq (f1, f2, bits)
|
|
|
218 #define bigfloat_ge(f1,f2) (mpf_cmp (f1, f2) >= 0)
|
|
|
219 #define bigfloat_gt(f1,f2) (mpf_cmp (f1, f2) > 0)
|
|
|
220
|
|
|
221 /***** Bigfloat: arithmetic *****/
|
|
|
222 #define bigfloat_neg(f,f2) mpf_neg (f, f2)
|
|
|
223 #define bigfloat_abs(f,f2) mpf_abs (f, f2)
|
|
|
224 #define bigfloat_add(res,f1,f2) mpf_add (res, f1, f2)
|
|
|
225 #define bigfloat_sub(res,f1,f2) mpf_sub (res, f1, f2)
|
|
|
226 #define bigfloat_mul(res,f1,f2) mpf_mul (res, f1, f2)
|
|
|
227 #define bigfloat_div(res,f1,f2) mpf_div (res, f1, f2)
|
|
|
228 #define bigfloat_ceil(res,f) mpf_ceil (res, f)
|
|
|
229 #define bigfloat_floor(res,f) mpf_floor (res, f)
|
|
|
230 #define bigfloat_trunc(res,f) mpf_trunc (res, f)
|
|
|
231 #define bigfloat_sqrt(res,f) mpf_sqrt (res, f)
|
|
|
232 #define bigfloat_pow(res,f,exp) mpf_pow_ui (res, f, exp)
|
|
|
233
|
|
|
234 #endif /* INCLUDED_number_gmp_h_ */
|
