Mercurial > hg > xemacs-beta
diff src/number-mp.h @ 1983:9c872f33ecbe
[xemacs-hg @ 2004-04-05 22:49:31 by james]
Add bignum, ratio, and bigfloat support.
author | james |
---|---|
date | Mon, 05 Apr 2004 22:50:11 +0000 |
parents | |
children | 2fc0e2f18322 |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/src/number-mp.h Mon Apr 05 22:50:11 2004 +0000 @@ -0,0 +1,161 @@ +/* Definitions of numeric types for XEmacs using the MP library. + Copyright (C) 2004 Jerry James. + +This file is part of XEmacs. + +XEmacs is free software; you can redistribute it and/or modify it +under the terms of the GNU General Public License as published by the +Free Software Foundation; either version 2, or (at your option) any +later version. + +XEmacs 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 General Public License +for more details. + +You should have received a copy of the GNU General Public License +along with XEmacs; see the file COPYING. If not, write to +the Free Software Foundation, Inc., 59 Temple Place - Suite 330, +Boston, MA 02111-1307, USA. */ + +/* Synched up with: Not in FSF. */ + +/* This library defines the following types: + bignum = MINT + + The MP library does not include support for ratios or bigfloats. +*/ + +#ifndef INCLUDED_number_mp_h_ +#define INCLUDED_number_mp_h_ + +/* BSD MP libraries without MP_PREFIX define a function named pow in mp.h that + has a different prototype from the one in math.h. We don't use that + function anyway, so we do this for safety purposes. However, this means + that number-mp.h must always be included before math.h. */ +#define pow mp_pow +#include <mp.h> +#undef pow + +#ifdef MP_PREFIX +#define MP_GCD mp_gcd +#define MP_ITOM mp_itom +#define MP_MADD mp_madd +#define MP_MCMP mp_mcmp +#define MP_MDIV mp_mdiv +#define MP_MFREE mp_mfree +#define MP_MSUB mp_msub +#define MP_MULT mp_mult +#define MP_SDIV mp_sdiv +#ifdef HAVE_MP_MOVE +#define MP_MOVE(x,y) mp_move (x, y) +#else +#define MP_MOVE(x,y) mp_madd (x, bignum_zero, y) +#endif +#else +#define MP_GCD gcd +#define MP_ITOM itom +#define MP_MADD madd +#define MP_MCMP mcmp +#define MP_MDIV mdiv +#define MP_MFREE mfree +#define MP_MSUB msub +#define MP_MULT mult +#define MP_SDIV sdiv +#ifdef HAVE_MP_MOVE +#define MP_MOVE(x,y) move (x, y) +#else +#define MP_MOVE(x,y) madd (x, bignum_zero, y) +#endif +#endif + +typedef MINT *bignum; + +extern void init_number_mp(void); + + +/********************************* Bignums **********************************/ + +#define HAVE_BIGNUM 1 + +extern MINT *bignum_zero, *intern_bignum; +extern MINT *bignum_min_int, *bignum_max_int, *bignum_max_uint; +extern MINT *bignum_min_long, *bignum_max_long, *bignum_max_ulong; +extern short div_rem; + +/***** Bignum: basic functions *****/ +#define bignum_init(b) (b = MP_ITOM (0)) +#define bignum_fini(b) MP_MFREE (b) +#define bignum_hashcode(b) bignum_to_uint (b) +#define bignum_sign(b) MP_MCMP (b, bignum_zero) +#define bignum_evenp(b) (MP_SDIV (b, 2, intern_bignum, &div_rem), \ + div_rem == 0) +#define bignum_oddp(b) (MP_SDIV (b, 2, intern_bignum, &div_rem), \ + div_rem != 0) + +/***** Bignum: size *****/ +#define bignum_fits_int_p(b) (MP_MCMP (b, bignum_min_int) >= 0 && \ + MP_MCMP (b, bignum_max_int) <= 0) +#define bignum_fits_uint_p(b) (MP_MCMP (b, bignum_zero) >= 0 && \ + MP_MCMP (b, bignum_max_uint) <= 0) +#define bignum_fits_long_p(b) (MP_MCMP (b, bignum_min_long) >= 0 && \ + MP_MCMP (b, bignum_max_long) <= 0) +#define bignum_fits_ulong_p(b) (MP_MCMP (b, bignum_zero) >= 0 && \ + MP_MCMP (b, bignum_max_ulong) <= 0) + +/***** Bignum: conversions *****/ +extern char *bignum_to_string(bignum, int); +extern int bignum_to_int(bignum); +extern unsigned int bignum_to_uint(bignum); +extern long bignum_to_long(bignum); +extern unsigned long bignum_to_ulong(bignum); +extern double bignum_to_double(bignum); + +/***** Bignum: converting assignments *****/ +#define bignum_set(b1, b2) MP_MOVE (b2, b1) +extern int bignum_set_string(bignum, const char *, int); +extern void bignum_set_long(bignum, long); +extern void bignum_set_ulong(bignum, unsigned long); +extern void bignum_set_double(bignum, double); + +/***** Bignum: comparisons *****/ +#define bignum_cmp(b1,b2) MP_MCMP (b1, b2) +#define bignum_lt(b1,b2) (MP_MCMP (b1, b2) < 0) +#define bignum_le(b1,b2) (MP_MCMP (b1, b2) <= 0) +#define bignum_eql(b1,b2) (MP_MCMP (b1, b2) == 0) +#define bignum_ge(b1,b2) (MP_MCMP (b1, b2) >= 0) +#define bignum_gt(b1,b2) (MP_MCMP (b1, b2) > 0) + +/***** Bignum: arithmetic *****/ +#define bignum_neg(b,b2) MP_MSUB (bignum_zero, b2, b) +#define bignum_abs(b,b2) (MP_MCMP (b2, bignum_zero) < 0 \ + ? MP_MSUB (bignum_zero, b2, b) \ + : MP_MADD (bignum_zero, b2, b)) +#define bignum_add(b,b1,b2) MP_MADD (b1, b2, b) +#define bignum_sub(b,b1,b2) MP_MSUB (b1, b2, b) +#define bignum_mul(b,b1,b2) MP_MULT (b1, b2, b) +extern int bignum_divisible_p(bignum, bignum); +#define bignum_div(b,b1,b2) MP_MDIV (b1, b2, b, intern_bignum) +extern void bignum_ceil(bignum, bignum, bignum); +extern void bignum_floor(bignum, bignum, bignum); +#define bignum_mod(b,b1,b2) MP_MDIV (b1, b2, intern_bignum, b) +extern void bignum_pow(bignum, bignum, unsigned long); +#define bignum_gcd(res,b1,b2) MP_GCD (b1, b2, res) +extern void bignum_lcm(bignum, bignum, bignum); + +/***** Bignum: bit manipulations *****/ +extern void bignum_and(bignum, bignum, bignum); +extern void bignum_ior(bignum, bignum, bignum); +extern void bignum_xor(bignum, bignum, bignum); +extern void bignum_not(bignum, bignum); +extern void bignum_setbit(bignum, unsigned long); +extern void bignum_clrbit(bignum, unsigned long); +extern int bignum_testbit(bignum, unsigned long); +extern void bignum_lshift(bignum, bignum, unsigned long); +extern void bignum_rshift(bignum, bignum, unsigned long); + +/***** Bignum: random numbers *****/ +extern void bignum_random_seed(unsigned long); +extern void bignum_random(bignum, bignum); + +#endif /* INCLUDED_number_mp_h_ */