Commit | Line | Data |
---|---|---|
0f2d19dd JB |
1 | /* classes: h_files */ |
2 | ||
3c9a524f DH |
3 | #ifndef SCM_NUMBERS_H |
4 | #define SCM_NUMBERS_H | |
dee01b01 | 5 | |
7facc08a | 6 | /* Copyright (C) 1995,1996,1998,2000,2001,2002,2003,2004,2005, 2006, |
de7aa61a | 7 | * 2008, 2009, 2010, 2011, 2013, 2014 Free Software Foundation, Inc. |
dee01b01 | 8 | * |
73be1d9e | 9 | * This library is free software; you can redistribute it and/or |
53befeb7 NJ |
10 | * modify it under the terms of the GNU Lesser General Public License |
11 | * as published by the Free Software Foundation; either version 3 of | |
12 | * the License, or (at your option) any later version. | |
dee01b01 | 13 | * |
53befeb7 NJ |
14 | * This library is distributed in the hope that it will be useful, but |
15 | * WITHOUT ANY WARRANTY; without even the implied warranty of | |
73be1d9e MV |
16 | * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
17 | * Lesser General Public License for more details. | |
dee01b01 | 18 | * |
73be1d9e MV |
19 | * You should have received a copy of the GNU Lesser General Public |
20 | * License along with this library; if not, write to the Free Software | |
53befeb7 NJ |
21 | * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA |
22 | * 02110-1301 USA | |
73be1d9e | 23 | */ |
d3a6bc94 | 24 | |
0f2d19dd JB |
25 | \f |
26 | ||
cd036260 MV |
27 | #include <gmp.h> |
28 | ||
b4309c3c | 29 | #include "libguile/__scm.h" |
34f74fe1 MD |
30 | #include "libguile/print.h" |
31 | ||
889975e5 MG |
32 | #ifndef SCM_T_WCHAR_DEFINED |
33 | typedef scm_t_int32 scm_t_wchar; | |
34 | #define SCM_T_WCHAR_DEFINED | |
35 | #endif /* SCM_T_WCHAR_DEFINED */ | |
36 | ||
0f2d19dd JB |
37 | \f |
38 | ||
14282d0f | 39 | /* Immediate Numbers, also known as fixnums |
0f2d19dd | 40 | * |
14282d0f | 41 | * Inums are exact integer data that fits within an SCM word. */ |
0f2d19dd | 42 | |
c0060c95 MV |
43 | /* SCM_T_SIGNED_MAX is (- (expt 2 n) 1), |
44 | * SCM_MOST_POSITIVE_FIXNUM should be (- (expt 2 (- n 2)) 1) | |
45 | * which is the same as (/ (- (expt 2 n) 4) 4) | |
46 | */ | |
47 | ||
004c0902 | 48 | #define SCM_I_FIXNUM_BIT (SCM_LONG_BIT - 2) |
c0060c95 | 49 | #define SCM_MOST_POSITIVE_FIXNUM ((SCM_T_SIGNED_BITS_MAX-3)/4) |
004c0902 | 50 | #define SCM_MOST_NEGATIVE_FIXNUM (-SCM_MOST_POSITIVE_FIXNUM-1) |
5c75b29f | 51 | |
e293c94c MW |
52 | /* SCM_SRS (X, Y) is signed right shift, defined as floor (X / 2^Y), |
53 | where Y must be non-negative and less than the width in bits of X. | |
54 | It's common for >> to do this, but the C standards do not specify | |
55 | what happens when X is negative. | |
0f2d19dd | 56 | |
e293c94c MW |
57 | NOTE: X must not perform side effects. */ |
58 | #if (-1 >> 2 == -1) && (-4 >> 2 == -1) && (-5 >> 2 == -2) && (-8 >> 2 == -2) | |
59 | # define SCM_SRS(x, y) ((x) >> (y)) | |
0f2d19dd | 60 | #else |
e293c94c MW |
61 | # define SCM_SRS(x, y) \ |
62 | ((x) < 0 \ | |
63 | ? -1 - (scm_t_signed_bits) (~(scm_t_bits)(x) >> (y)) \ | |
64 | : ((x) >> (y))) | |
65 | #endif | |
66 | ||
0f2d19dd | 67 | |
3aecd364 MW |
68 | /* The first implementation of SCM_I_INUM below depends on behavior that |
69 | is specified by GNU C but not by C standards, namely that when | |
70 | casting to a signed integer of width N, the value is reduced modulo | |
71 | 2^N to be within range of the type. The second implementation below | |
72 | should be portable to all conforming C implementations, but may be | |
73 | less efficient if the compiler is not sufficiently clever. | |
74 | ||
75 | NOTE: X must not perform side effects. */ | |
76 | #ifdef __GNUC__ | |
77 | # define SCM_I_INUM(x) (SCM_SRS ((scm_t_signed_bits) SCM_UNPACK (x), 2)) | |
78 | #else | |
79 | # define SCM_I_INUM(x) \ | |
80 | (SCM_UNPACK (x) > LONG_MAX \ | |
81 | ? -1 - (scm_t_signed_bits) (~SCM_UNPACK (x) >> 2) \ | |
82 | : (scm_t_signed_bits) (SCM_UNPACK (x) >> 2)) | |
83 | #endif | |
0f2d19dd | 84 | |
3aa13a05 MV |
85 | #define SCM_I_INUMP(x) (2 & SCM_UNPACK (x)) |
86 | #define SCM_I_NINUMP(x) (!SCM_I_INUMP (x)) | |
5a572ca1 | 87 | #define SCM_I_MAKINUM(x) \ |
de7aa61a | 88 | (SCM_PACK ((((scm_t_bits) (x)) << 2) + scm_tc2_int)) |
0f2d19dd | 89 | |
894a712b DH |
90 | /* SCM_FIXABLE is true if its long argument can be encoded in an SCM_INUM. */ |
91 | #define SCM_POSFIXABLE(n) ((n) <= SCM_MOST_POSITIVE_FIXNUM) | |
92 | #define SCM_NEGFIXABLE(n) ((n) >= SCM_MOST_NEGATIVE_FIXNUM) | |
14282d0f | 93 | #define SCM_FIXABLE(n) (SCM_POSFIXABLE (n) && SCM_NEGFIXABLE (n)) |
894a712b DH |
94 | |
95 | ||
cff5fa33 MW |
96 | #define SCM_INUM0 (SCM_I_MAKINUM (0)) /* A name for 0 */ |
97 | #define SCM_INUM1 (SCM_I_MAKINUM (1)) /* A name for 1 */ | |
98 | ||
0f2d19dd | 99 | |
c0060c95 | 100 | /* SCM_MAXEXP is the maximum double precision exponent |
88eb6852 MD |
101 | * SCM_FLTMAX is less than or scm_equal the largest single precision float |
102 | */ | |
103 | ||
f0ad7455 RB |
104 | #if SCM_HAVE_STDC_HEADERS |
105 | # ifndef GO32 | |
1e498fbd SJ |
106 | # include <float.h> |
107 | # ifdef __MINGW32__ | |
108 | # define copysign _copysign | |
1e498fbd SJ |
109 | # define finite _finite |
110 | # endif /* __MINGW32__ */ | |
f0ad7455 | 111 | # endif /* ndef GO32 */ |
88eb6852 | 112 | #endif /* def STDC_HEADERS */ |
f0ad7455 | 113 | |
88eb6852 | 114 | #ifdef DBL_MAX_10_EXP |
f0ad7455 | 115 | # define SCM_MAXEXP DBL_MAX_10_EXP |
88eb6852 | 116 | #else |
f0ad7455 | 117 | # define SCM_MAXEXP 308 /* IEEE doubles */ |
88eb6852 | 118 | #endif /* def DBL_MAX_10_EXP */ |
f0ad7455 | 119 | |
88eb6852 | 120 | #ifdef FLT_MAX |
f0ad7455 | 121 | # define SCM_FLTMAX FLT_MAX |
88eb6852 | 122 | #else |
f0ad7455 | 123 | # define SCM_FLTMAX 1e+23 |
88eb6852 MD |
124 | #endif /* def FLT_MAX */ |
125 | ||
0f2d19dd | 126 | |
2881e77b MV |
127 | /* SCM_INTBUFLEN is the maximum number of characters neccessary for |
128 | * the printed or scm_string representation of an scm_t_intmax in | |
129 | * radix 2. The buffer passed to scm_iint2str and scm_iuint2str must | |
130 | * be of this size, for example. | |
0f2d19dd | 131 | */ |
2881e77b | 132 | #define SCM_INTBUFLEN (5 + SCM_CHAR_BIT*sizeof(scm_t_intmax)) |
0f2d19dd JB |
133 | |
134 | \f | |
135 | ||
136 | /* Numbers | |
137 | */ | |
138 | ||
534c55a9 DH |
139 | |
140 | /* Note that scm_tc16_real and scm_tc16_complex are given tc16-codes that only | |
141 | * differ in one bit: This way, checking if an object is an inexact number can | |
142 | * be done quickly (using the TYP16S macro). */ | |
143 | ||
f03314f9 DH |
144 | /* Number subtype 1 to 3 (note the dependency on the predicates SCM_INEXACTP |
145 | * and SCM_NUMP) */ | |
534c55a9 DH |
146 | #define scm_tc16_big (scm_tc7_number + 1 * 256L) |
147 | #define scm_tc16_real (scm_tc7_number + 2 * 256L) | |
148 | #define scm_tc16_complex (scm_tc7_number + 3 * 256L) | |
f92e85f7 | 149 | #define scm_tc16_fraction (scm_tc7_number + 4 * 256L) |
534c55a9 | 150 | |
f03314f9 DH |
151 | #define SCM_INEXACTP(x) \ |
152 | (!SCM_IMP (x) && (0xfeff & SCM_CELL_TYPE (x)) == scm_tc16_real) | |
dc7da0be AW |
153 | #define SCM_REALP(x) (SCM_HAS_TYP16 (x, scm_tc16_real)) |
154 | #define SCM_COMPLEXP(x) (SCM_HAS_TYP16 (x, scm_tc16_complex)) | |
88eb6852 | 155 | |
0aed71aa AW |
156 | #define SCM_REAL_VALUE(x) (((scm_t_double *) SCM2PTR (x))->real) |
157 | #define SCM_COMPLEX_REAL(x) (((scm_t_complex *) SCM2PTR (x))->real) | |
158 | #define SCM_COMPLEX_IMAG(x) (((scm_t_complex *) SCM2PTR (x))->imag) | |
0f2d19dd | 159 | |
d3e7e88d | 160 | /* Each bignum is just an mpz_t stored in a double cell starting at word 1. */ |
fba8fdc0 | 161 | #define SCM_I_BIG_MPZ(x) (*((mpz_t *) (SCM_CELL_OBJECT_LOC((x),1)))) |
dc7da0be | 162 | #define SCM_BIGP(x) (SCM_HAS_TYP16 (x, scm_tc16_big)) |
0f2d19dd | 163 | |
3aa13a05 | 164 | #define SCM_NUMBERP(x) (SCM_I_INUMP(x) || SCM_NUMP(x)) |
dc7da0be | 165 | #define SCM_NUMP(x) (SCM_HAS_TYP7 (x, scm_tc7_number)) |
f92e85f7 | 166 | |
dc7da0be | 167 | #define SCM_FRACTIONP(x) (SCM_HAS_TYP16 (x, scm_tc16_fraction)) |
651f07f8 DH |
168 | #define SCM_FRACTION_NUMERATOR(x) (SCM_CELL_OBJECT_1 (x)) |
169 | #define SCM_FRACTION_DENOMINATOR(x) (SCM_CELL_OBJECT_2 (x)) | |
88eb6852 | 170 | |
0f2d19dd JB |
171 | \f |
172 | ||
92c2555f | 173 | typedef struct scm_t_double |
0f2d19dd JB |
174 | { |
175 | SCM type; | |
88eb6852 MD |
176 | SCM pad; |
177 | double real; | |
92c2555f | 178 | } scm_t_double; |
0f2d19dd | 179 | |
92c2555f | 180 | typedef struct scm_t_complex |
0f2d19dd | 181 | { |
03604fcf LC |
182 | SCM type; |
183 | SCM pad; | |
88eb6852 MD |
184 | double real; |
185 | double imag; | |
92c2555f | 186 | } scm_t_complex; |
0f2d19dd | 187 | |
904a78f1 | 188 | |
0f2d19dd | 189 | \f |
0f2d19dd | 190 | |
33b001fd | 191 | SCM_API SCM scm_exact_p (SCM x); |
022dda69 | 192 | SCM_API int scm_is_exact (SCM x); |
33b001fd MV |
193 | SCM_API SCM scm_odd_p (SCM n); |
194 | SCM_API SCM scm_even_p (SCM n); | |
2519490c MW |
195 | SCM_API SCM scm_finite_p (SCM x); |
196 | SCM_API SCM scm_inf_p (SCM x); | |
197 | SCM_API SCM scm_nan_p (SCM x); | |
232898d9 MV |
198 | SCM_API SCM scm_inf (void); |
199 | SCM_API SCM scm_nan (void); | |
33b001fd MV |
200 | SCM_API SCM scm_abs (SCM x); |
201 | SCM_API SCM scm_quotient (SCM x, SCM y); | |
202 | SCM_API SCM scm_remainder (SCM x, SCM y); | |
203 | SCM_API SCM scm_modulo (SCM x, SCM y); | |
5fbf680b | 204 | SCM_API void scm_euclidean_divide (SCM x, SCM y, SCM *q, SCM *r); |
ff62c168 MW |
205 | SCM_API SCM scm_euclidean_quotient (SCM x, SCM y); |
206 | SCM_API SCM scm_euclidean_remainder (SCM x, SCM y); | |
8f9da340 MW |
207 | SCM_API void scm_floor_divide (SCM x, SCM y, SCM *q, SCM *r); |
208 | SCM_API SCM scm_floor_quotient (SCM x, SCM y); | |
209 | SCM_API SCM scm_floor_remainder (SCM x, SCM y); | |
210 | SCM_API void scm_ceiling_divide (SCM x, SCM y, SCM *q, SCM *r); | |
211 | SCM_API SCM scm_ceiling_quotient (SCM x, SCM y); | |
212 | SCM_API SCM scm_ceiling_remainder (SCM x, SCM y); | |
213 | SCM_API void scm_truncate_divide (SCM x, SCM y, SCM *q, SCM *r); | |
214 | SCM_API SCM scm_truncate_quotient (SCM x, SCM y); | |
215 | SCM_API SCM scm_truncate_remainder (SCM x, SCM y); | |
5fbf680b | 216 | SCM_API void scm_centered_divide (SCM x, SCM y, SCM *q, SCM *r); |
ff62c168 MW |
217 | SCM_API SCM scm_centered_quotient (SCM x, SCM y); |
218 | SCM_API SCM scm_centered_remainder (SCM x, SCM y); | |
8f9da340 MW |
219 | SCM_API void scm_round_divide (SCM x, SCM y, SCM *q, SCM *r); |
220 | SCM_API SCM scm_round_quotient (SCM x, SCM y); | |
221 | SCM_API SCM scm_round_remainder (SCM x, SCM y); | |
33b001fd MV |
222 | SCM_API SCM scm_gcd (SCM x, SCM y); |
223 | SCM_API SCM scm_lcm (SCM n1, SCM n2); | |
224 | SCM_API SCM scm_logand (SCM n1, SCM n2); | |
225 | SCM_API SCM scm_logior (SCM n1, SCM n2); | |
226 | SCM_API SCM scm_logxor (SCM n1, SCM n2); | |
227 | SCM_API SCM scm_logtest (SCM n1, SCM n2); | |
228 | SCM_API SCM scm_logbit_p (SCM n1, SCM n2); | |
229 | SCM_API SCM scm_lognot (SCM n); | |
d885e204 | 230 | SCM_API SCM scm_modulo_expt (SCM n, SCM k, SCM m); |
33b001fd | 231 | SCM_API SCM scm_integer_expt (SCM z1, SCM z2); |
e08a12b5 MW |
232 | SCM_API SCM scm_ash (SCM n, SCM count); |
233 | SCM_API SCM scm_round_ash (SCM n, SCM count); | |
33b001fd MV |
234 | SCM_API SCM scm_bit_extract (SCM n, SCM start, SCM end); |
235 | SCM_API SCM scm_logcount (SCM n); | |
236 | SCM_API SCM scm_integer_length (SCM n); | |
1be6b49c | 237 | |
5fbf680b | 238 | SCM_INTERNAL SCM scm_i_euclidean_divide (SCM x, SCM y); |
8f9da340 MW |
239 | SCM_INTERNAL SCM scm_i_floor_divide (SCM x, SCM y); |
240 | SCM_INTERNAL SCM scm_i_ceiling_divide (SCM x, SCM y); | |
241 | SCM_INTERNAL SCM scm_i_truncate_divide (SCM x, SCM y); | |
5fbf680b | 242 | SCM_INTERNAL SCM scm_i_centered_divide (SCM x, SCM y); |
8f9da340 | 243 | SCM_INTERNAL SCM scm_i_round_divide (SCM x, SCM y); |
5fbf680b | 244 | |
78d3deb1 AW |
245 | SCM_INTERNAL SCM scm_i_gcd (SCM x, SCM y, SCM rest); |
246 | SCM_INTERNAL SCM scm_i_lcm (SCM x, SCM y, SCM rest); | |
247 | SCM_INTERNAL SCM scm_i_logand (SCM x, SCM y, SCM rest); | |
248 | SCM_INTERNAL SCM scm_i_logior (SCM x, SCM y, SCM rest); | |
249 | SCM_INTERNAL SCM scm_i_logxor (SCM x, SCM y, SCM rest); | |
250 | ||
2881e77b MV |
251 | SCM_API size_t scm_iint2str (scm_t_intmax num, int rad, char *p); |
252 | SCM_API size_t scm_iuint2str (scm_t_uintmax num, int rad, char *p); | |
33b001fd MV |
253 | SCM_API SCM scm_number_to_string (SCM x, SCM radix); |
254 | SCM_API int scm_print_real (SCM sexp, SCM port, scm_print_state *pstate); | |
255 | SCM_API int scm_print_complex (SCM sexp, SCM port, scm_print_state *pstate); | |
256 | SCM_API int scm_bigprint (SCM exp, SCM port, scm_print_state *pstate); | |
2b829bbb KR |
257 | SCM_API SCM scm_c_locale_stringn_to_number (const char *mem, size_t len, |
258 | unsigned int radix); | |
3f47e526 | 259 | SCM_INTERNAL SCM scm_i_string_to_number (SCM str, unsigned int radix); |
33b001fd | 260 | SCM_API SCM scm_string_to_number (SCM str, SCM radix); |
33b001fd MV |
261 | SCM_API SCM scm_bigequal (SCM x, SCM y); |
262 | SCM_API SCM scm_real_equalp (SCM x, SCM y); | |
263 | SCM_API SCM scm_complex_equalp (SCM x, SCM y); | |
264 | SCM_API SCM scm_number_p (SCM x); | |
8507ec80 | 265 | SCM_API SCM scm_complex_p (SCM x); |
33b001fd | 266 | SCM_API SCM scm_real_p (SCM x); |
8507ec80 | 267 | SCM_API SCM scm_rational_p (SCM z); |
33b001fd | 268 | SCM_API SCM scm_integer_p (SCM x); |
900a897c | 269 | SCM_API SCM scm_exact_integer_p (SCM x); |
33b001fd | 270 | SCM_API SCM scm_inexact_p (SCM x); |
022dda69 | 271 | SCM_API int scm_is_inexact (SCM x); |
33b001fd MV |
272 | SCM_API SCM scm_num_eq_p (SCM x, SCM y); |
273 | SCM_API SCM scm_less_p (SCM x, SCM y); | |
274 | SCM_API SCM scm_gr_p (SCM x, SCM y); | |
275 | SCM_API SCM scm_leq_p (SCM x, SCM y); | |
276 | SCM_API SCM scm_geq_p (SCM x, SCM y); | |
277 | SCM_API SCM scm_zero_p (SCM z); | |
278 | SCM_API SCM scm_positive_p (SCM x); | |
279 | SCM_API SCM scm_negative_p (SCM x); | |
280 | SCM_API SCM scm_max (SCM x, SCM y); | |
281 | SCM_API SCM scm_min (SCM x, SCM y); | |
282 | SCM_API SCM scm_sum (SCM x, SCM y); | |
40882e3d | 283 | SCM_API SCM scm_oneplus (SCM x); |
33b001fd | 284 | SCM_API SCM scm_difference (SCM x, SCM y); |
40882e3d | 285 | SCM_API SCM scm_oneminus (SCM x); |
33b001fd | 286 | SCM_API SCM scm_product (SCM x, SCM y); |
33b001fd | 287 | SCM_API SCM scm_divide (SCM x, SCM y); |
f92e85f7 MV |
288 | SCM_API SCM scm_floor (SCM x); |
289 | SCM_API SCM scm_ceiling (SCM x); | |
3101f40f MV |
290 | SCM_API double scm_c_truncate (double x); |
291 | SCM_API double scm_c_round (double x); | |
f92e85f7 MV |
292 | SCM_API SCM scm_truncate_number (SCM x); |
293 | SCM_API SCM scm_round_number (SCM x); | |
6fc4d012 | 294 | SCM_API SCM scm_expt (SCM z1, SCM z2); |
ad79736c AW |
295 | SCM_API SCM scm_sin (SCM z); |
296 | SCM_API SCM scm_cos (SCM z); | |
297 | SCM_API SCM scm_tan (SCM z); | |
298 | SCM_API SCM scm_sinh (SCM z); | |
299 | SCM_API SCM scm_cosh (SCM z); | |
300 | SCM_API SCM scm_tanh (SCM z); | |
301 | SCM_API SCM scm_asin (SCM z); | |
302 | SCM_API SCM scm_acos (SCM z); | |
303 | SCM_API SCM scm_atan (SCM x, SCM y); | |
304 | SCM_API SCM scm_sys_asinh (SCM z); | |
305 | SCM_API SCM scm_sys_acosh (SCM z); | |
306 | SCM_API SCM scm_sys_atanh (SCM z); | |
33b001fd MV |
307 | SCM_API SCM scm_make_rectangular (SCM z1, SCM z2); |
308 | SCM_API SCM scm_make_polar (SCM z1, SCM z2); | |
309 | SCM_API SCM scm_real_part (SCM z); | |
310 | SCM_API SCM scm_imag_part (SCM z); | |
311 | SCM_API SCM scm_magnitude (SCM z); | |
312 | SCM_API SCM scm_angle (SCM z); | |
313 | SCM_API SCM scm_exact_to_inexact (SCM z); | |
314 | SCM_API SCM scm_inexact_to_exact (SCM z); | |
315 | SCM_API SCM scm_trunc (SCM x); | |
8ab3d8a0 KR |
316 | SCM_API SCM scm_log (SCM z); |
317 | SCM_API SCM scm_log10 (SCM z); | |
318 | SCM_API SCM scm_exp (SCM z); | |
319 | SCM_API SCM scm_sqrt (SCM z); | |
882c8963 | 320 | SCM_API void scm_exact_integer_sqrt (SCM k, SCM *s, SCM *r); |
33b001fd | 321 | |
78d3deb1 AW |
322 | SCM_INTERNAL SCM scm_i_min (SCM x, SCM y, SCM rest); |
323 | SCM_INTERNAL SCM scm_i_max (SCM x, SCM y, SCM rest); | |
324 | SCM_INTERNAL SCM scm_i_sum (SCM x, SCM y, SCM rest); | |
325 | SCM_INTERNAL SCM scm_i_difference (SCM x, SCM y, SCM rest); | |
326 | SCM_INTERNAL SCM scm_i_product (SCM x, SCM y, SCM rest); | |
327 | SCM_INTERNAL SCM scm_i_divide (SCM x, SCM y, SCM rest); | |
882c8963 | 328 | SCM_INTERNAL SCM scm_i_exact_integer_sqrt (SCM k); |
78d3deb1 | 329 | |
d3e7e88d | 330 | /* bignum internal functions */ |
102dbb6f LC |
331 | SCM_INTERNAL SCM scm_i_mkbig (void); |
332 | SCM_API /* FIXME: not internal */ SCM scm_i_normbig (SCM x); | |
333 | SCM_INTERNAL int scm_i_bigcmp (SCM a, SCM b); | |
334 | SCM_INTERNAL SCM scm_i_dbl2big (double d); | |
335 | SCM_INTERNAL SCM scm_i_dbl2num (double d); | |
336 | SCM_API /* FIXME: not internal */ double scm_i_big2dbl (SCM b); | |
337 | SCM_API /* FIXME: not internal */ SCM scm_i_long2big (long n); | |
338 | SCM_API /* FIXME: not internal */ SCM scm_i_ulong2big (unsigned long n); | |
339 | SCM_API /* FIXME: not internal */ SCM scm_i_clonebig (SCM src_big, int same_sign_p); | |
d3e7e88d | 340 | |
f92e85f7 | 341 | /* ratio functions */ |
f92e85f7 MV |
342 | SCM_API SCM scm_rationalize (SCM x, SCM err); |
343 | SCM_API SCM scm_numerator (SCM z); | |
344 | SCM_API SCM scm_denominator (SCM z); | |
f92e85f7 MV |
345 | |
346 | /* fraction internal functions */ | |
102dbb6f LC |
347 | SCM_INTERNAL double scm_i_fraction2double (SCM z); |
348 | SCM_INTERNAL SCM scm_i_fraction_equalp (SCM x, SCM y); | |
349 | SCM_INTERNAL int scm_i_print_fraction (SCM sexp, SCM port, scm_print_state *pstate); | |
f92e85f7 | 350 | |
b479fe9a | 351 | /* general internal functions */ |
102dbb6f LC |
352 | SCM_INTERNAL void scm_i_print_double (double val, SCM port); |
353 | SCM_INTERNAL void scm_i_print_complex (double real, double imag, SCM port); | |
b479fe9a | 354 | |
d1b5b400 | 355 | /* conversion functions for integers */ |
73e4de09 MV |
356 | |
357 | SCM_API int scm_is_integer (SCM val); | |
900a897c | 358 | SCM_API int scm_is_exact_integer (SCM val); |
73e4de09 MV |
359 | SCM_API int scm_is_signed_integer (SCM val, |
360 | scm_t_intmax min, scm_t_intmax max); | |
361 | SCM_API int scm_is_unsigned_integer (SCM val, | |
362 | scm_t_uintmax min, scm_t_uintmax max); | |
363 | ||
364 | SCM_API SCM scm_from_signed_integer (scm_t_intmax val); | |
365 | SCM_API SCM scm_from_unsigned_integer (scm_t_uintmax val); | |
366 | ||
367 | SCM_API scm_t_intmax scm_to_signed_integer (SCM val, | |
368 | scm_t_intmax min, | |
369 | scm_t_intmax max); | |
370 | SCM_API scm_t_uintmax scm_to_unsigned_integer (SCM val, | |
371 | scm_t_uintmax min, | |
372 | scm_t_uintmax max); | |
373 | ||
bfd7932e MV |
374 | SCM_API scm_t_int8 scm_to_int8 (SCM x); |
375 | SCM_API SCM scm_from_int8 (scm_t_int8 x); | |
376 | ||
377 | SCM_API scm_t_uint8 scm_to_uint8 (SCM x); | |
378 | SCM_API SCM scm_from_uint8 (scm_t_uint8 x); | |
379 | ||
380 | SCM_API scm_t_int16 scm_to_int16 (SCM x); | |
381 | SCM_API SCM scm_from_int16 (scm_t_int16 x); | |
382 | ||
383 | SCM_API scm_t_uint16 scm_to_uint16 (SCM x); | |
384 | SCM_API SCM scm_from_uint16 (scm_t_uint16 x); | |
385 | ||
386 | SCM_API scm_t_int32 scm_to_int32 (SCM x); | |
387 | SCM_API SCM scm_from_int32 (scm_t_int32 x); | |
388 | ||
389 | SCM_API scm_t_uint32 scm_to_uint32 (SCM x); | |
390 | SCM_API SCM scm_from_uint32 (scm_t_uint32 x); | |
391 | ||
904a78f1 MG |
392 | SCM_API scm_t_wchar scm_to_wchar (SCM x); |
393 | SCM_API SCM scm_from_wchar (scm_t_wchar x); | |
394 | ||
bfd7932e MV |
395 | SCM_API scm_t_int64 scm_to_int64 (SCM x); |
396 | SCM_API SCM scm_from_int64 (scm_t_int64 x); | |
397 | ||
398 | SCM_API scm_t_uint64 scm_to_uint64 (SCM x); | |
399 | SCM_API SCM scm_from_uint64 (scm_t_uint64 x); | |
400 | ||
cd036260 MV |
401 | SCM_API void scm_to_mpz (SCM x, mpz_t rop); |
402 | SCM_API SCM scm_from_mpz (mpz_t rop); | |
403 | ||
404 | ||
d1b5b400 MV |
405 | /* The conversion functions for other types are aliased to the |
406 | appropriate ones from above. We pick the right one based on the | |
407 | size of the type. | |
408 | ||
409 | Not each and every possibility is covered by the code below, and | |
410 | while it is trivial to complete the tests, it might be better to | |
411 | just test for the 'sane' possibilities. When one of the tests | |
412 | below fails, chances are good that some silent assumption somewhere | |
413 | else will also fail. | |
414 | */ | |
415 | ||
416 | #if SCM_SIZEOF_CHAR == 1 | |
417 | #define scm_to_schar scm_to_int8 | |
418 | #define scm_from_schar scm_from_int8 | |
419 | #define scm_to_uchar scm_to_uint8 | |
420 | #define scm_from_uchar scm_from_uint8 | |
73e4de09 | 421 | #if CHAR_MIN == 0 |
d1b5b400 MV |
422 | #define scm_to_char scm_to_uint8 |
423 | #define scm_from_char scm_from_uint8 | |
73e4de09 | 424 | #else |
d1b5b400 MV |
425 | #define scm_to_char scm_to_int8 |
426 | #define scm_from_char scm_from_int8 | |
73e4de09 | 427 | #endif |
d1b5b400 MV |
428 | #else |
429 | #error sizeof(char) is not 1. | |
22008a78 MV |
430 | #endif |
431 | ||
d1b5b400 MV |
432 | #if SCM_SIZEOF_SHORT == 1 |
433 | #define scm_to_short scm_to_int8 | |
434 | #define scm_from_short scm_from_int8 | |
435 | #define scm_to_ushort scm_to_uint8 | |
436 | #define scm_from_ushort scm_from_uint8 | |
437 | #else | |
438 | #if SCM_SIZEOF_SHORT == 2 | |
439 | #define scm_to_short scm_to_int16 | |
440 | #define scm_from_short scm_from_int16 | |
441 | #define scm_to_ushort scm_to_uint16 | |
442 | #define scm_from_ushort scm_from_uint16 | |
443 | #else | |
444 | #if SCM_SIZEOF_SHORT == 4 | |
445 | #define scm_to_short scm_to_int32 | |
446 | #define scm_from_short scm_from_int32 | |
447 | #define scm_to_ushort scm_to_uint32 | |
448 | #define scm_from_ushort scm_from_uint32 | |
73e4de09 | 449 | #else |
d1b5b400 MV |
450 | #error sizeof(short) is not 1, 2, or 4. |
451 | #endif | |
452 | #endif | |
73e4de09 MV |
453 | #endif |
454 | ||
d1b5b400 MV |
455 | #if SCM_SIZEOF_INT == 4 |
456 | #define scm_to_int scm_to_int32 | |
457 | #define scm_from_int scm_from_int32 | |
458 | #define scm_to_uint scm_to_uint32 | |
459 | #define scm_from_uint scm_from_uint32 | |
460 | #else | |
461 | #if SCM_SIZEOF_INT == 8 | |
462 | #define scm_to_int scm_to_int64 | |
463 | #define scm_from_int scm_from_int64 | |
464 | #define scm_to_uint scm_to_uint64 | |
465 | #define scm_from_uint scm_from_uint64 | |
466 | #else | |
467 | #error sizeof(int) is not 4 or 8. | |
468 | #endif | |
469 | #endif | |
73e4de09 | 470 | |
d1b5b400 MV |
471 | #if SCM_SIZEOF_LONG == 4 |
472 | #define scm_to_long scm_to_int32 | |
473 | #define scm_from_long scm_from_int32 | |
474 | #define scm_to_ulong scm_to_uint32 | |
475 | #define scm_from_ulong scm_from_uint32 | |
476 | #else | |
477 | #if SCM_SIZEOF_LONG == 8 | |
478 | #define scm_to_long scm_to_int64 | |
479 | #define scm_from_long scm_from_int64 | |
480 | #define scm_to_ulong scm_to_uint64 | |
481 | #define scm_from_ulong scm_from_uint64 | |
482 | #else | |
483 | #error sizeof(long) is not 4 or 8. | |
484 | #endif | |
485 | #endif | |
73e4de09 | 486 | |
d1b5b400 MV |
487 | #if SCM_SIZEOF_INTMAX == 4 |
488 | #define scm_to_intmax scm_to_int32 | |
489 | #define scm_from_intmax scm_from_int32 | |
490 | #define scm_to_uintmax scm_to_uint32 | |
491 | #define scm_from_uintmax scm_from_uint32 | |
492 | #else | |
493 | #if SCM_SIZEOF_INTMAX == 8 | |
494 | #define scm_to_intmax scm_to_int64 | |
495 | #define scm_from_intmax scm_from_int64 | |
496 | #define scm_to_uintmax scm_to_uint64 | |
497 | #define scm_from_uintmax scm_from_uint64 | |
498 | #else | |
55f26379 | 499 | #error sizeof(scm_t_intmax) is not 4 or 8. |
d1b5b400 MV |
500 | #endif |
501 | #endif | |
73e4de09 | 502 | |
d1b5b400 MV |
503 | #if SCM_SIZEOF_LONG_LONG == 0 |
504 | #else | |
505 | #if SCM_SIZEOF_LONG_LONG == 8 | |
506 | #define scm_to_long_long scm_to_int64 | |
507 | #define scm_from_long_long scm_from_int64 | |
508 | #define scm_to_ulong_long scm_to_uint64 | |
509 | #define scm_from_ulong_long scm_from_uint64 | |
510 | #else | |
511 | #error sizeof(long long) is not 8. | |
512 | #endif | |
513 | #endif | |
73e4de09 | 514 | |
d1b5b400 MV |
515 | #if SCM_SIZEOF_SIZE_T == 4 |
516 | #define scm_to_ssize_t scm_to_int32 | |
517 | #define scm_from_ssize_t scm_from_int32 | |
518 | #define scm_to_size_t scm_to_uint32 | |
519 | #define scm_from_size_t scm_from_uint32 | |
520 | #else | |
521 | #if SCM_SIZEOF_SIZE_T == 8 | |
522 | #define scm_to_ssize_t scm_to_int64 | |
523 | #define scm_from_ssize_t scm_from_int64 | |
524 | #define scm_to_size_t scm_to_uint64 | |
525 | #define scm_from_size_t scm_from_uint64 | |
526 | #else | |
527 | #error sizeof(size_t) is not 4 or 8. | |
528 | #endif | |
22008a78 MV |
529 | #endif |
530 | ||
7facc08a MW |
531 | #if SCM_SIZEOF_SCM_T_PTRDIFF == 4 |
532 | #define scm_to_ptrdiff_t scm_to_int32 | |
533 | #define scm_from_ptrdiff_t scm_from_int32 | |
534 | #else | |
535 | #if SCM_SIZEOF_SCM_T_PTRDIFF == 8 | |
536 | #define scm_to_ptrdiff_t scm_to_int64 | |
537 | #define scm_from_ptrdiff_t scm_from_int64 | |
538 | #else | |
539 | #error sizeof(scm_t_ptrdiff) is not 4 or 8. | |
540 | #endif | |
541 | #endif | |
542 | ||
0b3b7369 AW |
543 | #if SCM_SIZEOF_INTPTR_T == 0 |
544 | /* No intptr_t; use size_t functions. */ | |
545 | #define scm_to_intptr_t scm_to_ssize_t | |
546 | #define scm_from_intptr_t scm_from_ssize_t | |
547 | #elif SCM_SIZEOF_INTPTR_T == 4 | |
548 | #define scm_to_intptr_t scm_to_int32 | |
549 | #define scm_from_intptr_t scm_from_int32 | |
550 | #elif SCM_SIZEOF_INTPTR_T == 8 | |
551 | #define scm_to_intptr_t scm_to_int64 | |
552 | #define scm_from_intptr_t scm_from_int64 | |
553 | #else | |
554 | #error sizeof(intptr_t) is not 4 or 8. | |
555 | #endif | |
556 | ||
557 | #if SCM_SIZEOF_UINTPTR_T == 0 | |
558 | /* No uintptr_t; use size_t functions. */ | |
559 | #define scm_to_uintptr_t scm_to_size_t | |
560 | #define scm_from_uintptr_t scm_from_size_t | |
561 | #elif SCM_SIZEOF_UINTPTR_T == 4 | |
562 | #define scm_to_uintptr_t scm_to_uint32 | |
563 | #define scm_from_uintptr_t scm_from_uint32 | |
564 | #elif SCM_SIZEOF_UINTPTR_T == 8 | |
565 | #define scm_to_uintptr_t scm_to_uint64 | |
566 | #define scm_from_uintptr_t scm_from_uint64 | |
567 | #else | |
568 | #error sizeof(uintptr_t) is not 4 or 8. | |
569 | #endif | |
570 | ||
8507ec80 | 571 | /* conversion functions for double */ |
22008a78 | 572 | |
73e4de09 | 573 | SCM_API int scm_is_real (SCM val); |
55f26379 | 574 | SCM_API int scm_is_rational (SCM val); |
73e4de09 MV |
575 | SCM_API double scm_to_double (SCM val); |
576 | SCM_API SCM scm_from_double (double val); | |
577 | ||
8507ec80 MV |
578 | /* conversion functions for complex */ |
579 | ||
580 | SCM_API int scm_is_complex (SCM val); | |
581 | SCM_API SCM scm_c_make_rectangular (double re, double im); | |
582 | SCM_API SCM scm_c_make_polar (double mag, double ang); | |
583 | SCM_API double scm_c_real_part (SCM z); | |
584 | SCM_API double scm_c_imag_part (SCM z); | |
585 | SCM_API double scm_c_magnitude (SCM z); | |
586 | SCM_API double scm_c_angle (SCM z); | |
587 | ||
588 | SCM_API int scm_is_number (SCM val); | |
589 | ||
b57bf272 AW |
590 | /* If nonzero, tell gmp to use GC_malloc for its allocations. */ |
591 | SCM_API int scm_install_gmp_memory_functions; | |
592 | ||
102dbb6f | 593 | SCM_INTERNAL void scm_init_numbers (void); |
0f2d19dd | 594 | |
3c9a524f | 595 | #endif /* SCM_NUMBERS_H */ |
89e00824 ML |
596 | |
597 | /* | |
598 | Local Variables: | |
599 | c-file-style: "gnu" | |
600 | End: | |
601 | */ |