[bpt/guile.git] / libguile / numbers.h

/* classes: h_files */

#ifndef SCM_NUMBERS_H
#define SCM_NUMBERS_H

/* Copyright (C) 1995,1996,1998,2000,2001,2002,2003,2004,2005, 2006, 2008, 2009, 2010 Free Software Foundation, Inc.
 *
 * This library is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public License
 * as published by the Free Software Foundation; either version 3 of
 * the License, or (at your option) any later version.
 *
 * This library is distributed in the hope that it will be useful, but
 * WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with this library; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
 * 02110-1301 USA
 */

\f

#include <gmp.h>

#include "libguile/__scm.h"
#include "libguile/print.h"

#ifndef SCM_T_WCHAR_DEFINED
typedef scm_t_int32 scm_t_wchar;
#define SCM_T_WCHAR_DEFINED
#endif /* SCM_T_WCHAR_DEFINED */

\f

/* Immediate Numbers, also known as fixnums
 *
 * Inums are exact integer data that fits within an SCM word.  */

/* SCM_T_SIGNED_MAX is                 (- (expt 2 n) 1),
 * SCM_MOST_POSITIVE_FIXNUM should be  (- (expt 2 (- n 2)) 1)
 * which is the same as                (/ (- (expt 2 n) 4) 4)
 */

#define SCM_I_FIXNUM_BIT         (SCM_LONG_BIT - 2)
#define SCM_MOST_POSITIVE_FIXNUM ((SCM_T_SIGNED_BITS_MAX-3)/4)
#define SCM_MOST_NEGATIVE_FIXNUM (-SCM_MOST_POSITIVE_FIXNUM-1)

/* SCM_SRS is signed right shift */
#if (-1 == (((-1) << 2) + 2) >> 2)
# define SCM_SRS(x, y) ((x) >> (y))
#else
# define SCM_SRS(x, y) ((x) < 0 ? ~((~(x)) >> (y)) : ((x) >> (y)))
#endif /* (-1 == (((-1) << 2) + 2) >> 2) */


#define SCM_I_INUMP(x)	(2 & SCM_UNPACK (x))
#define SCM_I_NINUMP(x) (!SCM_I_INUMP (x))
#define SCM_I_MAKINUM(x) \
  (SCM_PACK ((((scm_t_signed_bits) (x)) << 2) + scm_tc2_int))
#define SCM_I_INUM(x)   (SCM_SRS ((scm_t_signed_bits) SCM_UNPACK (x), 2))

/* SCM_FIXABLE is true if its long argument can be encoded in an SCM_INUM. */
#define SCM_POSFIXABLE(n) ((n) <= SCM_MOST_POSITIVE_FIXNUM)
#define SCM_NEGFIXABLE(n) ((n) >= SCM_MOST_NEGATIVE_FIXNUM)
#define SCM_FIXABLE(n) (SCM_POSFIXABLE (n) && SCM_NEGFIXABLE (n))


/* A name for 0. */
#define SCM_INUM0 (SCM_I_MAKINUM (0))

/* SCM_MAXEXP is the maximum double precision exponent
 * SCM_FLTMAX is less than or scm_equal the largest single precision float
 */

#if SCM_HAVE_STDC_HEADERS
# ifndef GO32
#  include <float.h>
#  ifdef __MINGW32__
#   define copysign _copysign
#   define finite _finite
#  endif /* __MINGW32__ */
# endif /* ndef GO32 */
#endif /* def STDC_HEADERS */

#ifdef DBL_MAX_10_EXP
# define SCM_MAXEXP DBL_MAX_10_EXP
#else
# define SCM_MAXEXP 308		/* IEEE doubles */
#endif /* def DBL_MAX_10_EXP */

#ifdef FLT_MAX
# define SCM_FLTMAX FLT_MAX
#else
# define SCM_FLTMAX 1e+23
#endif /* def FLT_MAX */


/* SCM_INTBUFLEN is the maximum number of characters neccessary for
 * the printed or scm_string representation of an scm_t_intmax in
 * radix 2.  The buffer passed to scm_iint2str and scm_iuint2str must
 * be of this size, for example.
 */
#define SCM_INTBUFLEN (5 + SCM_CHAR_BIT*sizeof(scm_t_intmax))

\f

/* Numbers 
 */


/* Note that scm_tc16_real and scm_tc16_complex are given tc16-codes that only
 * differ in one bit: This way, checking if an object is an inexact number can
 * be done quickly (using the TYP16S macro).  */

/* Number subtype 1 to 3 (note the dependency on the predicates SCM_INEXACTP
 * and SCM_NUMP)  */
#define scm_tc16_big		(scm_tc7_number + 1 * 256L)
#define scm_tc16_real           (scm_tc7_number + 2 * 256L)
#define scm_tc16_complex        (scm_tc7_number + 3 * 256L)
#define scm_tc16_fraction       (scm_tc7_number + 4 * 256L)

#define SCM_INEXACTP(x) \
  (!SCM_IMP (x) && (0xfeff & SCM_CELL_TYPE (x)) == scm_tc16_real)
#define SCM_REALP(x) (!SCM_IMP (x) && SCM_TYP16 (x) == scm_tc16_real)
#define SCM_COMPLEXP(x) (!SCM_IMP (x) && SCM_TYP16 (x) == scm_tc16_complex)

#define SCM_REAL_VALUE(x) (((scm_t_double *) SCM2PTR (x))->real)
#define SCM_COMPLEX_REAL(x) (((scm_t_complex *) SCM2PTR (x))->real)
#define SCM_COMPLEX_IMAG(x) (((scm_t_complex *) SCM2PTR (x))->imag)

/* Each bignum is just an mpz_t stored in a double cell starting at word 1. */
#define SCM_I_BIG_MPZ(x) (*((mpz_t *) (SCM_CELL_OBJECT_LOC((x),1))))
#define SCM_BIGP(x) (!SCM_IMP (x) && SCM_TYP16 (x) == scm_tc16_big)

#define SCM_NUMBERP(x) (SCM_I_INUMP(x) || SCM_NUMP(x))
#define SCM_NUMP(x) (!SCM_IMP(x) \
  && (((0xfcff & SCM_CELL_TYPE (x)) == scm_tc7_number) \
      || ((0xfbff & SCM_CELL_TYPE (x)) == scm_tc7_number)))
