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