|
1983
|
1 /* Definitions of numeric types for XEmacs using the 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 = MINT
|
|
|
25
|
|
|
26 The MP library does not include support for ratios or bigfloats.
|
|
|
27 */
|
|
|
28
|
|
|
29 #ifndef INCLUDED_number_mp_h_
|
|
|
30 #define INCLUDED_number_mp_h_
|
|
|
31
|
|
|
32 /* BSD MP libraries without MP_PREFIX define a function named pow in mp.h that
|
|
|
33 has a different prototype from the one in math.h. We don't use that
|
|
|
34 function anyway, so we do this for safety purposes. However, this means
|
|
|
35 that number-mp.h must always be included before math.h. */
|
|
|
36 #define pow mp_pow
|
|
|
37 #include <mp.h>
|
|
|
38 #undef pow
|
|
|
39
|
|
|
40 #ifdef MP_PREFIX
|
|
|
41 #define MP_GCD mp_gcd
|
|
|
42 #define MP_ITOM mp_itom
|
|
|
43 #define MP_MADD mp_madd
|
|
|
44 #define MP_MCMP mp_mcmp
|
|
|
45 #define MP_MDIV mp_mdiv
|
|
|
46 #define MP_MFREE mp_mfree
|
|
|
47 #define MP_MSUB mp_msub
|
|
|
48 #define MP_MULT mp_mult
|
|
|
49 #define MP_SDIV mp_sdiv
|
|
|
50 #ifdef HAVE_MP_MOVE
|
|
|
51 #define MP_MOVE(x,y) mp_move (x, y)
|
|
|
52 #else
|
|
|
53 #define MP_MOVE(x,y) mp_madd (x, bignum_zero, y)
|
|
|
54 #endif
|
|
|
55 #else
|
|
|
56 #define MP_GCD gcd
|
|
|
57 #define MP_ITOM itom
|
|
|
58 #define MP_MADD madd
|
|
|
59 #define MP_MCMP mcmp
|
|
|
60 #define MP_MDIV mdiv
|
|
|
61 #define MP_MFREE mfree
|
|
|
62 #define MP_MSUB msub
|
|
|
63 #define MP_MULT mult
|
|
|
64 #define MP_SDIV sdiv
|
|
|
65 #ifdef HAVE_MP_MOVE
|
|
|
66 #define MP_MOVE(x,y) move (x, y)
|
|
|
67 #else
|
|
|
68 #define MP_MOVE(x,y) madd (x, bignum_zero, y)
|
|
|
69 #endif
|
|
|
70 #endif
|
|
|
71
|
|
|
72 typedef MINT *bignum;
|
|
|
73
|
|
|
74 extern void init_number_mp(void);
|
|
|
75
|
|
|
76
|
|
|
77 /********************************* Bignums **********************************/
|
|
|
78
|
|
|
79 #define HAVE_BIGNUM 1
|
|
|
80
|
|
|
81 extern MINT *bignum_zero, *intern_bignum;
|
|
|
82 extern MINT *bignum_min_int, *bignum_max_int, *bignum_max_uint;
|
|
|
83 extern MINT *bignum_min_long, *bignum_max_long, *bignum_max_ulong;
|
|
|
84 extern short div_rem;
|
|
|
85
|
|
|
86 /***** Bignum: basic functions *****/
|
|
|
87 #define bignum_init(b) (b = MP_ITOM (0))
|
|
|
88 #define bignum_fini(b) MP_MFREE (b)
|
|
|
89 #define bignum_hashcode(b) bignum_to_uint (b)
|
|
|
90 #define bignum_sign(b) MP_MCMP (b, bignum_zero)
|
|
|
91 #define bignum_evenp(b) (MP_SDIV (b, 2, intern_bignum, &div_rem), \
|
|
|
92 div_rem == 0)
|
|
|
93 #define bignum_oddp(b) (MP_SDIV (b, 2, intern_bignum, &div_rem), \
|
|
|
94 div_rem != 0)
|
|
|
95
|
|
|
96 /***** Bignum: size *****/
|
|
|
97 #define bignum_fits_int_p(b) (MP_MCMP (b, bignum_min_int) >= 0 && \
|
|
|
98 MP_MCMP (b, bignum_max_int) <= 0)
|
|
|
99 #define bignum_fits_uint_p(b) (MP_MCMP (b, bignum_zero) >= 0 && \
|
|
|
100 MP_MCMP (b, bignum_max_uint) <= 0)
|
|
|
101 #define bignum_fits_long_p(b) (MP_MCMP (b, bignum_min_long) >= 0 && \
|
|
|
102 MP_MCMP (b, bignum_max_long) <= 0)
|
|
|
103 #define bignum_fits_ulong_p(b) (MP_MCMP (b, bignum_zero) >= 0 && \
|
|
|
104 MP_MCMP (b, bignum_max_ulong) <= 0)
|
|
|
105
|
|
|
106 /***** Bignum: conversions *****/
|
|
|
107 extern char *bignum_to_string(bignum, int);
|
|
|
108 extern int bignum_to_int(bignum);
|
|
|
109 extern unsigned int bignum_to_uint(bignum);
|
|
|
110 extern long bignum_to_long(bignum);
|
|
|
111 extern unsigned long bignum_to_ulong(bignum);
|
|
|
112 extern double bignum_to_double(bignum);
|
|
|
113
|
|
|
114 /***** Bignum: converting assignments *****/
|
|
|
115 #define bignum_set(b1, b2) MP_MOVE (b2, b1)
|
|
|
116 extern int bignum_set_string(bignum, const char *, int);
|
|
|
117 extern void bignum_set_long(bignum, long);
|
|
|
118 extern void bignum_set_ulong(bignum, unsigned long);
|
|
|
119 extern void bignum_set_double(bignum, double);
|
|
|
120
|
|
|
121 /***** Bignum: comparisons *****/
|
|
|
122 #define bignum_cmp(b1,b2) MP_MCMP (b1, b2)
|
|
|
123 #define bignum_lt(b1,b2) (MP_MCMP (b1, b2) < 0)
|
|
|
124 #define bignum_le(b1,b2) (MP_MCMP (b1, b2) <= 0)
|
|
|
125 #define bignum_eql(b1,b2) (MP_MCMP (b1, b2) == 0)
|
|
|
126 #define bignum_ge(b1,b2) (MP_MCMP (b1, b2) >= 0)
|
|
|
127 #define bignum_gt(b1,b2) (MP_MCMP (b1, b2) > 0)
|
|
|
128
|
|
|
129 /***** Bignum: arithmetic *****/
|
|
|
130 #define bignum_neg(b,b2) MP_MSUB (bignum_zero, b2, b)
|
|
|
131 #define bignum_abs(b,b2) (MP_MCMP (b2, bignum_zero) < 0 \
|
|
|
132 ? MP_MSUB (bignum_zero, b2, b) \
|
|
|
133 : MP_MADD (bignum_zero, b2, b))
|
|
|
134 #define bignum_add(b,b1,b2) MP_MADD (b1, b2, b)
|
|
|
135 #define bignum_sub(b,b1,b2) MP_MSUB (b1, b2, b)
|
|
|
136 #define bignum_mul(b,b1,b2) MP_MULT (b1, b2, b)
|
|
|
137 extern int bignum_divisible_p(bignum, bignum);
|
|
|
138 #define bignum_div(b,b1,b2) MP_MDIV (b1, b2, b, intern_bignum)
|
|
|
139 extern void bignum_ceil(bignum, bignum, bignum);
|
|
|
140 extern void bignum_floor(bignum, bignum, bignum);
|
|
|
141 #define bignum_mod(b,b1,b2) MP_MDIV (b1, b2, intern_bignum, b)
|
|
|
142 extern void bignum_pow(bignum, bignum, unsigned long);
|
|
|
143 #define bignum_gcd(res,b1,b2) MP_GCD (b1, b2, res)
|
|
|
144 extern void bignum_lcm(bignum, bignum, bignum);
|
|
|
145
|
|
|
146 /***** Bignum: bit manipulations *****/
|
|
|
147 extern void bignum_and(bignum, bignum, bignum);
|
|
|
148 extern void bignum_ior(bignum, bignum, bignum);
|
|
|
149 extern void bignum_xor(bignum, bignum, bignum);
|
|
|
150 extern void bignum_not(bignum, bignum);
|
|
|
151 extern void bignum_setbit(bignum, unsigned long);
|
|
|
152 extern void bignum_clrbit(bignum, unsigned long);
|
|
|
153 extern int bignum_testbit(bignum, unsigned long);
|
|
|
154 extern void bignum_lshift(bignum, bignum, unsigned long);
|
|
|
155 extern void bignum_rshift(bignum, bignum, unsigned long);
|
|
|
156
|
|
|
157 /***** Bignum: random numbers *****/
|
|
|
158 extern void bignum_random_seed(unsigned long);
|
|
|
159 extern void bignum_random(bignum, bignum);
|
|
|
160
|
|
|
161 #endif /* INCLUDED_number_mp_h_ */
|