/* 0xfcff (#b1100) for 0 free, 1 big, 2 real, 3 complex, then 0xfbff (#b1011) for 4 fraction */

#define SCM_FRACTIONP(x) (!SCM_IMP (x) && SCM_TYP16 (x) == scm_tc16_fraction)
#define SCM_FRACTION_NUMERATOR(x) (SCM_CELL_OBJECT_1 (x))
#define SCM_FRACTION_DENOMINATOR(x) (SCM_CELL_OBJECT_2 (x))

\f

typedef struct scm_t_double
{
  SCM type;
  SCM pad;
  double real;
} scm_t_double;

typedef struct scm_t_complex
{
  SCM type;
  SCM pad;
  double real;
  double imag;
} scm_t_complex;


\f

SCM_API SCM scm_exact_p (SCM x);
SCM_API SCM scm_odd_p (SCM n);
SCM_API SCM scm_even_p (SCM n);
SCM_API SCM scm_inf_p (SCM n);
SCM_API SCM scm_nan_p (SCM n);
SCM_API SCM scm_inf (void);
SCM_API SCM scm_nan (void);
SCM_API SCM scm_abs (SCM x);
SCM_API SCM scm_quotient (SCM x, SCM y);
SCM_API SCM scm_remainder (SCM x, SCM y);
SCM_API SCM scm_modulo (SCM x, SCM y);
SCM_API SCM scm_gcd (SCM x, SCM y);
SCM_API SCM scm_lcm (SCM n1, SCM n2);
SCM_API SCM scm_logand (SCM n1, SCM n2);
SCM_API SCM scm_logior (SCM n1, SCM n2);
SCM_API SCM scm_logxor (SCM n1, SCM n2);
SCM_API SCM scm_logtest (SCM n1, SCM n2);
SCM_API SCM scm_logbit_p (SCM n1, SCM n2);
SCM_API SCM scm_lognot (SCM n);
SCM_API SCM scm_modulo_expt (SCM n, SCM k, SCM m);
SCM_API SCM scm_integer_expt (SCM z1, SCM z2);
SCM_API SCM scm_ash (SCM n, SCM cnt);
SCM_API SCM scm_bit_extract (SCM n, SCM start, SCM end);
SCM_API SCM scm_logcount (SCM n);
SCM_API SCM scm_integer_length (SCM n);

SCM_INTERNAL SCM scm_i_gcd (SCM x, SCM y, SCM rest);
SCM_INTERNAL SCM scm_i_lcm (SCM x, SCM y, SCM rest);
SCM_INTERNAL SCM scm_i_logand (SCM x, SCM y, SCM rest);
SCM_INTERNAL SCM scm_i_logior (SCM x, SCM y, SCM rest);
SCM_INTERNAL SCM scm_i_logxor (SCM x, SCM y, SCM rest);

SCM_API size_t scm_iint2str (scm_t_intmax num, int rad, char *p);
SCM_API size_t scm_iuint2str (scm_t_uintmax num, int rad, char *p);
SCM_API SCM scm_number_to_string (SCM x, SCM radix);
SCM_API int scm_print_real (SCM sexp, SCM port, scm_print_state *pstate);
SCM_API int scm_print_complex (SCM sexp, SCM port, scm_print_state *pstate);
SCM_API int scm_bigprint (SCM exp, SCM port, scm_print_state *pstate);
SCM_API SCM scm_c_locale_stringn_to_number (const char *mem, size_t len,
					    unsigned int radix);
SCM_INTERNAL SCM scm_i_string_to_number (SCM str, unsigned int radix);
SCM_API SCM scm_string_to_number (SCM str, SCM radix);
SCM_API SCM scm_bigequal (SCM x, SCM y);
SCM_API SCM scm_real_equalp (SCM x, SCM y);
SCM_API SCM scm_complex_equalp (SCM x, SCM y);
SCM_API SCM scm_number_p (SCM x);
SCM_API SCM scm_complex_p (SCM x);
SCM_API SCM scm_real_p (SCM x);
SCM_API SCM scm_rational_p (SCM z);
SCM_API SCM scm_integer_p (SCM x);
SCM_API SCM scm_inexact_p (SCM x);
SCM_API SCM scm_num_eq_p (SCM x, SCM y);
SCM_API SCM scm_less_p (SCM x, SCM y);
SCM_API SCM scm_gr_p (SCM x, SCM y);
SCM_API SCM scm_leq_p (SCM x, SCM y);
SCM_API SCM scm_geq_p (SCM x, SCM y);
SCM_API SCM scm_zero_p (SCM z);
SCM_API SCM scm_positive_p (SCM x);
SCM_API SCM scm_negative_p (SCM x);
SCM_API SCM scm_max (SCM x, SCM y);
SCM_API SCM scm_min (SCM x, SCM y);
SCM_API SCM scm_sum (SCM x, SCM y);
SCM_API SCM scm_oneplus (SCM x);
SCM_API SCM scm_difference (SCM x, SCM y);
SCM_API SCM scm_oneminus (SCM x);
SCM_API SCM scm_product (SCM x, SCM y);
SCM_API SCM scm_divide (SCM x, SCM y);
SCM_API SCM scm_floor (SCM x);
SCM_API SCM scm_ceiling (SCM x);
SCM_API double scm_c_truncate (double x);
SCM_API double scm_c_round (double x);
SCM_API SCM scm_truncate_number (SCM x);
SCM_API SCM scm_round_number (SCM x);
SCM_API SCM scm_expt (SCM z1, SCM z2);
SCM_API SCM scm_sin (SCM z);
SCM_API SCM scm_cos (SCM z);
SCM_API SCM scm_tan (SCM z);
SCM_API SCM scm_sinh (SCM z);
SCM_API SCM scm_cosh (SCM z);
SCM_API SCM scm_tanh (SCM z);
SCM_API SCM scm_asin (SCM z);
SCM_API SCM scm_acos (SCM z);
SCM_API SCM scm_atan (SCM x, SCM y);
SCM_API SCM scm_sys_asinh (SCM z);
SCM_API SCM scm_sys_acosh (SCM z);
SCM_API SCM scm_sys_atanh (SCM z);
SCM_API SCM scm_make_rectangular (SCM z1, SCM z2);
SCM_API SCM scm_make_polar (SCM z1, SCM z2);
SCM_API SCM scm_real_part (SCM z);
SCM_API SCM scm_imag_part (SCM z);
SCM_API SCM scm_magnitude (SCM z);
SCM_API SCM scm_angle (SCM z);
SCM_API SCM scm_exact_to_inexact (SCM z);
SCM_API SCM scm_inexact_to_exact (SCM z);
SCM_API SCM scm_trunc (SCM x);
SCM_API SCM scm_log (SCM z);
SCM_API SCM scm_log10 (SCM z);
SCM_API SCM scm_exp (SCM z);
SCM_API SCM scm_sqrt (SCM z);

