# HG changeset patch # User Aidan Kehoe # Date 1322330354 0 # Node ID c9e5612f542425f03f595a55cb03dce1a2b7886a # Parent 3e5d5e8e4bb7fd0b9f176f35196523fc089bd975 Support the MP library on recent FreeBSD, have it pass relevant tests. src/ChangeLog addition: 2011-11-26 Aidan Kehoe * 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 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. diff -r 3e5d5e8e4bb7 -r c9e5612f5424 src/ChangeLog --- 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 + + * 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 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 * .dbxrc.in: Move from etc/.dbxrc.in. diff -r 3e5d5e8e4bb7 -r c9e5612f5424 src/number-mp.c --- 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 */ diff -r 3e5d5e8e4bb7 -r c9e5612f5424 src/number-mp.h --- 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 +END_C_DECLS #undef pow #ifdef MP_PREFIX diff -r 3e5d5e8e4bb7 -r c9e5612f5424 src/number.c --- 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); }