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