Mercurial > hg > xemacs-beta
changeset 5602:c9e5612f5424
Support the MP library on recent FreeBSD, have it pass relevant tests.
src/ChangeLog addition:
2011-11-26 Aidan Kehoe <kehoea@parhasard.net>
* number-mp.c (bignum_to_string):
Don't overwrite the accumulator we've just set up for this
function.
* number-mp.c (BIGNUM_TO_TYPE):
mp_itom() doesn't necessarily do what this code used to think with
negative numbers, it can treat them as unsigned ints. Subtract
numbers from bignum_zero instead of multiplying them by -1 to
convert them to their negative equivalents.
* number-mp.c (bignum_to_int):
* number-mp.c (bignum_to_uint):
* number-mp.c (bignum_to_long):
* number-mp.c (bignum_to_ulong):
* number-mp.c (bignum_to_double):
Use the changed BIGNUM_TO_TYPE() in these functions.
* number-mp.c (bignum_ceil):
* number-mp.c (bignum_floor):
In these functions, be more careful about rounding to positive and
negative infinity, respectively. Don't use the sign of QUOTIENT
when working out out whether to add or subtract one, rather use
the sign QUOTIENT would have if arbitrary-precision division were
done.
* number-mp.h:
* number-mp.h (MP_GCD):
Wrap #include <mp.h> in BEGIN_C_DECLS/END_C_DECLS.
* number.c (Fbigfloat_get_precision):
* number.c (Fbigfloat_set_precision):
Don't attempt to call XBIGFLOAT_GET_PREC if this build doesn't
support big floats.
author | Aidan Kehoe <kehoea@parhasard.net> |
---|---|
date | Sat, 26 Nov 2011 17:59:14 +0000 |
parents | 3e5d5e8e4bb7 |
children | 2b2f37e84792 |
files | src/ChangeLog src/number-mp.c src/number-mp.h src/number.c |
diffstat | 4 files changed, 117 insertions(+), 23 deletions(-) [+] |
line wrap: on
line diff
--- a/src/ChangeLog Mon Nov 21 19:46:04 2011 +0100 +++ b/src/ChangeLog Sat Nov 26 17:59:14 2011 +0000 @@ -1,3 +1,39 @@ +2011-11-26 Aidan Kehoe <kehoea@parhasard.net> + + * number-mp.c (bignum_to_string): + Don't overwrite the accumulator we've just set up for this + function. + + * number-mp.c (BIGNUM_TO_TYPE): + mp_itom() doesn't necessarily do what this code used to think with + negative numbers, it can treat them as unsigned ints. Subtract + numbers from bignum_zero instead of multiplying them by -1 to + convert them to their negative equivalents. + + * number-mp.c (bignum_to_int): + * number-mp.c (bignum_to_uint): + * number-mp.c (bignum_to_long): + * number-mp.c (bignum_to_ulong): + * number-mp.c (bignum_to_double): + Use the changed BIGNUM_TO_TYPE() in these functions. + + * number-mp.c (bignum_ceil): + * number-mp.c (bignum_floor): + In these functions, be more careful about rounding to positive and + negative infinity, respectively. Don't use the sign of QUOTIENT + when working out out whether to add or subtract one, rather use + the sign QUOTIENT would have if arbitrary-precision division were + done. + + * number-mp.h: + * number-mp.h (MP_GCD): + Wrap #include <mp.h> in BEGIN_C_DECLS/END_C_DECLS. + + * number.c (Fbigfloat_get_precision): + * number.c (Fbigfloat_set_precision): + Don't attempt to call XBIGFLOAT_GET_PREC if this build doesn't + support big floats. + 2011-11-21 Marcus Crestani <crestani@informatik.uni-tuebingen.de> * .dbxrc.in: Move from etc/.dbxrc.in.
--- a/src/number-mp.c Mon Nov 21 19:46:04 2011 +0100 +++ b/src/number-mp.c Sat Nov 26 17:59:14 2011 +0000 @@ -42,7 +42,7 @@ /* FIXME: signal something if base is < 2 or doesn't fit into a short. */ /* Save the sign for later */ - sign = MP_MCMP (b, bignum_zero); + sign = bignum_sign (b); if (sign == 0) { @@ -57,8 +57,6 @@ else MP_MOVE (b, quo); - quo = MP_ITOM (0); - /* Loop over the digits of b (in BASE) and place each one into buffer */ for (i = 0U; MP_MCMP(quo, bignum_zero) > 0; i++) { @@ -90,13 +88,28 @@ } #define BIGNUM_TO_TYPE(type,accumtype) do { \ - MP_MULT (b, quo, quo); \ + if (0 == sign) \ + { \ + return (type)0; \ + } \ + \ + bignum_init (quo); \ + \ + if (sign < 0) \ + { \ + MP_MSUB (bignum_zero, b, quo); \ + } \ + else \ + { \ + MP_MOVE (b, quo); \ + } \ + \ for (i = 0U; i < sizeof(type); i++) \ { \ MP_SDIV (quo, 256, quo, &rem); \ retval |= ((accumtype) rem) << (8 * i); \ } \ - MP_MFREE (quo); \ + bignum_fini (quo); \ } while (0) int @@ -107,8 +120,7 @@ REGISTER unsigned int i; MINT *quo; - sign = MP_MCMP (b, bignum_zero) < 0 ? -1 : 1; - quo = MP_ITOM (sign); + sign = bignum_sign (b); BIGNUM_TO_TYPE (int, unsigned int); return ((int) retval) * sign; } @@ -116,12 +128,12 @@ unsigned int bignum_to_uint (bignum b) { - short rem; + short rem, sign; unsigned int retval = 0U; REGISTER unsigned int i; MINT *quo; - quo = MP_ITOM (MP_MCMP (b, bignum_zero) < 0 ? -1 : 1); + sign = bignum_sign (b); BIGNUM_TO_TYPE (unsigned int, unsigned int); return retval; } @@ -134,8 +146,7 @@ REGISTER unsigned int i; MINT *quo; - sign = MP_MCMP (b, bignum_zero) < 0 ? -1 : 1; - quo = MP_ITOM (sign); + sign = bignum_sign (b); BIGNUM_TO_TYPE (long, unsigned long); return ((long) retval) * sign; } @@ -143,12 +154,12 @@ unsigned long bignum_to_ulong (bignum b) { - short rem; + short rem, sign; unsigned long retval = 0UL; REGISTER unsigned int i; MINT *quo; - quo = MP_ITOM (MP_MCMP (b, bignum_zero) < 0 ? -1 : 1); + sign = bignum_sign (b); BIGNUM_TO_TYPE (unsigned long, unsigned long); return retval; } @@ -161,9 +172,17 @@ REGISTER unsigned int i; MINT *quo; - sign = MP_MCMP (b, bignum_zero) < 0 ? -1 : 1; - quo = MP_ITOM (sign); - MP_MULT (b, quo, quo); + sign = bignum_sign (b); + bignum_init (quo); + if (sign < 0) + { + MP_MSUB (bignum_zero, b, quo); + } + else + { + MP_MOVE (b, quo); + } + for (i = 0U; MP_MCMP (quo, bignum_zero) > 0; i++) { MP_SDIV (quo, 256, quo, &rem); @@ -303,17 +322,51 @@ void bignum_ceil (bignum quotient, bignum N, bignum D) { MP_MDIV (N, D, quotient, intern_bignum); - if (MP_MCMP (intern_bignum, bignum_zero) > 0 && - MP_MCMP (quotient, bignum_zero) > 0) - MP_MADD (quotient, bignum_one, quotient); + MP_MDIV (N, D, quotient, intern_bignum); + if (MP_MCMP (intern_bignum, bignum_zero) != 0) + { + short signN = MP_MCMP (N, bignum_zero); + short signD = MP_MCMP (D, bignum_zero); + + /* If the quotient is positive, add one, since we're rounding to + positive infinity. */ + if (signD < 0) + { + if (signN <= 0) + { + MP_MADD (quotient, bignum_one, quotient); + } + } + else if (signN >= 0) + { + MP_MADD (quotient, bignum_one, quotient); + } + } } void bignum_floor (bignum quotient, bignum N, bignum D) { MP_MDIV (N, D, quotient, intern_bignum); - if (MP_MCMP (intern_bignum, bignum_zero) > 0 && - MP_MCMP (quotient, bignum_zero) < 0) - MP_MSUB (quotient, bignum_one, quotient); + + if (MP_MCMP (intern_bignum, bignum_zero) != 0) + { + short signN = MP_MCMP (N, bignum_zero); + short signD = MP_MCMP (D, bignum_zero); + + /* If the quotient is negative, subtract one, we're rounding to minus + infinity. */ + if (signD < 0) + { + if (signN >= 0) + { + MP_MSUB (quotient, bignum_one, quotient); + } + } + else if (signN < 0) + { + MP_MSUB (quotient, bignum_one, quotient); + } + } } /* RESULT = N to the POWth power */
--- a/src/number-mp.h Mon Nov 21 19:46:04 2011 +0100 +++ b/src/number-mp.h Sat Nov 26 17:59:14 2011 +0000 @@ -32,7 +32,9 @@ 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 +BEGIN_C_DECLS #include <mp.h> +END_C_DECLS #undef pow #ifdef MP_PREFIX
--- a/src/number.c Mon Nov 21 19:46:04 2011 +0100 +++ b/src/number.c Sat Nov 26 17:59:14 2011 +0000 @@ -328,12 +328,14 @@ (f)) { CHECK_BIGFLOAT (f); +#ifdef HAVE_BIGFLOAT #ifdef HAVE_BIGNUM bignum_set_ulong (scratch_bignum, XBIGFLOAT_GET_PREC (f)); return Fcanonicalize_number (make_bignum_bg (scratch_bignum)); #else return make_fixnum ((int) XBIGFLOAT_GET_PREC (f)); #endif +#endif } DEFUN ("bigfloat-set-precision", Fbigfloat_set_precision, 2, 2, 0, /* @@ -364,8 +366,9 @@ dead_wrong_type_argument (Qintegerp, f); return Qnil; } - +#ifdef HAVE_BIGFLOAT XBIGFLOAT_SET_PREC (f, prec); +#endif return Fbigfloat_get_precision (f); }