SCM_INTERNAL SCM scm_i_min (SCM x, SCM y, SCM rest);
SCM_INTERNAL SCM scm_i_max (SCM x, SCM y, SCM rest);
SCM_INTERNAL SCM scm_i_sum (SCM x, SCM y, SCM rest);
SCM_INTERNAL SCM scm_i_difference (SCM x, SCM y, SCM rest);
SCM_INTERNAL SCM scm_i_product (SCM x, SCM y, SCM rest);
SCM_INTERNAL SCM scm_i_divide (SCM x, SCM y, SCM rest);

/* bignum internal functions */
SCM_INTERNAL SCM scm_i_mkbig (void);
SCM_API /* FIXME: not internal */ SCM scm_i_normbig (SCM x);
SCM_INTERNAL int scm_i_bigcmp (SCM a, SCM b);
SCM_INTERNAL SCM scm_i_dbl2big (double d);
SCM_INTERNAL SCM scm_i_dbl2num (double d);
SCM_API /* FIXME: not internal */ double scm_i_big2dbl (SCM b);
SCM_API /* FIXME: not internal */ SCM scm_i_long2big (long n);
SCM_API /* FIXME: not internal */ SCM scm_i_ulong2big (unsigned long n);
SCM_API /* FIXME: not internal */ SCM scm_i_clonebig (SCM src_big, int same_sign_p);

/* ratio functions */
SCM_API SCM scm_rationalize (SCM x, SCM err);
SCM_API SCM scm_numerator (SCM z);
SCM_API SCM scm_denominator (SCM z);

/* fraction internal functions */
SCM_INTERNAL double scm_i_fraction2double (SCM z);
SCM_INTERNAL SCM scm_i_fraction_equalp (SCM x, SCM y);
SCM_INTERNAL int scm_i_print_fraction (SCM sexp, SCM port, scm_print_state *pstate);

/* general internal functions */
SCM_INTERNAL void scm_i_print_double (double val, SCM port);
SCM_INTERNAL void scm_i_print_complex (double real, double imag, SCM port);

/* conversion functions for integers */

SCM_API int scm_is_integer (SCM val);
SCM_API int scm_is_signed_integer (SCM val,
				   scm_t_intmax min, scm_t_intmax max);
SCM_API int scm_is_unsigned_integer (SCM val,
				     scm_t_uintmax min, scm_t_uintmax max);

SCM_API SCM scm_from_signed_integer (scm_t_intmax val);
SCM_API SCM scm_from_unsigned_integer (scm_t_uintmax val);

SCM_API scm_t_intmax scm_to_signed_integer (SCM val,
					    scm_t_intmax min,
					    scm_t_intmax max);
SCM_API scm_t_uintmax scm_to_unsigned_integer (SCM val,
					       scm_t_uintmax min,
					       scm_t_uintmax max);

SCM_API scm_t_int8   scm_to_int8     (SCM x);
SCM_API SCM          scm_from_int8   (scm_t_int8 x);

SCM_API scm_t_uint8  scm_to_uint8    (SCM x);
SCM_API SCM          scm_from_uint8  (scm_t_uint8 x);

SCM_API scm_t_int16  scm_to_int16    (SCM x);
SCM_API SCM          scm_from_int16  (scm_t_int16 x);

SCM_API scm_t_uint16 scm_to_uint16   (SCM x);
SCM_API SCM          scm_from_uint16 (scm_t_uint16 x);

SCM_API scm_t_int32  scm_to_int32    (SCM x);
SCM_API SCM          scm_from_int32  (scm_t_int32 x);

SCM_API scm_t_uint32 scm_to_uint32   (SCM x);
SCM_API SCM          scm_from_uint32 (scm_t_uint32 x);

SCM_API scm_t_wchar  scm_to_wchar    (SCM x);
SCM_API SCM          scm_from_wchar  (scm_t_wchar x);

SCM_API scm_t_int64  scm_to_int64    (SCM x);
SCM_API SCM          scm_from_int64  (scm_t_int64 x);

SCM_API scm_t_uint64 scm_to_uint64   (SCM x);
SCM_API SCM          scm_from_uint64 (scm_t_uint64 x);

SCM_API void scm_to_mpz (SCM x, mpz_t rop);
SCM_API SCM  scm_from_mpz (mpz_t rop);


/* The conversion functions for other types are aliased to the
   appropriate ones from above.  We pick the right one based on the
   size of the type.

   Not each and every possibility is covered by the code below, and
   while it is trivial to complete the tests, it might be better to
   just test for the 'sane' possibilities.  When one of the tests
   below fails, chances are good that some silent assumption somewhere
   else will also fail.
*/

#if SCM_SIZEOF_CHAR == 1
#define scm_to_schar   scm_to_int8
#define scm_from_schar scm_from_int8
#define scm_to_uchar   scm_to_uint8
#define scm_from_uchar scm_from_uint8
#if CHAR_MIN == 0
#define scm_to_char    scm_to_uint8
#define scm_from_char  scm_from_uint8
#else
#define scm_to_char    scm_to_int8
#define scm_from_char  scm_from_int8
#endif
#else
#error sizeof(char) is not 1.
#endif

#if SCM_SIZEOF_SHORT == 1
#define scm_to_short    scm_to_int8
#define scm_from_short  scm_from_int8
#define scm_to_ushort   scm_to_uint8
#define scm_from_ushort scm_from_uint8
#else
#if SCM_SIZEOF_SHORT == 2
#define scm_to_short    scm_to_int16
#define scm_from_short  scm_from_int16
#define scm_to_ushort   scm_to_uint16
#define scm_from_ushort scm_from_uint16
#else
#if SCM_SIZEOF_SHORT == 4
#define scm_to_short    scm_to_int32
#define scm_from_short  scm_from_int32
#define scm_to_ushort   scm_to_uint32
#define scm_from_ushort scm_from_uint32
#else
#error sizeof(short) is not 1, 2, or 4.
#endif
#endif
#endif

