[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 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 2.1 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., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
 */

\f

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

#if SCM_HAVE_FLOATINGPOINT_H
# include <floatingpoint.h>
#endif

#if SCM_HAVE_IEEEFP_H
# include <ieeefp.h>
#endif

#if SCM_HAVE_NAN_H
# if defined (SCO)
#   define _IEEE 1
# endif
# include <nan.h>
# if defined (SCO)
#   undef _IEEE
# endif
#endif /* SCM_HAVE_NAN_H */

\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_INUMP(x)	(2 & SCM_UNPACK (x))
#define SCM_NINUMP(x) 	(!SCM_INUMP (x))
#define SCM_MAKINUM(x) \
  (SCM_PACK ((((scm_t_signed_bits) (x)) << 2) + scm_tc2_int))
#define SCM_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_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 isnan _isnan
#   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 exact immediate.
 */
#define SCM_INTBUFLEN (5 + SCM_LONG_BIT)

\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_MEM(x) ((scm_t_complex *) SCM_CELL_WORD_1 (x))
#define SCM_COMPLEX_REAL(x) (SCM_COMPLEX_MEM (x)->real)
#define SCM_COMPLEX_IMAG(x) (SCM_COMPLEX_MEM (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_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_SLOPPY_FRACTIONP(x) (SCM_TYP16 (x) == scm_tc16_fraction)
#define SCM_FRACTION_NUMERATOR(x)   ((SCM) (SCM_CELL_WORD_1 (x)))
#define SCM_FRACTION_DENOMINATOR(x) ((SCM) (SCM_CELL_WORD_2 (x)))
#define SCM_FRACTION_SET_NUMERATOR(x, v)   ((SCM) (SCM_SET_CELL_WORD_1 ((x), (v))))
#define SCM_FRACTION_SET_DENOMINATOR(x, v) ((SCM) (SCM_SET_CELL_WORD_2 ((x), (v))))
  /* I think the left half word is free in the type, so I'll use bit 17 */
#define SCM_FRACTION_REDUCED_BIT 0x10000
#define SCM_FRACTION_REDUCED_SET(x) (SCM_SET_CELL_TYPE((x), (SCM_CELL_TYPE (x) | SCM_FRACTION_REDUCED_BIT)))
#define SCM_FRACTION_REDUCED_CLEAR(x) (SCM_SET_CELL_TYPE((x), (SCM_CELL_TYPE (x) & ~SCM_FRACTION_REDUCED_BIT)))
#define SCM_FRACTION_REDUCED(x)     (0x10000 & SCM_CELL_TYPE (x))

\f

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

typedef struct scm_t_complex
{
  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_API size_t scm_iint2str (long 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_i_mem2number (const char *mem, size_t len, unsigned int radix);
SCM_API SCM scm_string_to_number (SCM str, SCM radix);
SCM_API SCM scm_make_real (double x);
SCM_API SCM scm_make_complex (double x, double y);
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_real_p (SCM x);
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_difference (SCM x, SCM y);
SCM_API SCM scm_product (SCM x, SCM y);
SCM_API double scm_num2dbl (SCM a, const char * why);
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_asinh (double x);
SCM_API double scm_acosh (double x);
SCM_API double scm_atanh (double x);
SCM_API double scm_truncate (double x);
SCM_API double scm_round (double x);
SCM_API SCM scm_truncate_number (SCM x);
SCM_API SCM scm_round_number (SCM x);
SCM_API SCM scm_sys_expt (SCM z1, SCM z2);
SCM_API SCM scm_sys_atan2 (SCM z1, SCM z2);
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_short2num (short n);
SCM_API SCM scm_ushort2num (unsigned short n);
SCM_API SCM scm_int2num (int n);
SCM_API SCM scm_uint2num (unsigned int n);
SCM_API SCM scm_long2num (long n);
SCM_API SCM scm_ulong2num (unsigned long n);
SCM_API SCM scm_size2num (size_t n);
SCM_API SCM scm_ptrdiff2num (scm_t_ptrdiff n);
SCM_API short scm_num2short (SCM num, unsigned long int pos,
			     const char *s_caller);
SCM_API unsigned short scm_num2ushort (SCM num, unsigned long int pos,
				       const char *s_caller);
SCM_API int scm_num2int (SCM num, unsigned long int pos,
			 const char *s_caller);
SCM_API unsigned int scm_num2uint (SCM num, unsigned long int pos,
				   const char *s_caller);
SCM_API long scm_num2long (SCM num, unsigned long int pos,
			   const char *s_caller);
SCM_API unsigned long scm_num2ulong (SCM num, unsigned long int pos,
				     const char *s_caller);
SCM_API scm_t_ptrdiff scm_num2ptrdiff (SCM num, unsigned long int pos,
                                       const char *s_caller);
SCM_API size_t scm_num2size (SCM num, unsigned long int pos,
			     const char *s_caller);
#if SCM_SIZEOF_LONG_LONG != 0
SCM_API SCM scm_long_long2num (long long sl);
SCM_API SCM scm_ulong_long2num (unsigned long long sl);
SCM_API long long scm_num2long_long (SCM num, unsigned long int pos,
				     const char *s_caller);
SCM_API unsigned long long scm_num2ulong_long (SCM num, unsigned long int pos,
					       const char *s_caller);
#endif

SCM_API SCM scm_float2num (float n);
SCM_API SCM scm_double2num (double n);
SCM_API float scm_num2float (SCM num, unsigned long int pos,
			     const char *s_caller);
SCM_API double scm_num2double (SCM num, unsigned long int pos,
			       const char *s_caller);


/* bignum internal functions */
SCM_API SCM scm_i_mkbig (void);
SCM_API SCM scm_i_normbig (SCM x);
SCM_API int scm_i_bigcmp (SCM a, SCM b);
SCM_API SCM scm_i_dbl2big (double d);
SCM_API SCM scm_i_dbl2num (double d);
SCM_API double scm_i_big2dbl (SCM b);
SCM_API SCM scm_i_short2big (short n);
SCM_API SCM scm_i_ushort2big (unsigned short n);
SCM_API SCM scm_i_int2big (int n);
SCM_API SCM scm_i_uint2big (unsigned int n);
SCM_API SCM scm_i_long2big (long n);
SCM_API SCM scm_i_ulong2big (unsigned long n);
SCM_API SCM scm_i_size2big (size_t n);
SCM_API SCM scm_i_ptrdiff2big (scm_t_ptrdiff n);

#if SCM_SIZEOF_LONG_LONG != 0
SCM_API SCM scm_i_long_long2big (long long n);
SCM_API SCM scm_i_ulong_long2big (unsigned long long n);
#endif


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

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


#ifdef GUILE_DEBUG
SCM_API SCM scm_sys_check_number_conversions (void);
#endif

SCM_API 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
32c91687	6	/* Copyright (C) 1995,1996,1998,2000,2001,2002,2003,2004 Free Software Foundation, Inc.
dee01b01	7	*
73be1d9e MV	8	* This library is free software; you can redistribute it and/or
	9	* modify it under the terms of the GNU Lesser General Public
	10	* License as published by the Free Software Foundation; either
	11	* version 2.1 of the License, or (at your option) any later version.
dee01b01	12	*
73be1d9e	13	* This library is distributed in the hope that it will be useful,
0f2d19dd	14	* but 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
	19	* License along with this library; if not, write to the Free Software
	20	* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
	21	*/
d3a6bc94	22
0f2d19dd JB	23	\f
0f2d19dd JB	24
b4309c3c	25	#include "libguile/__scm.h"
34f74fe1 MD	26	#include "libguile/print.h"
34f74fe1 MD	27
f0ad7455 RB	28	#if SCM_HAVE_FLOATINGPOINT_H
f0ad7455 RB	29	# include <floatingpoint.h>
232898d9 MV	30	#endif
232898d9 MV	31
f0ad7455 RB	32	#if SCM_HAVE_IEEEFP_H
f0ad7455 RB	33	# include <ieeefp.h>
232898d9 MV	34	#endif
232898d9 MV	35
f0ad7455 RB	36	#if SCM_HAVE_NAN_H
	37	# if defined (SCO)
	38	# define _IEEE 1
	39	# endif
	40	# include <nan.h>
	41	# if defined (SCO)
	42	# undef _IEEE
	43	# endif
	44	#endif /* SCM_HAVE_NAN_H */
232898d9	45
0f2d19dd JB	46	\f
0f2d19dd JB	47
14282d0f	48	/* Immediate Numbers, also known as fixnums
0f2d19dd	49	*
14282d0f	50	* Inums are exact integer data that fits within an SCM word. */
0f2d19dd	51
c0060c95 MV	52	/* SCM_T_SIGNED_MAX is (- (expt 2 n) 1),
	53	* SCM_MOST_POSITIVE_FIXNUM should be (- (expt 2 (- n 2)) 1)
	54	* which is the same as (/ (- (expt 2 n) 4) 4)
	55	*/
	56
004c0902	57	#define SCM_I_FIXNUM_BIT (SCM_LONG_BIT - 2)
c0060c95	58	#define SCM_MOST_POSITIVE_FIXNUM ((SCM_T_SIGNED_BITS_MAX-3)/4)
004c0902	59	#define SCM_MOST_NEGATIVE_FIXNUM (-SCM_MOST_POSITIVE_FIXNUM-1)
5c75b29f	60
56100716 DH	61	/* SCM_SRS is signed right shift */
56100716 DH	62	#if (-1 == (((-1) << 2) + 2) >> 2)
14282d0f	63	# define SCM_SRS(x, y) ((x) >> (y))
0f2d19dd	64	#else
14282d0f	65	# define SCM_SRS(x, y) ((x) < 0 ? ~((~(x)) >> (y)) : ((x) >> (y)))
56100716	66	#endif /* (-1 == (((-1) << 2) + 2) >> 2) */
0f2d19dd JB	67
0f2d19dd JB	68
56100716 DH	69	#define SCM_INUMP(x) (2 & SCM_UNPACK (x))
56100716 DH	70	#define SCM_NINUMP(x) (!SCM_INUMP (x))
5d7d39ff DH	71	#define SCM_MAKINUM(x) \
5d7d39ff DH	72	(SCM_PACK ((((scm_t_signed_bits) (x)) << 2) + scm_tc2_int))
14282d0f	73	#define SCM_INUM(x) (SCM_SRS ((scm_t_signed_bits) SCM_UNPACK (x), 2))
0f2d19dd JB	74
0f2d19dd JB	75
894a712b DH	76	/* SCM_FIXABLE is true if its long argument can be encoded in an SCM_INUM. */
	77	#define SCM_POSFIXABLE(n) ((n) <= SCM_MOST_POSITIVE_FIXNUM)
	78	#define SCM_NEGFIXABLE(n) ((n) >= SCM_MOST_NEGATIVE_FIXNUM)
14282d0f	79	#define SCM_FIXABLE(n) (SCM_POSFIXABLE (n) && SCM_NEGFIXABLE (n))
894a712b DH	80
894a712b DH	81
56100716	82	/* A name for 0. */
b82c6ce0	83	#define SCM_INUM0 (SCM_MAKINUM (0))
0f2d19dd JB	84
0f2d19dd JB	85
c0060c95	86	/* SCM_MAXEXP is the maximum double precision exponent
88eb6852 MD	87	* SCM_FLTMAX is less than or scm_equal the largest single precision float
	88	*/
	89
f0ad7455 RB	90	#if SCM_HAVE_STDC_HEADERS
f0ad7455 RB	91	# ifndef GO32
1e498fbd SJ	92	# include <float.h>
	93	# ifdef __MINGW32__
	94	# define copysign _copysign
	95	# define isnan _isnan
	96	# define finite _finite
	97	# endif /* __MINGW32__ */
f0ad7455	98	# endif /* ndef GO32 */
88eb6852	99	#endif /* def STDC_HEADERS */
f0ad7455	100
88eb6852	101	#ifdef DBL_MAX_10_EXP
f0ad7455	102	# define SCM_MAXEXP DBL_MAX_10_EXP
88eb6852	103	#else
f0ad7455	104	# define SCM_MAXEXP 308 /* IEEE doubles */
88eb6852	105	#endif /* def DBL_MAX_10_EXP */
f0ad7455	106
88eb6852	107	#ifdef FLT_MAX
f0ad7455	108	# define SCM_FLTMAX FLT_MAX
88eb6852	109	#else
f0ad7455	110	# define SCM_FLTMAX 1e+23
88eb6852 MD	111	#endif /* def FLT_MAX */
88eb6852 MD	112
0f2d19dd	113
0f2d19dd JB	114	/* SCM_INTBUFLEN is the maximum number of characters neccessary for the
	115	* printed or scm_string representation of an exact immediate.
	116	*/
c014a02e	117	#define SCM_INTBUFLEN (5 + SCM_LONG_BIT)
0f2d19dd JB	118
	119	\f
	120
	121	/* Numbers
	122	*/
	123
534c55a9 DH	124
	125	/* Note that scm_tc16_real and scm_tc16_complex are given tc16-codes that only
	126	* differ in one bit: This way, checking if an object is an inexact number can
	127	* be done quickly (using the TYP16S macro). */
	128
f03314f9 DH	129	/* Number subtype 1 to 3 (note the dependency on the predicates SCM_INEXACTP
f03314f9 DH	130	* and SCM_NUMP) */
534c55a9 DH	131	#define scm_tc16_big (scm_tc7_number + 1 * 256L)
	132	#define scm_tc16_real (scm_tc7_number + 2 * 256L)
	133	#define scm_tc16_complex (scm_tc7_number + 3 * 256L)
f92e85f7	134	#define scm_tc16_fraction (scm_tc7_number + 4 * 256L)
534c55a9	135
f03314f9 DH	136	#define SCM_INEXACTP(x) \
f03314f9 DH	137	(!SCM_IMP (x) && (0xfeff & SCM_CELL_TYPE (x)) == scm_tc16_real)
0d5e3480 DH	138	#define SCM_REALP(x) (!SCM_IMP (x) && SCM_TYP16 (x) == scm_tc16_real)
0d5e3480 DH	139	#define SCM_COMPLEXP(x) (!SCM_IMP (x) && SCM_TYP16 (x) == scm_tc16_complex)
88eb6852	140
92c2555f MV	141	#define SCM_REAL_VALUE(x) (((scm_t_double *) SCM2PTR (x))->real)
92c2555f MV	142	#define SCM_COMPLEX_MEM(x) ((scm_t_complex *) SCM_CELL_WORD_1 (x))
405aaef9 DH	143	#define SCM_COMPLEX_REAL(x) (SCM_COMPLEX_MEM (x)->real)
405aaef9 DH	144	#define SCM_COMPLEX_IMAG(x) (SCM_COMPLEX_MEM (x)->imag)
0f2d19dd	145
d3e7e88d	146	/* Each bignum is just an mpz_t stored in a double cell starting at word 1. */
fba8fdc0	147	#define SCM_I_BIG_MPZ(x) (((mpz_t ) (SCM_CELL_OBJECT_LOC((x),1))))
d3e7e88d	148	#define SCM_BIGP(x) (!SCM_IMP (x) && SCM_TYP16 (x) == scm_tc16_big)
0f2d19dd	149
7a710745	150	#define SCM_NUMBERP(x) (SCM_INUMP(x) \|\| SCM_NUMP(x))
d3e7e88d	151	#define SCM_NUMP(x) (!SCM_IMP(x) \
f92e85f7 MV	152	&& (((0xfcff & SCM_CELL_TYPE (x)) == scm_tc7_number) \
	153	\|\| ((0xfbff & SCM_CELL_TYPE (x)) == scm_tc7_number)))
	154	/* 0xfcff (#b1100) for 0 free, 1 big, 2 real, 3 complex, then 0xfbff (#b1011) for 4 fraction */
	155
	156	#define SCM_FRACTIONP(x) (!SCM_IMP (x) && SCM_TYP16 (x) == scm_tc16_fraction)
	157	#define SCM_SLOPPY_FRACTIONP(x) (SCM_TYP16 (x) == scm_tc16_fraction)
	158	#define SCM_FRACTION_NUMERATOR(x) ((SCM) (SCM_CELL_WORD_1 (x)))
	159	#define SCM_FRACTION_DENOMINATOR(x) ((SCM) (SCM_CELL_WORD_2 (x)))
	160	#define SCM_FRACTION_SET_NUMERATOR(x, v) ((SCM) (SCM_SET_CELL_WORD_1 ((x), (v))))
	161	#define SCM_FRACTION_SET_DENOMINATOR(x, v) ((SCM) (SCM_SET_CELL_WORD_2 ((x), (v))))
	162	/* I think the left half word is free in the type, so I'll use bit 17 */
	163	#define SCM_FRACTION_REDUCED_BIT 0x10000
	164	#define SCM_FRACTION_REDUCED_SET(x) (SCM_SET_CELL_TYPE((x), (SCM_CELL_TYPE (x) \| SCM_FRACTION_REDUCED_BIT)))
	165	#define SCM_FRACTION_REDUCED_CLEAR(x) (SCM_SET_CELL_TYPE((x), (SCM_CELL_TYPE (x) & ~SCM_FRACTION_REDUCED_BIT)))
	166	#define SCM_FRACTION_REDUCED(x) (0x10000 & SCM_CELL_TYPE (x))
88eb6852	167
0f2d19dd JB	168	\f
0f2d19dd JB	169
92c2555f	170	typedef struct scm_t_double
0f2d19dd JB	171	{
0f2d19dd JB	172	SCM type;
88eb6852 MD	173	SCM pad;
88eb6852 MD	174	double real;
92c2555f	175	} scm_t_double;
0f2d19dd	176
92c2555f	177	typedef struct scm_t_complex
0f2d19dd	178	{
88eb6852 MD	179	double real;
88eb6852 MD	180	double imag;
92c2555f	181	} scm_t_complex;
0f2d19dd	182
0f2d19dd	183	\f
0f2d19dd	184
33b001fd MV	185	SCM_API SCM scm_exact_p (SCM x);
	186	SCM_API SCM scm_odd_p (SCM n);
	187	SCM_API SCM scm_even_p (SCM n);
232898d9 MV	188	SCM_API SCM scm_inf_p (SCM n);
	189	SCM_API SCM scm_nan_p (SCM n);
	190	SCM_API SCM scm_inf (void);
	191	SCM_API SCM scm_nan (void);
33b001fd MV	192	SCM_API SCM scm_abs (SCM x);
	193	SCM_API SCM scm_quotient (SCM x, SCM y);
	194	SCM_API SCM scm_remainder (SCM x, SCM y);
	195	SCM_API SCM scm_modulo (SCM x, SCM y);
	196	SCM_API SCM scm_gcd (SCM x, SCM y);
	197	SCM_API SCM scm_lcm (SCM n1, SCM n2);
	198	SCM_API SCM scm_logand (SCM n1, SCM n2);
	199	SCM_API SCM scm_logior (SCM n1, SCM n2);
	200	SCM_API SCM scm_logxor (SCM n1, SCM n2);
	201	SCM_API SCM scm_logtest (SCM n1, SCM n2);
	202	SCM_API SCM scm_logbit_p (SCM n1, SCM n2);
	203	SCM_API SCM scm_lognot (SCM n);
d885e204	204	SCM_API SCM scm_modulo_expt (SCM n, SCM k, SCM m);
33b001fd MV	205	SCM_API SCM scm_integer_expt (SCM z1, SCM z2);
	206	SCM_API SCM scm_ash (SCM n, SCM cnt);
	207	SCM_API SCM scm_bit_extract (SCM n, SCM start, SCM end);
	208	SCM_API SCM scm_logcount (SCM n);
	209	SCM_API SCM scm_integer_length (SCM n);
1be6b49c	210
33b001fd MV	211	SCM_API size_t scm_iint2str (long num, int rad, char *p);
	212	SCM_API SCM scm_number_to_string (SCM x, SCM radix);
	213	SCM_API int scm_print_real (SCM sexp, SCM port, scm_print_state *pstate);
	214	SCM_API int scm_print_complex (SCM sexp, SCM port, scm_print_state *pstate);
	215	SCM_API int scm_bigprint (SCM exp, SCM port, scm_print_state *pstate);
	216	SCM_API SCM scm_i_mem2number (const char *mem, size_t len, unsigned int radix);
	217	SCM_API SCM scm_string_to_number (SCM str, SCM radix);
	218	SCM_API SCM scm_make_real (double x);
	219	SCM_API SCM scm_make_complex (double x, double y);
	220	SCM_API SCM scm_bigequal (SCM x, SCM y);
	221	SCM_API SCM scm_real_equalp (SCM x, SCM y);
	222	SCM_API SCM scm_complex_equalp (SCM x, SCM y);
	223	SCM_API SCM scm_number_p (SCM x);
	224	SCM_API SCM scm_real_p (SCM x);
	225	SCM_API SCM scm_integer_p (SCM x);
	226	SCM_API SCM scm_inexact_p (SCM x);
	227	SCM_API SCM scm_num_eq_p (SCM x, SCM y);
	228	SCM_API SCM scm_less_p (SCM x, SCM y);
	229	SCM_API SCM scm_gr_p (SCM x, SCM y);
	230	SCM_API SCM scm_leq_p (SCM x, SCM y);
	231	SCM_API SCM scm_geq_p (SCM x, SCM y);
	232	SCM_API SCM scm_zero_p (SCM z);
	233	SCM_API SCM scm_positive_p (SCM x);
	234	SCM_API SCM scm_negative_p (SCM x);
	235	SCM_API SCM scm_max (SCM x, SCM y);
	236	SCM_API SCM scm_min (SCM x, SCM y);
	237	SCM_API SCM scm_sum (SCM x, SCM y);
	238	SCM_API SCM scm_difference (SCM x, SCM y);
	239	SCM_API SCM scm_product (SCM x, SCM y);
	240	SCM_API double scm_num2dbl (SCM a, const char * why);
	241	SCM_API SCM scm_divide (SCM x, SCM y);
f92e85f7 MV	242	SCM_API SCM scm_floor (SCM x);
f92e85f7 MV	243	SCM_API SCM scm_ceiling (SCM x);
33b001fd MV	244	SCM_API double scm_asinh (double x);
	245	SCM_API double scm_acosh (double x);
	246	SCM_API double scm_atanh (double x);
	247	SCM_API double scm_truncate (double x);
	248	SCM_API double scm_round (double x);
f92e85f7 MV	249	SCM_API SCM scm_truncate_number (SCM x);
f92e85f7 MV	250	SCM_API SCM scm_round_number (SCM x);
33b001fd MV	251	SCM_API SCM scm_sys_expt (SCM z1, SCM z2);
	252	SCM_API SCM scm_sys_atan2 (SCM z1, SCM z2);
	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);
33b001fd MV	262
	263	SCM_API SCM scm_short2num (short n);
	264	SCM_API SCM scm_ushort2num (unsigned short n);
	265	SCM_API SCM scm_int2num (int n);
	266	SCM_API SCM scm_uint2num (unsigned int n);
	267	SCM_API SCM scm_long2num (long n);
	268	SCM_API SCM scm_ulong2num (unsigned long n);
	269	SCM_API SCM scm_size2num (size_t n);
22185d5f	270	SCM_API SCM scm_ptrdiff2num (scm_t_ptrdiff n);
33b001fd MV	271	SCM_API short scm_num2short (SCM num, unsigned long int pos,
	272	const char *s_caller);
	273	SCM_API unsigned short scm_num2ushort (SCM num, unsigned long int pos,
	274	const char *s_caller);
	275	SCM_API int scm_num2int (SCM num, unsigned long int pos,
	276	const char *s_caller);
	277	SCM_API unsigned int scm_num2uint (SCM num, unsigned long int pos,
	278	const char *s_caller);
	279	SCM_API long scm_num2long (SCM num, unsigned long int pos,
	280	const char *s_caller);
	281	SCM_API unsigned long scm_num2ulong (SCM num, unsigned long int pos,
	282	const char *s_caller);
22185d5f RB	283	SCM_API scm_t_ptrdiff scm_num2ptrdiff (SCM num, unsigned long int pos,
22185d5f RB	284	const char *s_caller);
33b001fd MV	285	SCM_API size_t scm_num2size (SCM num, unsigned long int pos,
33b001fd MV	286	const char *s_caller);
22185d5f	287	#if SCM_SIZEOF_LONG_LONG != 0
33b001fd MV	288	SCM_API SCM scm_long_long2num (long long sl);
	289	SCM_API SCM scm_ulong_long2num (unsigned long long sl);
	290	SCM_API long long scm_num2long_long (SCM num, unsigned long int pos,
	291	const char *s_caller);
	292	SCM_API unsigned long long scm_num2ulong_long (SCM num, unsigned long int pos,
	293	const char *s_caller);
9ea8cdcb	294	#endif
1be6b49c	295
33b001fd MV	296	SCM_API SCM scm_float2num (float n);
	297	SCM_API SCM scm_double2num (double n);
	298	SCM_API float scm_num2float (SCM num, unsigned long int pos,
	299	const char *s_caller);
	300	SCM_API double scm_num2double (SCM num, unsigned long int pos,
	301	const char *s_caller);
581ded70	302
d3e7e88d RB	303
	304	/* bignum internal functions */
	305	SCM_API SCM scm_i_mkbig (void);
	306	SCM_API SCM scm_i_normbig (SCM x);
	307	SCM_API int scm_i_bigcmp (SCM a, SCM b);
	308	SCM_API SCM scm_i_dbl2big (double d);
f92e85f7	309	SCM_API SCM scm_i_dbl2num (double d);
d3e7e88d RB	310	SCM_API double scm_i_big2dbl (SCM b);
	311	SCM_API SCM scm_i_short2big (short n);
	312	SCM_API SCM scm_i_ushort2big (unsigned short n);
	313	SCM_API SCM scm_i_int2big (int n);
	314	SCM_API SCM scm_i_uint2big (unsigned int n);
	315	SCM_API SCM scm_i_long2big (long n);
	316	SCM_API SCM scm_i_ulong2big (unsigned long n);
	317	SCM_API SCM scm_i_size2big (size_t n);
	318	SCM_API SCM scm_i_ptrdiff2big (scm_t_ptrdiff n);
	319
	320	#if SCM_SIZEOF_LONG_LONG != 0
	321	SCM_API SCM scm_i_long_long2big (long long n);
	322	SCM_API SCM scm_i_ulong_long2big (unsigned long long n);
	323	#endif
	324
	325
f92e85f7 MV	326	/* ratio functions */
	327	SCM_API SCM scm_make_ratio (SCM num, SCM den);
	328	SCM_API SCM scm_rationalize (SCM x, SCM err);
	329	SCM_API SCM scm_numerator (SCM z);
	330	SCM_API SCM scm_denominator (SCM z);
	331	SCM_API SCM scm_rational_p (SCM z);
	332
	333	/* fraction internal functions */
	334	SCM_API double scm_i_fraction2double (SCM z);
	335	SCM_API SCM scm_i_fraction_equalp (SCM x, SCM y);
	336	SCM_API int scm_i_print_fraction (SCM sexp, SCM port, scm_print_state *pstate);
	337
d3e7e88d	338
b4e15479 SJ	339	#ifdef GUILE_DEBUG
	340	SCM_API SCM scm_sys_check_number_conversions (void);
	341	#endif
	342
33b001fd	343	SCM_API void scm_init_numbers (void);
0f2d19dd	344
3c9a524f	345	#endif /* SCM_NUMBERS_H */
89e00824 ML	346
	347	/*
	348	Local Variables:
	349	c-file-style: "gnu"
	350	End:
	351	*/