Fix infinite loop in expander
[bpt/guile.git] / libguile / numbers.h
CommitLineData
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
33typedef 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 173typedef 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 180typedef 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 191SCM_API SCM scm_exact_p (SCM x);
022dda69 192SCM_API int scm_is_exact (SCM x);
33b001fd
MV
193SCM_API SCM scm_odd_p (SCM n);
194SCM_API SCM scm_even_p (SCM n);
2519490c
MW
195SCM_API SCM scm_finite_p (SCM x);
196SCM_API SCM scm_inf_p (SCM x);
197SCM_API SCM scm_nan_p (SCM x);
232898d9
MV
198SCM_API SCM scm_inf (void);
199SCM_API SCM scm_nan (void);
33b001fd
MV
200SCM_API SCM scm_abs (SCM x);
201SCM_API SCM scm_quotient (SCM x, SCM y);
202SCM_API SCM scm_remainder (SCM x, SCM y);
203SCM_API SCM scm_modulo (SCM x, SCM y);
5fbf680b 204SCM_API void scm_euclidean_divide (SCM x, SCM y, SCM *q, SCM *r);
ff62c168
MW
205SCM_API SCM scm_euclidean_quotient (SCM x, SCM y);
206SCM_API SCM scm_euclidean_remainder (SCM x, SCM y);
8f9da340
MW
207SCM_API void scm_floor_divide (SCM x, SCM y, SCM *q, SCM *r);
208SCM_API SCM scm_floor_quotient (SCM x, SCM y);
209SCM_API SCM scm_floor_remainder (SCM x, SCM y);
210SCM_API void scm_ceiling_divide (SCM x, SCM y, SCM *q, SCM *r);
211SCM_API SCM scm_ceiling_quotient (SCM x, SCM y);
212SCM_API SCM scm_ceiling_remainder (SCM x, SCM y);
213SCM_API void scm_truncate_divide (SCM x, SCM y, SCM *q, SCM *r);
214SCM_API SCM scm_truncate_quotient (SCM x, SCM y);
215SCM_API SCM scm_truncate_remainder (SCM x, SCM y);
5fbf680b 216SCM_API void scm_centered_divide (SCM x, SCM y, SCM *q, SCM *r);
ff62c168
MW
217SCM_API SCM scm_centered_quotient (SCM x, SCM y);
218SCM_API SCM scm_centered_remainder (SCM x, SCM y);
8f9da340
MW
219SCM_API void scm_round_divide (SCM x, SCM y, SCM *q, SCM *r);
220SCM_API SCM scm_round_quotient (SCM x, SCM y);
221SCM_API SCM scm_round_remainder (SCM x, SCM y);
33b001fd
MV
222SCM_API SCM scm_gcd (SCM x, SCM y);
223SCM_API SCM scm_lcm (SCM n1, SCM n2);
224SCM_API SCM scm_logand (SCM n1, SCM n2);
225SCM_API SCM scm_logior (SCM n1, SCM n2);
226SCM_API SCM scm_logxor (SCM n1, SCM n2);
227SCM_API SCM scm_logtest (SCM n1, SCM n2);
228SCM_API SCM scm_logbit_p (SCM n1, SCM n2);
229SCM_API SCM scm_lognot (SCM n);
d885e204 230SCM_API SCM scm_modulo_expt (SCM n, SCM k, SCM m);
33b001fd 231SCM_API SCM scm_integer_expt (SCM z1, SCM z2);
e08a12b5
MW
232SCM_API SCM scm_ash (SCM n, SCM count);
233SCM_API SCM scm_round_ash (SCM n, SCM count);
33b001fd
MV
234SCM_API SCM scm_bit_extract (SCM n, SCM start, SCM end);
235SCM_API SCM scm_logcount (SCM n);
236SCM_API SCM scm_integer_length (SCM n);
1be6b49c 237
5fbf680b 238SCM_INTERNAL SCM scm_i_euclidean_divide (SCM x, SCM y);
8f9da340
MW
239SCM_INTERNAL SCM scm_i_floor_divide (SCM x, SCM y);
240SCM_INTERNAL SCM scm_i_ceiling_divide (SCM x, SCM y);
241SCM_INTERNAL SCM scm_i_truncate_divide (SCM x, SCM y);
5fbf680b 242SCM_INTERNAL SCM scm_i_centered_divide (SCM x, SCM y);
8f9da340 243SCM_INTERNAL SCM scm_i_round_divide (SCM x, SCM y);
5fbf680b 244
78d3deb1
AW
245SCM_INTERNAL SCM scm_i_gcd (SCM x, SCM y, SCM rest);
246SCM_INTERNAL SCM scm_i_lcm (SCM x, SCM y, SCM rest);
247SCM_INTERNAL SCM scm_i_logand (SCM x, SCM y, SCM rest);
248SCM_INTERNAL SCM scm_i_logior (SCM x, SCM y, SCM rest);
249SCM_INTERNAL SCM scm_i_logxor (SCM x, SCM y, SCM rest);
250
2881e77b
MV
251SCM_API size_t scm_iint2str (scm_t_intmax num, int rad, char *p);
252SCM_API size_t scm_iuint2str (scm_t_uintmax num, int rad, char *p);
33b001fd
MV
253SCM_API SCM scm_number_to_string (SCM x, SCM radix);
254SCM_API int scm_print_real (SCM sexp, SCM port, scm_print_state *pstate);
255SCM_API int scm_print_complex (SCM sexp, SCM port, scm_print_state *pstate);
256SCM_API int scm_bigprint (SCM exp, SCM port, scm_print_state *pstate);
2b829bbb
KR
257SCM_API SCM scm_c_locale_stringn_to_number (const char *mem, size_t len,
258 unsigned int radix);
3f47e526 259SCM_INTERNAL SCM scm_i_string_to_number (SCM str, unsigned int radix);
33b001fd 260SCM_API SCM scm_string_to_number (SCM str, SCM radix);
33b001fd
MV
261SCM_API SCM scm_bigequal (SCM x, SCM y);
262SCM_API SCM scm_real_equalp (SCM x, SCM y);
263SCM_API SCM scm_complex_equalp (SCM x, SCM y);
264SCM_API SCM scm_number_p (SCM x);
8507ec80 265SCM_API SCM scm_complex_p (SCM x);
33b001fd 266SCM_API SCM scm_real_p (SCM x);
8507ec80 267SCM_API SCM scm_rational_p (SCM z);
33b001fd 268SCM_API SCM scm_integer_p (SCM x);
900a897c 269SCM_API SCM scm_exact_integer_p (SCM x);
33b001fd 270SCM_API SCM scm_inexact_p (SCM x);
022dda69 271SCM_API int scm_is_inexact (SCM x);
33b001fd
MV
272SCM_API SCM scm_num_eq_p (SCM x, SCM y);
273SCM_API SCM scm_less_p (SCM x, SCM y);
274SCM_API SCM scm_gr_p (SCM x, SCM y);
275SCM_API SCM scm_leq_p (SCM x, SCM y);
276SCM_API SCM scm_geq_p (SCM x, SCM y);
277SCM_API SCM scm_zero_p (SCM z);
278SCM_API SCM scm_positive_p (SCM x);
279SCM_API SCM scm_negative_p (SCM x);
280SCM_API SCM scm_max (SCM x, SCM y);
281SCM_API SCM scm_min (SCM x, SCM y);
282SCM_API SCM scm_sum (SCM x, SCM y);
40882e3d 283SCM_API SCM scm_oneplus (SCM x);
33b001fd 284SCM_API SCM scm_difference (SCM x, SCM y);
40882e3d 285SCM_API SCM scm_oneminus (SCM x);
33b001fd 286SCM_API SCM scm_product (SCM x, SCM y);
33b001fd 287SCM_API SCM scm_divide (SCM x, SCM y);
f92e85f7
MV
288SCM_API SCM scm_floor (SCM x);
289SCM_API SCM scm_ceiling (SCM x);
3101f40f
MV
290SCM_API double scm_c_truncate (double x);
291SCM_API double scm_c_round (double x);
f92e85f7
MV
292SCM_API SCM scm_truncate_number (SCM x);
293SCM_API SCM scm_round_number (SCM x);
6fc4d012 294SCM_API SCM scm_expt (SCM z1, SCM z2);
ad79736c
AW
295SCM_API SCM scm_sin (SCM z);
296SCM_API SCM scm_cos (SCM z);
297SCM_API SCM scm_tan (SCM z);
298SCM_API SCM scm_sinh (SCM z);
299SCM_API SCM scm_cosh (SCM z);
300SCM_API SCM scm_tanh (SCM z);
301SCM_API SCM scm_asin (SCM z);
302SCM_API SCM scm_acos (SCM z);
303SCM_API SCM scm_atan (SCM x, SCM y);
304SCM_API SCM scm_sys_asinh (SCM z);
305SCM_API SCM scm_sys_acosh (SCM z);
306SCM_API SCM scm_sys_atanh (SCM z);
33b001fd
MV
307SCM_API SCM scm_make_rectangular (SCM z1, SCM z2);
308SCM_API SCM scm_make_polar (SCM z1, SCM z2);
309SCM_API SCM scm_real_part (SCM z);
310SCM_API SCM scm_imag_part (SCM z);
311SCM_API SCM scm_magnitude (SCM z);
312SCM_API SCM scm_angle (SCM z);
313SCM_API SCM scm_exact_to_inexact (SCM z);
314SCM_API SCM scm_inexact_to_exact (SCM z);
315SCM_API SCM scm_trunc (SCM x);
8ab3d8a0
KR
316SCM_API SCM scm_log (SCM z);
317SCM_API SCM scm_log10 (SCM z);
318SCM_API SCM scm_exp (SCM z);
319SCM_API SCM scm_sqrt (SCM z);
882c8963 320SCM_API void scm_exact_integer_sqrt (SCM k, SCM *s, SCM *r);
33b001fd 321
78d3deb1
AW
322SCM_INTERNAL SCM scm_i_min (SCM x, SCM y, SCM rest);
323SCM_INTERNAL SCM scm_i_max (SCM x, SCM y, SCM rest);
324SCM_INTERNAL SCM scm_i_sum (SCM x, SCM y, SCM rest);
325SCM_INTERNAL SCM scm_i_difference (SCM x, SCM y, SCM rest);
326SCM_INTERNAL SCM scm_i_product (SCM x, SCM y, SCM rest);
327SCM_INTERNAL SCM scm_i_divide (SCM x, SCM y, SCM rest);
882c8963 328SCM_INTERNAL SCM scm_i_exact_integer_sqrt (SCM k);
78d3deb1 329
d3e7e88d 330/* bignum internal functions */
102dbb6f
LC
331SCM_INTERNAL SCM scm_i_mkbig (void);
332SCM_API /* FIXME: not internal */ SCM scm_i_normbig (SCM x);
333SCM_INTERNAL int scm_i_bigcmp (SCM a, SCM b);
334SCM_INTERNAL SCM scm_i_dbl2big (double d);
335SCM_INTERNAL SCM scm_i_dbl2num (double d);
336SCM_API /* FIXME: not internal */ double scm_i_big2dbl (SCM b);
337SCM_API /* FIXME: not internal */ SCM scm_i_long2big (long n);
338SCM_API /* FIXME: not internal */ SCM scm_i_ulong2big (unsigned long n);
339SCM_API /* FIXME: not internal */ SCM scm_i_clonebig (SCM src_big, int same_sign_p);
d3e7e88d 340
f92e85f7 341/* ratio functions */
f92e85f7
MV
342SCM_API SCM scm_rationalize (SCM x, SCM err);
343SCM_API SCM scm_numerator (SCM z);
344SCM_API SCM scm_denominator (SCM z);
f92e85f7
MV
345
346/* fraction internal functions */
102dbb6f
LC
347SCM_INTERNAL double scm_i_fraction2double (SCM z);
348SCM_INTERNAL SCM scm_i_fraction_equalp (SCM x, SCM y);
349SCM_INTERNAL int scm_i_print_fraction (SCM sexp, SCM port, scm_print_state *pstate);
f92e85f7 350
b479fe9a 351/* general internal functions */
102dbb6f
LC
352SCM_INTERNAL void scm_i_print_double (double val, SCM port);
353SCM_INTERNAL void scm_i_print_complex (double real, double imag, SCM port);
b479fe9a 354
d1b5b400 355/* conversion functions for integers */
73e4de09
MV
356
357SCM_API int scm_is_integer (SCM val);
900a897c 358SCM_API int scm_is_exact_integer (SCM val);
73e4de09
MV
359SCM_API int scm_is_signed_integer (SCM val,
360 scm_t_intmax min, scm_t_intmax max);
361SCM_API int scm_is_unsigned_integer (SCM val,
362 scm_t_uintmax min, scm_t_uintmax max);
363
364SCM_API SCM scm_from_signed_integer (scm_t_intmax val);
365SCM_API SCM scm_from_unsigned_integer (scm_t_uintmax val);
366
367SCM_API scm_t_intmax scm_to_signed_integer (SCM val,
368 scm_t_intmax min,
369 scm_t_intmax max);
370SCM_API scm_t_uintmax scm_to_unsigned_integer (SCM val,
371 scm_t_uintmax min,
372 scm_t_uintmax max);
373
bfd7932e
MV
374SCM_API scm_t_int8 scm_to_int8 (SCM x);
375SCM_API SCM scm_from_int8 (scm_t_int8 x);
376
377SCM_API scm_t_uint8 scm_to_uint8 (SCM x);
378SCM_API SCM scm_from_uint8 (scm_t_uint8 x);
379
380SCM_API scm_t_int16 scm_to_int16 (SCM x);
381SCM_API SCM scm_from_int16 (scm_t_int16 x);
382
383SCM_API scm_t_uint16 scm_to_uint16 (SCM x);
384SCM_API SCM scm_from_uint16 (scm_t_uint16 x);
385
386SCM_API scm_t_int32 scm_to_int32 (SCM x);
387SCM_API SCM scm_from_int32 (scm_t_int32 x);
388
389SCM_API scm_t_uint32 scm_to_uint32 (SCM x);
390SCM_API SCM scm_from_uint32 (scm_t_uint32 x);
391
904a78f1
MG
392SCM_API scm_t_wchar scm_to_wchar (SCM x);
393SCM_API SCM scm_from_wchar (scm_t_wchar x);
394
bfd7932e
MV
395SCM_API scm_t_int64 scm_to_int64 (SCM x);
396SCM_API SCM scm_from_int64 (scm_t_int64 x);
397
398SCM_API scm_t_uint64 scm_to_uint64 (SCM x);
399SCM_API SCM scm_from_uint64 (scm_t_uint64 x);
400
cd036260
MV
401SCM_API void scm_to_mpz (SCM x, mpz_t rop);
402SCM_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 573SCM_API int scm_is_real (SCM val);
55f26379 574SCM_API int scm_is_rational (SCM val);
73e4de09
MV
575SCM_API double scm_to_double (SCM val);
576SCM_API SCM scm_from_double (double val);
577
8507ec80
MV
578/* conversion functions for complex */
579
580SCM_API int scm_is_complex (SCM val);
581SCM_API SCM scm_c_make_rectangular (double re, double im);
582SCM_API SCM scm_c_make_polar (double mag, double ang);
583SCM_API double scm_c_real_part (SCM z);
584SCM_API double scm_c_imag_part (SCM z);
585SCM_API double scm_c_magnitude (SCM z);
586SCM_API double scm_c_angle (SCM z);
587
588SCM_API int scm_is_number (SCM val);
589
b57bf272
AW
590/* If nonzero, tell gmp to use GC_malloc for its allocations. */
591SCM_API int scm_install_gmp_memory_functions;
592
102dbb6f 593SCM_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*/