#if SCM_SIZEOF_INT == 4
#define scm_to_int    scm_to_int32
#define scm_from_int  scm_from_int32
#define scm_to_uint   scm_to_uint32
#define scm_from_uint scm_from_uint32
#else
#if SCM_SIZEOF_INT == 8
#define scm_to_int    scm_to_int64
#define scm_from_int  scm_from_int64
#define scm_to_uint   scm_to_uint64
#define scm_from_uint scm_from_uint64
#else
#error sizeof(int) is not 4 or 8.
#endif
#endif

#if SCM_SIZEOF_LONG == 4
#define scm_to_long    scm_to_int32
#define scm_from_long  scm_from_int32
#define scm_to_ulong   scm_to_uint32
#define scm_from_ulong scm_from_uint32
#else
#if SCM_SIZEOF_LONG == 8
#define scm_to_long    scm_to_int64
#define scm_from_long  scm_from_int64
#define scm_to_ulong   scm_to_uint64
#define scm_from_ulong scm_from_uint64
#else
#error sizeof(long) is not 4 or 8.
#endif
#endif

#if SCM_SIZEOF_INTMAX == 4
#define scm_to_intmax    scm_to_int32
#define scm_from_intmax  scm_from_int32
#define scm_to_uintmax   scm_to_uint32
#define scm_from_uintmax scm_from_uint32
#else
#if SCM_SIZEOF_INTMAX == 8
#define scm_to_intmax    scm_to_int64
#define scm_from_intmax  scm_from_int64
#define scm_to_uintmax   scm_to_uint64
#define scm_from_uintmax scm_from_uint64
#else
#error sizeof(scm_t_intmax) is not 4 or 8.
#endif
#endif

#if SCM_SIZEOF_LONG_LONG == 0
#else
#if SCM_SIZEOF_LONG_LONG == 8
#define scm_to_long_long    scm_to_int64
#define scm_from_long_long  scm_from_int64
#define scm_to_ulong_long   scm_to_uint64
#define scm_from_ulong_long scm_from_uint64
#else
#error sizeof(long long) is not 8.
#endif
#endif

#if SCM_SIZEOF_SIZE_T == 4
#define scm_to_ssize_t    scm_to_int32
#define scm_from_ssize_t  scm_from_int32
#define scm_to_size_t     scm_to_uint32
#define scm_from_size_t   scm_from_uint32
#else
#if SCM_SIZEOF_SIZE_T == 8
#define scm_to_ssize_t    scm_to_int64
#define scm_from_ssize_t  scm_from_int64
#define scm_to_size_t     scm_to_uint64
#define scm_from_size_t   scm_from_uint64
#else
#error sizeof(size_t) is not 4 or 8.
#endif
#endif

/* conversion functions for double */

SCM_API int scm_is_real (SCM val);
SCM_API int scm_is_rational (SCM val);
SCM_API double scm_to_double (SCM val);
SCM_API SCM scm_from_double (double val);

/* conversion functions for complex */

SCM_API int scm_is_complex (SCM val);
SCM_API SCM scm_c_make_rectangular (double re, double im);
SCM_API SCM scm_c_make_polar (double mag, double ang);
SCM_API double scm_c_real_part (SCM z);
SCM_API double scm_c_imag_part (SCM z);
SCM_API double scm_c_magnitude (SCM z);
SCM_API double scm_c_angle (SCM z);

SCM_API int scm_is_number (SCM val);

SCM_INTERNAL void scm_init_numbers (void);

#endif  /* SCM_NUMBERS_H */

/*
  Local Variables:
  c-file-style: "gnu"
  End:
*/
Commit	Line	Data
0f2d19dd JB	1	/* classes: h_files */
0f2d19dd JB	2
3c9a524f DH	3	#ifndef SCM_NUMBERS_H
3c9a524f DH	4	#define SCM_NUMBERS_H
dee01b01	5
4ca48269	6	/* Copyright (C) 1995,1996,1998,2000,2001,2002,2003,2004,2005, 2006, 2008, 2009, 2010 Free Software Foundation, Inc.
dee01b01	7	*
73be1d9e	8	* This library is free software; you can redistribute it and/or
53befeb7 NJ	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.
dee01b01	12	*
53befeb7 NJ	13	* This library is distributed in the hope that it will be useful, but
53befeb7 NJ	14	* WITHOUT ANY WARRANTY; without even the implied warranty of
73be1d9e MV	15	* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
73be1d9e MV	16	* Lesser General Public License for more details.
dee01b01	17	*
73be1d9e MV	18	* You should have received a copy of the GNU Lesser General Public
73be1d9e MV	19	* License along with this library; if not, write to the Free Software
53befeb7 NJ	20	* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
53befeb7 NJ	21	* 02110-1301 USA
73be1d9e	22	*/
d3a6bc94	23
0f2d19dd JB	24	\f
0f2d19dd JB	25
cd036260 MV	26	#include <gmp.h>
cd036260 MV	27
b4309c3c	28	#include "libguile/__scm.h"
34f74fe1 MD	29	#include "libguile/print.h"
34f74fe1 MD	30
889975e5 MG	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
0f2d19dd JB	36	\f
0f2d19dd JB	37
14282d0f	38	/* Immediate Numbers, also known as fixnums
0f2d19dd	39	*
14282d0f	40	* Inums are exact integer data that fits within an SCM word. */
0f2d19dd	41
c0060c95 MV	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
004c0902	47	#define SCM_I_FIXNUM_BIT (SCM_LONG_BIT - 2)
c0060c95	48	#define SCM_MOST_POSITIVE_FIXNUM ((SCM_T_SIGNED_BITS_MAX-3)/4)
004c0902	49	#define SCM_MOST_NEGATIVE_FIXNUM (-SCM_MOST_POSITIVE_FIXNUM-1)
5c75b29f	50
56100716 DH	51	/* SCM_SRS is signed right shift */
56100716 DH	52	#if (-1 == (((-1) << 2) + 2) >> 2)
14282d0f	53	# define SCM_SRS(x, y) ((x) >> (y))
0f2d19dd	54	#else
14282d0f	55	# define SCM_SRS(x, y) ((x) < 0 ? ~((~(x)) >> (y)) : ((x) >> (y)))
56100716	56	#endif /* (-1 == (((-1) << 2) + 2) >> 2) */
0f2d19dd JB	57
0f2d19dd JB	58
3aa13a05 MV	59	#define SCM_I_INUMP(x) (2 & SCM_UNPACK (x))
3aa13a05 MV	60	#define SCM_I_NINUMP(x) (!SCM_I_INUMP (x))
5a572ca1	61	#define SCM_I_MAKINUM(x) \
5d7d39ff	62	(SCM_PACK ((((scm_t_signed_bits) (x)) << 2) + scm_tc2_int))
3aa13a05	63	#define SCM_I_INUM(x) (SCM_SRS ((scm_t_signed_bits) SCM_UNPACK (x), 2))
0f2d19dd	64
894a712b DH	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)
14282d0f	68	#define SCM_FIXABLE(n) (SCM_POSFIXABLE (n) && SCM_NEGFIXABLE (n))
894a712b DH	69
894a712b DH	70
56100716	71	/* A name for 0. */
5a572ca1	72	#define SCM_INUM0 (SCM_I_MAKINUM (0))
0f2d19dd	73
c0060c95	74	/* SCM_MAXEXP is the maximum double precision exponent
88eb6852 MD	75	* SCM_FLTMAX is less than or scm_equal the largest single precision float
	76	*/
	77
