diff options
Diffstat (limited to 'src/floatfns.c')
| -rw-r--r-- | src/floatfns.c | 44 |
1 files changed, 13 insertions, 31 deletions
diff --git a/src/floatfns.c b/src/floatfns.c index 77e20d5640b..2f33b8652b2 100644 --- a/src/floatfns.c +++ b/src/floatfns.c | |||
| @@ -270,11 +270,8 @@ DEFUN ("abs", Fabs, Sabs, 1, 1, 0, | |||
| 270 | { | 270 | { |
| 271 | if (mpz_sgn (XBIGNUM (arg)->value) < 0) | 271 | if (mpz_sgn (XBIGNUM (arg)->value) < 0) |
| 272 | { | 272 | { |
| 273 | mpz_t val; | 273 | mpz_neg (mpz[0], XBIGNUM (arg)->value); |
| 274 | mpz_init (val); | 274 | arg = make_integer_mpz (); |
| 275 | mpz_neg (val, XBIGNUM (arg)->value); | ||
| 276 | arg = make_integer (val); | ||
| 277 | mpz_clear (val); | ||
| 278 | } | 275 | } |
| 279 | } | 276 | } |
| 280 | 277 | ||
| @@ -360,20 +357,10 @@ rounding_driver (Lisp_Object arg, Lisp_Object divisor, | |||
| 360 | { | 357 | { |
| 361 | if (EQ (divisor, make_fixnum (0))) | 358 | if (EQ (divisor, make_fixnum (0))) |
| 362 | xsignal0 (Qarith_error); | 359 | xsignal0 (Qarith_error); |
| 363 | mpz_t d, q; | 360 | int_divide (mpz[0], |
| 364 | mpz_init (d); | 361 | *bignum_integer (&mpz[0], arg), |
| 365 | mpz_init (q); | 362 | *bignum_integer (&mpz[1], divisor)); |
| 366 | int_divide (q, | 363 | return make_integer_mpz (); |
| 367 | (FIXNUMP (arg) | ||
| 368 | ? (mpz_set_intmax (q, XFIXNUM (arg)), q) | ||
| 369 | : XBIGNUM (arg)->value), | ||
| 370 | (FIXNUMP (divisor) | ||
| 371 | ? (mpz_set_intmax (d, XFIXNUM (divisor)), d) | ||
| 372 | : XBIGNUM (divisor)->value)); | ||
| 373 | Lisp_Object result = make_integer (q); | ||
| 374 | mpz_clear (d); | ||
| 375 | mpz_clear (q); | ||
| 376 | return result; | ||
| 377 | } | 364 | } |
| 378 | 365 | ||
| 379 | double f1 = FLOATP (arg) ? XFLOAT_DATA (arg) : XFIXNUM (arg); | 366 | double f1 = FLOATP (arg) ? XFLOAT_DATA (arg) : XFIXNUM (arg); |
| @@ -417,20 +404,15 @@ rounddiv_q (mpz_t q, mpz_t const n, mpz_t const d) | |||
| 417 | if (abs_r1 < abs_r + (q & 1)) | 404 | if (abs_r1 < abs_r + (q & 1)) |
| 418 | q += neg_d == neg_r ? 1 : -1; */ | 405 | q += neg_d == neg_r ? 1 : -1; */ |
| 419 | 406 | ||
| 420 | mpz_t r, abs_r1; | 407 | mpz_t *r = &mpz[2], *abs_r = r, *abs_r1 = &mpz[3]; |
| 421 | mpz_init (r); | 408 | mpz_tdiv_qr (q, *r, n, d); |
| 422 | mpz_init (abs_r1); | ||
| 423 | mpz_tdiv_qr (q, r, n, d); | ||
| 424 | bool neg_d = mpz_sgn (d) < 0; | 409 | bool neg_d = mpz_sgn (d) < 0; |
| 425 | bool neg_r = mpz_sgn (r) < 0; | 410 | bool neg_r = mpz_sgn (*r) < 0; |
| 426 | mpz_t *abs_r = &r; | 411 | mpz_abs (*abs_r, *r); |
| 427 | mpz_abs (*abs_r, r); | 412 | mpz_abs (*abs_r1, d); |
| 428 | mpz_abs (abs_r1, d); | 413 | mpz_sub (*abs_r1, *abs_r1, *abs_r); |
| 429 | mpz_sub (abs_r1, abs_r1, *abs_r); | 414 | if (mpz_cmp (*abs_r1, *abs_r) < (mpz_odd_p (q) != 0)) |
| 430 | if (mpz_cmp (abs_r1, *abs_r) < (mpz_odd_p (q) != 0)) | ||
| 431 | (neg_d == neg_r ? mpz_add_ui : mpz_sub_ui) (q, q, 1); | 415 | (neg_d == neg_r ? mpz_add_ui : mpz_sub_ui) (q, q, 1); |
| 432 | mpz_clear (r); | ||
| 433 | mpz_clear (abs_r1); | ||
| 434 | } | 416 | } |
| 435 | 417 | ||
| 436 | /* The code uses emacs_rint, so that it works to undefine HAVE_RINT | 418 | /* The code uses emacs_rint, so that it works to undefine HAVE_RINT |