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