f0ad7455 RB	78	#if SCM_HAVE_STDC_HEADERS
f0ad7455 RB	79	# ifndef GO32
1e498fbd SJ	80	# include <float.h>
	81	# ifdef __MINGW32__
	82	# define copysign _copysign
1e498fbd SJ	83	# define finite _finite
1e498fbd SJ	84	# endif /* __MINGW32__ */
f0ad7455	85	# endif /* ndef GO32 */
88eb6852	86	#endif /* def STDC_HEADERS */
f0ad7455	87
88eb6852	88	#ifdef DBL_MAX_10_EXP
f0ad7455	89	# define SCM_MAXEXP DBL_MAX_10_EXP
88eb6852	90	#else
f0ad7455	91	# define SCM_MAXEXP 308 /* IEEE doubles */
88eb6852	92	#endif /* def DBL_MAX_10_EXP */
f0ad7455	93
88eb6852	94	#ifdef FLT_MAX
f0ad7455	95	# define SCM_FLTMAX FLT_MAX
88eb6852	96	#else
f0ad7455	97	# define SCM_FLTMAX 1e+23
88eb6852 MD	98	#endif /* def FLT_MAX */
88eb6852 MD	99
0f2d19dd	100
2881e77b MV	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.
0f2d19dd	105	*/
2881e77b	106	#define SCM_INTBUFLEN (5 + SCM_CHAR_BIT*sizeof(scm_t_intmax))
0f2d19dd JB	107
	108	\f
	109
	110	/* Numbers
	111	*/
	112
534c55a9 DH	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
f03314f9 DH	118	/* Number subtype 1 to 3 (note the dependency on the predicates SCM_INEXACTP
f03314f9 DH	119	* and SCM_NUMP) */
534c55a9 DH	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)
f92e85f7	123	#define scm_tc16_fraction (scm_tc7_number + 4 * 256L)
534c55a9	124
f03314f9 DH	125	#define SCM_INEXACTP(x) \
f03314f9 DH	126	(!SCM_IMP (x) && (0xfeff & SCM_CELL_TYPE (x)) == scm_tc16_real)
0d5e3480 DH	127	#define SCM_REALP(x) (!SCM_IMP (x) && SCM_TYP16 (x) == scm_tc16_real)
0d5e3480 DH	128	#define SCM_COMPLEXP(x) (!SCM_IMP (x) && SCM_TYP16 (x) == scm_tc16_complex)
88eb6852	129
92c2555f	130	#define SCM_REAL_VALUE(x) (((scm_t_double *) SCM2PTR (x))->real)
03604fcf LC	131	#define SCM_COMPLEX_REAL(x) (((scm_t_complex *) SCM2PTR (x))->real)
03604fcf LC	132	#define SCM_COMPLEX_IMAG(x) (((scm_t_complex *) SCM2PTR (x))->imag)
0f2d19dd	133
d3e7e88d	134	/* Each bignum is just an mpz_t stored in a double cell starting at word 1. */
fba8fdc0	135	#define SCM_I_BIG_MPZ(x) (((mpz_t ) (SCM_CELL_OBJECT_LOC((x),1))))
d3e7e88d	136	#define SCM_BIGP(x) (!SCM_IMP (x) && SCM_TYP16 (x) == scm_tc16_big)
0f2d19dd	137
3aa13a05	138	#define SCM_NUMBERP(x) (SCM_I_INUMP(x) \|\| SCM_NUMP(x))
d3e7e88d	139	#define SCM_NUMP(x) (!SCM_IMP(x) \
f92e85f7 MV	140	&& (((0xfcff & SCM_CELL_TYPE (x)) == scm_tc7_number) \
	141	\|\| ((0xfbff & SCM_CELL_TYPE (x)) == scm_tc7_number)))
	142	/* 0xfcff (#b1100) for 0 free, 1 big, 2 real, 3 complex, then 0xfbff (#b1011) for 4 fraction */
	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))
651f07f8 DH	146	#define SCM_FRACTION_DENOMINATOR(x) (SCM_CELL_OBJECT_2 (x))
88eb6852	147
0f2d19dd JB	148	\f
0f2d19dd JB	149
92c2555f	150	typedef struct scm_t_double
0f2d19dd JB	151	{
0f2d19dd JB	152	SCM type;
88eb6852 MD	153	SCM pad;
88eb6852 MD	154	double real;
92c2555f	155	} scm_t_double;
0f2d19dd	156
92c2555f	157	typedef struct scm_t_complex
0f2d19dd	158	{
03604fcf LC	159	SCM type;
03604fcf LC	160	SCM pad;
88eb6852 MD	161	double real;
88eb6852 MD	162	double imag;
92c2555f	163	} scm_t_complex;
0f2d19dd	164
904a78f1	165
0f2d19dd	166	\f
0f2d19dd	167
33b001fd MV	168	SCM_API SCM scm_exact_p (SCM x);
	169	SCM_API SCM scm_odd_p (SCM n);
	170	SCM_API SCM scm_even_p (SCM n);
232898d9 MV	171	SCM_API SCM scm_inf_p (SCM n);
	172	SCM_API SCM scm_nan_p (SCM n);
	173	SCM_API SCM scm_inf (void);
	174	SCM_API SCM scm_nan (void);
