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,
   7  *   2008, 2009, 2010, 2011, 2013, 2014 Free Software Foundation, Inc.
   8  *
   9  * This library is free software; you can redistribute it and/or
  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.
  13  *
  14  * This library is distributed in the hope that it will be useful, but
  15  * WITHOUT ANY WARRANTY; without even the implied warranty of
  16  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  17  * Lesser General Public License for more details.
  18  *
  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
  21  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
  22  * 02110-1301 USA
  23  */
  24
  25 \f
  26
  27 #include <gmp.h>
  28
  29 #include "libguile/__scm.h"
  30 #include "libguile/print.h"
  31
  32 #ifndef SCM_T_WCHAR_DEFINED
  33 typedef scm_t_int32 scm_t_wchar;
  34 #define SCM_T_WCHAR_DEFINED
  35 #endif /* SCM_T_WCHAR_DEFINED */
  36
  37 \f
  38
  39 /* Immediate Numbers, also known as fixnums
  40  *
  41  * Inums are exact integer data that fits within an SCM word.  */
  42
  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
  48 #define SCM_I_FIXNUM_BIT         (SCM_LONG_BIT - 2)
  49 #define SCM_MOST_POSITIVE_FIXNUM ((SCM_T_SIGNED_BITS_MAX-3)/4)
  50 #define SCM_MOST_NEGATIVE_FIXNUM (-SCM_MOST_POSITIVE_FIXNUM-1)
  51
  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.
  56
  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))
  60 #else
  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
  67
  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
  84
  85 #define SCM_I_INUMP(x)  (2 & SCM_UNPACK (x))
  86 #define SCM_I_NINUMP(x) (!SCM_I_INUMP (x))
  87 #define SCM_I_MAKINUM(x) \
  88   (SCM_PACK ((((scm_t_bits) (x)) << 2) + scm_tc2_int))
  89
  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)
  93 #define SCM_FIXABLE(n) (SCM_POSFIXABLE (n) && SCM_NEGFIXABLE (n))
  94
  95
  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
  99
 100 /* SCM_MAXEXP is the maximum double precision exponent
 101  * SCM_FLTMAX is less than or scm_equal the largest single precision float
 102  */
 103
 104 #if SCM_HAVE_STDC_HEADERS
 105 # ifndef GO32
 106 #  include <float.h>
 107 #  ifdef __MINGW32__
 108 #   define copysign _copysign
 109 #   define finite _finite
 110 #  endif /* __MINGW32__ */
 111 # endif /* ndef GO32 */
 112 #endif /* def STDC_HEADERS */
 113
 114 #ifdef DBL_MAX_10_EXP
 115 # define SCM_MAXEXP DBL_MAX_10_EXP
 116 #else
 117 # define SCM_MAXEXP 308         /* IEEE doubles */
 118 #endif /* def DBL_MAX_10_EXP */
 119
 120 #ifdef FLT_MAX
 121 # define SCM_FLTMAX FLT_MAX
 122 #else
 123 # define SCM_FLTMAX 1e+23
 124 #endif /* def FLT_MAX */
 125
 126
 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.
 131  */
 132 #define SCM_INTBUFLEN (5 + SCM_CHAR_BIT*sizeof(scm_t_intmax))
 133
 134 \f
 135
 136 /* Numbers
 137  */
 138
 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
 144 /* Number subtype 1 to 3 (note the dependency on the predicates SCM_INEXACTP
 145  * and SCM_NUMP)  */
 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)
 149 #define scm_tc16_fraction       (scm_tc7_number + 4 * 256L)
 150
 151 #define SCM_INEXACTP(x) \
 152   (!SCM_IMP (x) && (0xfeff & SCM_CELL_TYPE (x)) == scm_tc16_real)
 153 #define SCM_REALP(x) (SCM_HAS_TYP16 (x, scm_tc16_real))
 154 #define SCM_COMPLEXP(x) (SCM_HAS_TYP16 (x, scm_tc16_complex))
 155
 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)
 159
 160 /* Each bignum is just an mpz_t stored in a double cell starting at word 1. */
 161 #define SCM_I_BIG_MPZ(x) (*((mpz_t *) (SCM_CELL_OBJECT_LOC((x),1))))
 162 #define SCM_BIGP(x) (SCM_HAS_TYP16 (x, scm_tc16_big))
 163
 164 #define SCM_NUMBERP(x) (SCM_I_INUMP(x) || SCM_NUMP(x))
 165 #define SCM_NUMP(x) (SCM_HAS_TYP7 (x, scm_tc7_number))
 166
 167 #define SCM_FRACTIONP(x) (SCM_HAS_TYP16 (x, scm_tc16_fraction))
 168 #define SCM_FRACTION_NUMERATOR(x) (SCM_CELL_OBJECT_1 (x))
 169 #define SCM_FRACTION_DENOMINATOR(x) (SCM_CELL_OBJECT_2 (x))
 170
 171 \f
 172
 173 typedef struct scm_t_double
 174 {
 175   SCM type;
 176   SCM pad;
 177   double real;
 178 } scm_t_double;
 179
 180 typedef struct scm_t_complex
 181 {
 182   SCM type;
 183   SCM pad;
 184   double real;
 185   double imag;
 186 } scm_t_complex;
 187
 188
 189 \f
 190
 191 SCM_API SCM scm_exact_p (SCM x);
 192 SCM_API int scm_is_exact (SCM x);
 193 SCM_API SCM scm_odd_p (SCM n);
 194 SCM_API SCM scm_even_p (SCM n);
 195 SCM_API SCM scm_finite_p (SCM x);
 196 SCM_API SCM scm_inf_p (SCM x);
 197 SCM_API SCM scm_nan_p (SCM x);
 198 SCM_API SCM scm_inf (void);
 199 SCM_API SCM scm_nan (void);
 200 SCM_API SCM scm_abs (SCM x);
 201 SCM_API SCM scm_quotient (SCM x, SCM y);
 202 SCM_API SCM scm_remainder (SCM x, SCM y);
 203 SCM_API SCM scm_modulo (SCM x, SCM y);
 204 SCM_API void scm_euclidean_divide (SCM x, SCM y, SCM *q, SCM *r);
 205 SCM_API SCM scm_euclidean_quotient (SCM x, SCM y);
 206 SCM_API SCM scm_euclidean_remainder (SCM x, SCM y);
 207 SCM_API void scm_floor_divide (SCM x, SCM y, SCM *q, SCM *r);
 208 SCM_API SCM scm_floor_quotient (SCM x, SCM y);
 209 SCM_API SCM scm_floor_remainder (SCM x, SCM y);
 210 SCM_API void scm_ceiling_divide (SCM x, SCM y, SCM *q, SCM *r);
 211 SCM_API SCM scm_ceiling_quotient (SCM x, SCM y);
 212 SCM_API SCM scm_ceiling_remainder (SCM x, SCM y);
 213 SCM_API void scm_truncate_divide (SCM x, SCM y, SCM *q, SCM *r);
 214 SCM_API SCM scm_truncate_quotient (SCM x, SCM y);
 215 SCM_API SCM scm_truncate_remainder (SCM x, SCM y);
 216 SCM_API void scm_centered_divide (SCM x, SCM y, SCM *q, SCM *r);
 217 SCM_API SCM scm_centered_quotient (SCM x, SCM y);
 218 SCM_API SCM scm_centered_remainder (SCM x, SCM y);
 219 SCM_API void scm_round_divide (SCM x, SCM y, SCM *q, SCM *r);
 220 SCM_API SCM scm_round_quotient (SCM x, SCM y);
 221 SCM_API SCM scm_round_remainder (SCM x, SCM y);
 222 SCM_API SCM scm_gcd (SCM x, SCM y);
 223 SCM_API SCM scm_lcm (SCM n1, SCM n2);
 224 SCM_API SCM scm_logand (SCM n1, SCM n2);
 225 SCM_API SCM scm_logior (SCM n1, SCM n2);
 226 SCM_API SCM scm_logxor (SCM n1, SCM n2);
 227 SCM_API SCM scm_logtest (SCM n1, SCM n2);
 228 SCM_API SCM scm_logbit_p (SCM n1, SCM n2);
 229 SCM_API SCM scm_lognot (SCM n);
 230 SCM_API SCM scm_modulo_expt (SCM n, SCM k, SCM m);
 231 SCM_API SCM scm_integer_expt (SCM z1, SCM z2);
 232 SCM_API SCM scm_ash (SCM n, SCM count);
 233 SCM_API SCM scm_round_ash (SCM n, SCM count);
 234 SCM_API SCM scm_bit_extract (SCM n, SCM start, SCM end);
 235 SCM_API SCM scm_logcount (SCM n);
 236 SCM_API SCM scm_integer_length (SCM n);
 237
 238 SCM_INTERNAL SCM scm_i_euclidean_divide (SCM x, SCM y);
 239 SCM_INTERNAL SCM scm_i_floor_divide (SCM x, SCM y);
 240 SCM_INTERNAL SCM scm_i_ceiling_divide (SCM x, SCM y);
 241 SCM_INTERNAL SCM scm_i_truncate_divide (SCM x, SCM y);
 242 SCM_INTERNAL SCM scm_i_centered_divide (SCM x, SCM y);
 243 SCM_INTERNAL SCM scm_i_round_divide (SCM x, SCM y);
 244
 245 SCM_INTERNAL SCM scm_i_gcd (SCM x, SCM y, SCM rest);
 246 SCM_INTERNAL SCM scm_i_lcm (SCM x, SCM y, SCM rest);
 247 SCM_INTERNAL SCM scm_i_logand (SCM x, SCM y, SCM rest);
 248 SCM_INTERNAL SCM scm_i_logior (SCM x, SCM y, SCM rest);
 249 SCM_INTERNAL SCM scm_i_logxor (SCM x, SCM y, SCM rest);
 250
 251 SCM_API size_t scm_iint2str (scm_t_intmax num, int rad, char *p);
 252 SCM_API size_t scm_iuint2str (scm_t_uintmax num, int rad, char *p);
 253 SCM_API SCM scm_number_to_string (SCM x, SCM radix);
 254 SCM_API int scm_print_real (SCM sexp, SCM port, scm_print_state *pstate);
 255 SCM_API int scm_print_complex (SCM sexp, SCM port, scm_print_state *pstate);
 256 SCM_API int scm_bigprint (SCM exp, SCM port, scm_print_state *pstate);
 257 SCM_API SCM scm_c_locale_stringn_to_number (const char *mem, size_t len,
 258                                             unsigned int radix);
 259 SCM_INTERNAL SCM scm_i_string_to_number (SCM str, unsigned int radix);
 260 SCM_API SCM scm_string_to_number (SCM str, SCM radix);
 261 SCM_API SCM scm_bigequal (SCM x, SCM y);
 262 SCM_API SCM scm_real_equalp (SCM x, SCM y);
 263 SCM_API SCM scm_complex_equalp (SCM x, SCM y);
 264 SCM_API SCM scm_number_p (SCM x);
 265 SCM_API SCM scm_complex_p (SCM x);
 266 SCM_API SCM scm_real_p (SCM x);
 267 SCM_API SCM scm_rational_p (SCM z);
 268 SCM_API SCM scm_integer_p (SCM x);
 269 SCM_API SCM scm_exact_integer_p (SCM x);
 270 SCM_API SCM scm_inexact_p (SCM x);
 271 SCM_API int scm_is_inexact (SCM x);
 272 SCM_API SCM scm_num_eq_p (SCM x, SCM y);
 273 SCM_API SCM scm_less_p (SCM x, SCM y);
 274 SCM_API SCM scm_gr_p (SCM x, SCM y);
 275 SCM_API SCM scm_leq_p (SCM x, SCM y);
 276 SCM_API SCM scm_geq_p (SCM x, SCM y);
 277 SCM_API SCM scm_zero_p (SCM z);
 278 SCM_API SCM scm_positive_p (SCM x);
 279 SCM_API SCM scm_negative_p (SCM x);
 280 SCM_API SCM scm_max (SCM x, SCM y);
 281 SCM_API SCM scm_min (SCM x, SCM y);
 282 SCM_API SCM scm_sum (SCM x, SCM y);
 283 SCM_API SCM scm_oneplus (SCM x);
 284 SCM_API SCM scm_difference (SCM x, SCM y);
 285 SCM_API SCM scm_oneminus (SCM x);
 286 SCM_API SCM scm_product (SCM x, SCM y);
 287 SCM_API SCM scm_divide (SCM x, SCM y);
 288 SCM_API SCM scm_floor (SCM x);
 289 SCM_API SCM scm_ceiling (SCM x);
 290 SCM_API double scm_c_truncate (double x);
 291 SCM_API double scm_c_round (double x);
 292 SCM_API SCM scm_truncate_number (SCM x);
 293 SCM_API SCM scm_round_number (SCM x);
 294 SCM_API SCM scm_expt (SCM z1, SCM z2);
 295 SCM_API SCM scm_sin (SCM z);
 296 SCM_API SCM scm_cos (SCM z);
 297 SCM_API SCM scm_tan (SCM z);
 298 SCM_API SCM scm_sinh (SCM z);
 299 SCM_API SCM scm_cosh (SCM z);
 300 SCM_API SCM scm_tanh (SCM z);
 301 SCM_API SCM scm_asin (SCM z);
 302 SCM_API SCM scm_acos (SCM z);
 303 SCM_API SCM scm_atan (SCM x, SCM y);
 304 SCM_API SCM scm_sys_asinh (SCM z);
 305 SCM_API SCM scm_sys_acosh (SCM z);
 306 SCM_API SCM scm_sys_atanh (SCM z);
 307 SCM_API SCM scm_make_rectangular (SCM z1, SCM z2);
 308 SCM_API SCM scm_make_polar (SCM z1, SCM z2);
 309 SCM_API SCM scm_real_part (SCM z);
 310 SCM_API SCM scm_imag_part (SCM z);
 311 SCM_API SCM scm_magnitude (SCM z);
 312 SCM_API SCM scm_angle (SCM z);
 313 SCM_API SCM scm_exact_to_inexact (SCM z);
 314 SCM_API SCM scm_inexact_to_exact (SCM z);
 315 SCM_API SCM scm_trunc (SCM x);
 316 SCM_API SCM scm_log (SCM z);
 317 SCM_API SCM scm_log10 (SCM z);
 318 SCM_API SCM scm_exp (SCM z);
 319 SCM_API SCM scm_sqrt (SCM z);
 320 SCM_API void scm_exact_integer_sqrt (SCM k, SCM *s, SCM *r);
 321
 322 SCM_INTERNAL SCM scm_i_min (SCM x, SCM y, SCM rest);
 323 SCM_INTERNAL SCM scm_i_max (SCM x, SCM y, SCM rest);
 324 SCM_INTERNAL SCM scm_i_sum (SCM x, SCM y, SCM rest);
 325 SCM_INTERNAL SCM scm_i_difference (SCM x, SCM y, SCM rest);
 326 SCM_INTERNAL SCM scm_i_product (SCM x, SCM y, SCM rest);
 327 SCM_INTERNAL SCM scm_i_divide (SCM x, SCM y, SCM rest);
 328 SCM_INTERNAL SCM scm_i_exact_integer_sqrt (SCM k);
 329
 330 /* bignum internal functions */
 331 SCM_INTERNAL SCM scm_i_mkbig (void);
 332 SCM_API /* FIXME: not internal */ SCM scm_i_normbig (SCM x);
 333 SCM_INTERNAL int scm_i_bigcmp (SCM a, SCM b);
 334 SCM_INTERNAL SCM scm_i_dbl2big (double d);
 335 SCM_INTERNAL SCM scm_i_dbl2num (double d);
 336 SCM_API /* FIXME: not internal */ double scm_i_big2dbl (SCM b);
 337 SCM_API /* FIXME: not internal */ SCM scm_i_long2big (long n);
 338 SCM_API /* FIXME: not internal */ SCM scm_i_ulong2big (unsigned long n);
 339 SCM_API /* FIXME: not internal */ SCM scm_i_clonebig (SCM src_big, int same_sign_p);
 340
 341 /* ratio functions */
 342 SCM_API SCM scm_rationalize (SCM x, SCM err);
 343 SCM_API SCM scm_numerator (SCM z);
 344 SCM_API SCM scm_denominator (SCM z);
 345
 346 /* fraction internal functions */
 347 SCM_INTERNAL double scm_i_fraction2double (SCM z);
 348 SCM_INTERNAL SCM scm_i_fraction_equalp (SCM x, SCM y);
 349 SCM_INTERNAL int scm_i_print_fraction (SCM sexp, SCM port, scm_print_state *pstate);
 350
 351 /* general internal functions */
 352 SCM_INTERNAL void scm_i_print_double (double val, SCM port);
 353 SCM_INTERNAL void scm_i_print_complex (double real, double imag, SCM port);
 354
 355 /* conversion functions for integers */
 356
 357 SCM_API int scm_is_integer (SCM val);
 358 SCM_API int scm_is_exact_integer (SCM val);
 359 SCM_API int scm_is_signed_integer (SCM val,
 360                                    scm_t_intmax min, scm_t_intmax max);
 361 SCM_API int scm_is_unsigned_integer (SCM val,
 362                                      scm_t_uintmax min, scm_t_uintmax max);
 363
 364 SCM_API SCM scm_from_signed_integer (scm_t_intmax val);
 365 SCM_API SCM scm_from_unsigned_integer (scm_t_uintmax val);
 366
 367 SCM_API scm_t_intmax scm_to_signed_integer (SCM val,
 368                                             scm_t_intmax min,
 369                                             scm_t_intmax max);
 370 SCM_API scm_t_uintmax scm_to_unsigned_integer (SCM val,
 371                                                scm_t_uintmax min,
 372                                                scm_t_uintmax max);
 373
 374 SCM_API scm_t_int8   scm_to_int8     (SCM x);
 375 SCM_API SCM          scm_from_int8   (scm_t_int8 x);
 376
 377 SCM_API scm_t_uint8  scm_to_uint8    (SCM x);
 378 SCM_API SCM          scm_from_uint8  (scm_t_uint8 x);
 379
 380 SCM_API scm_t_int16  scm_to_int16    (SCM x);
 381 SCM_API SCM          scm_from_int16  (scm_t_int16 x);
 382
 383 SCM_API scm_t_uint16 scm_to_uint16   (SCM x);
 384 SCM_API SCM          scm_from_uint16 (scm_t_uint16 x);
 385
 386 SCM_API scm_t_int32  scm_to_int32    (SCM x);
 387 SCM_API SCM          scm_from_int32  (scm_t_int32 x);
 388
 389 SCM_API scm_t_uint32 scm_to_uint32   (SCM x);
 390 SCM_API SCM          scm_from_uint32 (scm_t_uint32 x);
 391
 392 SCM_API scm_t_wchar  scm_to_wchar    (SCM x);
 393 SCM_API SCM          scm_from_wchar  (scm_t_wchar x);
 394
 395 SCM_API scm_t_int64  scm_to_int64    (SCM x);
 396 SCM_API SCM          scm_from_int64  (scm_t_int64 x);
 397
 398 SCM_API scm_t_uint64 scm_to_uint64   (SCM x);
 399 SCM_API SCM          scm_from_uint64 (scm_t_uint64 x);
 400
 401 SCM_API void scm_to_mpz (SCM x, mpz_t rop);
 402 SCM_API SCM  scm_from_mpz (mpz_t rop);
 403
 404
 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
 421 #if CHAR_MIN == 0
 422 #define scm_to_char    scm_to_uint8
 423 #define scm_from_char  scm_from_uint8
 424 #else
 425 #define scm_to_char    scm_to_int8
 426 #define scm_from_char  scm_from_int8
 427 #endif
 428 #else
 429 #error sizeof(char) is not 1.
 430 #endif
 431
 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
 449 #else
 450 #error sizeof(short) is not 1, 2, or 4.
 451 #endif
 452 #endif
 453 #endif
 454
 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
 470
 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
 486
 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
 499 #error sizeof(scm_t_intmax) is not 4 or 8.
 500 #endif
 501 #endif
 502
 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
 514
 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
 529 #endif
 530
 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
 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
 571 /* conversion functions for double */
 572
 573 SCM_API int scm_is_real (SCM val);
 574 SCM_API int scm_is_rational (SCM val);
 575 SCM_API double scm_to_double (SCM val);
 576 SCM_API SCM scm_from_double (double val);
 577
 578 /* conversion functions for complex */
 579
 580 SCM_API int scm_is_complex (SCM val);
 581 SCM_API SCM scm_c_make_rectangular (double re, double im);
 582 SCM_API SCM scm_c_make_polar (double mag, double ang);
 583 SCM_API double scm_c_real_part (SCM z);
 584 SCM_API double scm_c_imag_part (SCM z);
 585 SCM_API double scm_c_magnitude (SCM z);
 586 SCM_API double scm_c_angle (SCM z);
 587
 588 SCM_API int scm_is_number (SCM val);
 589
 590 /* If nonzero, tell gmp to use GC_malloc for its allocations.  */
 591 SCM_API int scm_install_gmp_memory_functions;
 592
 593 SCM_INTERNAL void scm_init_numbers (void);
 594
 595 #endif  /* SCM_NUMBERS_H */
 596
 597 /*
 598   Local Variables:
 599   c-file-style: "gnu"
 600   End:
 601 */