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