33b001fd MV	175	SCM_API SCM scm_abs (SCM x);
	176	SCM_API SCM scm_quotient (SCM x, SCM y);
	177	SCM_API SCM scm_remainder (SCM x, SCM y);
	178	SCM_API SCM scm_modulo (SCM x, SCM y);
	179	SCM_API SCM scm_gcd (SCM x, SCM y);
	180	SCM_API SCM scm_lcm (SCM n1, SCM n2);
	181	SCM_API SCM scm_logand (SCM n1, SCM n2);
	182	SCM_API SCM scm_logior (SCM n1, SCM n2);
	183	SCM_API SCM scm_logxor (SCM n1, SCM n2);
	184	SCM_API SCM scm_logtest (SCM n1, SCM n2);
	185	SCM_API SCM scm_logbit_p (SCM n1, SCM n2);
	186	SCM_API SCM scm_lognot (SCM n);
d885e204	187	SCM_API SCM scm_modulo_expt (SCM n, SCM k, SCM m);
33b001fd MV	188	SCM_API SCM scm_integer_expt (SCM z1, SCM z2);
	189	SCM_API SCM scm_ash (SCM n, SCM cnt);
	190	SCM_API SCM scm_bit_extract (SCM n, SCM start, SCM end);
	191	SCM_API SCM scm_logcount (SCM n);
	192	SCM_API SCM scm_integer_length (SCM n);
1be6b49c	193
78d3deb1 AW	194	SCM_INTERNAL SCM scm_i_gcd (SCM x, SCM y, SCM rest);
	195	SCM_INTERNAL SCM scm_i_lcm (SCM x, SCM y, SCM rest);
	196	SCM_INTERNAL SCM scm_i_logand (SCM x, SCM y, SCM rest);
	197	SCM_INTERNAL SCM scm_i_logior (SCM x, SCM y, SCM rest);
	198	SCM_INTERNAL SCM scm_i_logxor (SCM x, SCM y, SCM rest);
	199
2881e77b MV	200	SCM_API size_t scm_iint2str (scm_t_intmax num, int rad, char *p);
2881e77b MV	201	SCM_API size_t scm_iuint2str (scm_t_uintmax num, int rad, char *p);
33b001fd MV	202	SCM_API SCM scm_number_to_string (SCM x, SCM radix);
	203	SCM_API int scm_print_real (SCM sexp, SCM port, scm_print_state *pstate);
	204	SCM_API int scm_print_complex (SCM sexp, SCM port, scm_print_state *pstate);
	205	SCM_API int scm_bigprint (SCM exp, SCM port, scm_print_state *pstate);
2b829bbb KR	206	SCM_API SCM scm_c_locale_stringn_to_number (const char *mem, size_t len,
2b829bbb KR	207	unsigned int radix);
3f47e526	208	SCM_INTERNAL SCM scm_i_string_to_number (SCM str, unsigned int radix);
33b001fd	209	SCM_API SCM scm_string_to_number (SCM str, SCM radix);
33b001fd MV	210	SCM_API SCM scm_bigequal (SCM x, SCM y);
	211	SCM_API SCM scm_real_equalp (SCM x, SCM y);
	212	SCM_API SCM scm_complex_equalp (SCM x, SCM y);
	213	SCM_API SCM scm_number_p (SCM x);
8507ec80	214	SCM_API SCM scm_complex_p (SCM x);
33b001fd	215	SCM_API SCM scm_real_p (SCM x);
8507ec80	216	SCM_API SCM scm_rational_p (SCM z);
33b001fd MV	217	SCM_API SCM scm_integer_p (SCM x);
	218	SCM_API SCM scm_inexact_p (SCM x);
	219	SCM_API SCM scm_num_eq_p (SCM x, SCM y);
	220	SCM_API SCM scm_less_p (SCM x, SCM y);
	221	SCM_API SCM scm_gr_p (SCM x, SCM y);
	222	SCM_API SCM scm_leq_p (SCM x, SCM y);
	223	SCM_API SCM scm_geq_p (SCM x, SCM y);
	224	SCM_API SCM scm_zero_p (SCM z);
	225	SCM_API SCM scm_positive_p (SCM x);
	226	SCM_API SCM scm_negative_p (SCM x);
	227	SCM_API SCM scm_max (SCM x, SCM y);
	228	SCM_API SCM scm_min (SCM x, SCM y);
	229	SCM_API SCM scm_sum (SCM x, SCM y);
40882e3d	230	SCM_API SCM scm_oneplus (SCM x);
33b001fd	231	SCM_API SCM scm_difference (SCM x, SCM y);
40882e3d	232	SCM_API SCM scm_oneminus (SCM x);
33b001fd	233	SCM_API SCM scm_product (SCM x, SCM y);
33b001fd	234	SCM_API SCM scm_divide (SCM x, SCM y);
f92e85f7 MV	235	SCM_API SCM scm_floor (SCM x);
f92e85f7 MV	236	SCM_API SCM scm_ceiling (SCM x);
3101f40f MV	237	SCM_API double scm_c_truncate (double x);
3101f40f MV	238	SCM_API double scm_c_round (double x);
f92e85f7 MV	239	SCM_API SCM scm_truncate_number (SCM x);
f92e85f7 MV	240	SCM_API SCM scm_round_number (SCM x);
6fc4d012	241	SCM_API SCM scm_expt (SCM z1, SCM z2);
ad79736c AW	242	SCM_API SCM scm_sin (SCM z);
	243	SCM_API SCM scm_cos (SCM z);
	244	SCM_API SCM scm_tan (SCM z);
	245	SCM_API SCM scm_sinh (SCM z);
	246	SCM_API SCM scm_cosh (SCM z);
	247	SCM_API SCM scm_tanh (SCM z);
	248	SCM_API SCM scm_asin (SCM z);
	249	SCM_API SCM scm_acos (SCM z);
	250	SCM_API SCM scm_atan (SCM x, SCM y);
	251	SCM_API SCM scm_sys_asinh (SCM z);
	252	SCM_API SCM scm_sys_acosh (SCM z);
	253	SCM_API SCM scm_sys_atanh (SCM z);
