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