33b001fd MV	254	SCM_API SCM scm_make_rectangular (SCM z1, SCM z2);
	255	SCM_API SCM scm_make_polar (SCM z1, SCM z2);
	256	SCM_API SCM scm_real_part (SCM z);
	257	SCM_API SCM scm_imag_part (SCM z);
	258	SCM_API SCM scm_magnitude (SCM z);
	259	SCM_API SCM scm_angle (SCM z);
	260	SCM_API SCM scm_exact_to_inexact (SCM z);
	261	SCM_API SCM scm_inexact_to_exact (SCM z);
	262	SCM_API SCM scm_trunc (SCM x);
8ab3d8a0 KR	263	SCM_API SCM scm_log (SCM z);
	264	SCM_API SCM scm_log10 (SCM z);
	265	SCM_API SCM scm_exp (SCM z);
	266	SCM_API SCM scm_sqrt (SCM z);
33b001fd	267
78d3deb1 AW	268	SCM_INTERNAL SCM scm_i_min (SCM x, SCM y, SCM rest);
	269	SCM_INTERNAL SCM scm_i_max (SCM x, SCM y, SCM rest);
	270	SCM_INTERNAL SCM scm_i_sum (SCM x, SCM y, SCM rest);
	271	SCM_INTERNAL SCM scm_i_difference (SCM x, SCM y, SCM rest);
	272	SCM_INTERNAL SCM scm_i_product (SCM x, SCM y, SCM rest);
	273	SCM_INTERNAL SCM scm_i_divide (SCM x, SCM y, SCM rest);
	274
d3e7e88d	275	/* bignum internal functions */
102dbb6f LC	276	SCM_INTERNAL SCM scm_i_mkbig (void);
	277	SCM_API /* FIXME: not internal */ SCM scm_i_normbig (SCM x);
	278	SCM_INTERNAL int scm_i_bigcmp (SCM a, SCM b);
	279	SCM_INTERNAL SCM scm_i_dbl2big (double d);
	280	SCM_INTERNAL SCM scm_i_dbl2num (double d);
	281	SCM_API /* FIXME: not internal */ double scm_i_big2dbl (SCM b);
	282	SCM_API /* FIXME: not internal */ SCM scm_i_long2big (long n);
	283	SCM_API /* FIXME: not internal */ SCM scm_i_ulong2big (unsigned long n);
	284	SCM_API /* FIXME: not internal */ SCM scm_i_clonebig (SCM src_big, int same_sign_p);
d3e7e88d	285
f92e85f7	286	/* ratio functions */
f92e85f7 MV	287	SCM_API SCM scm_rationalize (SCM x, SCM err);
	288	SCM_API SCM scm_numerator (SCM z);
	289	SCM_API SCM scm_denominator (SCM z);
f92e85f7 MV	290
f92e85f7 MV	291	/* fraction internal functions */
102dbb6f LC	292	SCM_INTERNAL double scm_i_fraction2double (SCM z);
	293	SCM_INTERNAL SCM scm_i_fraction_equalp (SCM x, SCM y);
	294	SCM_INTERNAL int scm_i_print_fraction (SCM sexp, SCM port, scm_print_state *pstate);
f92e85f7	295
b479fe9a	296	/* general internal functions */
102dbb6f LC	297	SCM_INTERNAL void scm_i_print_double (double val, SCM port);
102dbb6f LC	298	SCM_INTERNAL void scm_i_print_complex (double real, double imag, SCM port);
b479fe9a	299
d1b5b400	300	/* conversion functions for integers */
73e4de09 MV	301
	302	SCM_API int scm_is_integer (SCM val);
	303	SCM_API int scm_is_signed_integer (SCM val,
	304	scm_t_intmax min, scm_t_intmax max);
	305	SCM_API int scm_is_unsigned_integer (SCM val,
	306	scm_t_uintmax min, scm_t_uintmax max);
	307
	308	SCM_API SCM scm_from_signed_integer (scm_t_intmax val);
	309	SCM_API SCM scm_from_unsigned_integer (scm_t_uintmax val);
	310
	311	SCM_API scm_t_intmax scm_to_signed_integer (SCM val,
	312	scm_t_intmax min,
	313	scm_t_intmax max);
	314	SCM_API scm_t_uintmax scm_to_unsigned_integer (SCM val,
	315	scm_t_uintmax min,
	316	scm_t_uintmax max);
	317
bfd7932e MV	318	SCM_API scm_t_int8 scm_to_int8 (SCM x);
	319	SCM_API SCM scm_from_int8 (scm_t_int8 x);
	320
	321	SCM_API scm_t_uint8 scm_to_uint8 (SCM x);
	322	SCM_API SCM scm_from_uint8 (scm_t_uint8 x);
	323
	324	SCM_API scm_t_int16 scm_to_int16 (SCM x);
	325	SCM_API SCM scm_from_int16 (scm_t_int16 x);
	326
	327	SCM_API scm_t_uint16 scm_to_uint16 (SCM x);
	328	SCM_API SCM scm_from_uint16 (scm_t_uint16 x);
	329
	330	SCM_API scm_t_int32 scm_to_int32 (SCM x);
	331	SCM_API SCM scm_from_int32 (scm_t_int32 x);
	332
	333	SCM_API scm_t_uint32 scm_to_uint32 (SCM x);
	334	SCM_API SCM scm_from_uint32 (scm_t_uint32 x);
	335
904a78f1 MG	336	SCM_API scm_t_wchar scm_to_wchar (SCM x);
	337	SCM_API SCM scm_from_wchar (scm_t_wchar x);
	338
bfd7932e MV	339	SCM_API scm_t_int64 scm_to_int64 (SCM x);
	340	SCM_API SCM scm_from_int64 (scm_t_int64 x);
	341
	342	SCM_API scm_t_uint64 scm_to_uint64 (SCM x);
	343	SCM_API SCM scm_from_uint64 (scm_t_uint64 x);
	344
cd036260 MV	345	SCM_API void scm_to_mpz (SCM x, mpz_t rop);
	346	SCM_API SCM scm_from_mpz (mpz_t rop);
	347
	348
d1b5b400 MV	349	/* The conversion functions for other types are aliased to the
	350	appropriate ones from above. We pick the right one based on the
	351	size of the type.
	352
	353	Not each and every possibility is covered by the code below, and
	354	while it is trivial to complete the tests, it might be better to
	355	just test for the 'sane' possibilities. When one of the tests
	356	below fails, chances are good that some silent assumption somewhere
	357	else will also fail.
	358	*/
	359
	360	#if SCM_SIZEOF_CHAR == 1
	361	#define scm_to_schar scm_to_int8
	362	#define scm_from_schar scm_from_int8
	363	#define scm_to_uchar scm_to_uint8
	364	#define scm_from_uchar scm_from_uint8
73e4de09	365	#if CHAR_MIN == 0
d1b5b400 MV	366	#define scm_to_char scm_to_uint8
d1b5b400 MV	367	#define scm_from_char scm_from_uint8
73e4de09	368	#else
d1b5b400 MV	369	#define scm_to_char scm_to_int8
d1b5b400 MV	370	#define scm_from_char scm_from_int8
73e4de09	371	#endif
d1b5b400 MV	372	#else
d1b5b400 MV	373	#error sizeof(char) is not 1.
22008a78 MV	374	#endif
22008a78 MV	375
d1b5b400 MV	376	#if SCM_SIZEOF_SHORT == 1
	377	#define scm_to_short scm_to_int8
	378	#define scm_from_short scm_from_int8
	379	#define scm_to_ushort scm_to_uint8
	380	#define scm_from_ushort scm_from_uint8
	381	#else
	382	#if SCM_SIZEOF_SHORT == 2
	383	#define scm_to_short scm_to_int16
	384	#define scm_from_short scm_from_int16
	385	#define scm_to_ushort scm_to_uint16
	386	#define scm_from_ushort scm_from_uint16
	387	#else
	388	#if SCM_SIZEOF_SHORT == 4
	389	#define scm_to_short scm_to_int32
	390	#define scm_from_short scm_from_int32
	391	#define scm_to_ushort scm_to_uint32
	392	#define scm_from_ushort scm_from_uint32
73e4de09	393	#else
d1b5b400 MV	394	#error sizeof(short) is not 1, 2, or 4.
	395	#endif
	396	#endif
73e4de09 MV	397	#endif
73e4de09 MV	398
d1b5b400 MV	399	#if SCM_SIZEOF_INT == 4
	400	#define scm_to_int scm_to_int32
	401	#define scm_from_int scm_from_int32
	402	#define scm_to_uint scm_to_uint32
	403	#define scm_from_uint scm_from_uint32
	404	#else
	405	#if SCM_SIZEOF_INT == 8
	406	#define scm_to_int scm_to_int64
	407	#define scm_from_int scm_from_int64
	408	#define scm_to_uint scm_to_uint64
	409	#define scm_from_uint scm_from_uint64
	410	#else
	411	#error sizeof(int) is not 4 or 8.
	412	#endif
	413	#endif
73e4de09	414
d1b5b400 MV	415	#if SCM_SIZEOF_LONG == 4
	416	#define scm_to_long scm_to_int32
	417	#define scm_from_long scm_from_int32
	418	#define scm_to_ulong scm_to_uint32
	419	#define scm_from_ulong scm_from_uint32
	420	#else
	421	#if SCM_SIZEOF_LONG == 8
	422	#define scm_to_long scm_to_int64
	423	#define scm_from_long scm_from_int64
	424	#define scm_to_ulong scm_to_uint64
	425	#define scm_from_ulong scm_from_uint64
	426	#else
	427	#error sizeof(long) is not 4 or 8.
	428	#endif
	429	#endif
73e4de09	430
d1b5b400 MV	431	#if SCM_SIZEOF_INTMAX == 4
	432	#define scm_to_intmax scm_to_int32
	433	#define scm_from_intmax scm_from_int32
	434	#define scm_to_uintmax scm_to_uint32
	435	#define scm_from_uintmax scm_from_uint32
	436	#else
	437	#if SCM_SIZEOF_INTMAX == 8
	438	#define scm_to_intmax scm_to_int64
	439	#define scm_from_intmax scm_from_int64
	440	#define scm_to_uintmax scm_to_uint64
	441	#define scm_from_uintmax scm_from_uint64
	442	#else
55f26379	443	#error sizeof(scm_t_intmax) is not 4 or 8.
d1b5b400 MV	444	#endif
d1b5b400 MV	445	#endif
73e4de09	446
d1b5b400 MV	447	#if SCM_SIZEOF_LONG_LONG == 0
	448	#else
	449	#if SCM_SIZEOF_LONG_LONG == 8
	450	#define scm_to_long_long scm_to_int64
	451	#define scm_from_long_long scm_from_int64
	452	#define scm_to_ulong_long scm_to_uint64
	453	#define scm_from_ulong_long scm_from_uint64
	454	#else
	455	#error sizeof(long long) is not 8.
	456	#endif
	457	#endif
73e4de09	458
d1b5b400 MV	459	#if SCM_SIZEOF_SIZE_T == 4
	460	#define scm_to_ssize_t scm_to_int32
	461	#define scm_from_ssize_t scm_from_int32
	462	#define scm_to_size_t scm_to_uint32
	463	#define scm_from_size_t scm_from_uint32
	464	#else
	465	#if SCM_SIZEOF_SIZE_T == 8
	466	#define scm_to_ssize_t scm_to_int64
	467	#define scm_from_ssize_t scm_from_int64
	468	#define scm_to_size_t scm_to_uint64
	469	#define scm_from_size_t scm_from_uint64
	470	#else
	471	#error sizeof(size_t) is not 4 or 8.
	472	#endif
22008a78 MV	473	#endif
22008a78 MV	474
8507ec80	475	/* conversion functions for double */
22008a78	476
73e4de09	477	SCM_API int scm_is_real (SCM val);
55f26379	478	SCM_API int scm_is_rational (SCM val);
73e4de09 MV	479	SCM_API double scm_to_double (SCM val);
	480	SCM_API SCM scm_from_double (double val);
	481
8507ec80 MV	482	/* conversion functions for complex */
	483
	484	SCM_API int scm_is_complex (SCM val);
	485	SCM_API SCM scm_c_make_rectangular (double re, double im);
	486	SCM_API SCM scm_c_make_polar (double mag, double ang);
	487	SCM_API double scm_c_real_part (SCM z);
	488	SCM_API double scm_c_imag_part (SCM z);
	489	SCM_API double scm_c_magnitude (SCM z);
	490	SCM_API double scm_c_angle (SCM z);
	491
	492	SCM_API int scm_is_number (SCM val);
	493
102dbb6f	494	SCM_INTERNAL void scm_init_numbers (void);
0f2d19dd	495
3c9a524f	496	#endif /* SCM_NUMBERS_H */
89e00824 ML	497
	498	/*
	499	Local Variables:
	500	c-file-style: "gnu"
	501	End:
	502	*/