diff options
Diffstat (limited to 'mpfr/set_q.c')
-rw-r--r-- | mpfr/set_q.c | 52 |
1 files changed, 22 insertions, 30 deletions
diff --git a/mpfr/set_q.c b/mpfr/set_q.c index c6c364d6e..886a4c09f 100644 --- a/mpfr/set_q.c +++ b/mpfr/set_q.c @@ -1,20 +1,20 @@ /* mpfr_set_q -- set a floating-point number from a multiple-precision rational -Copyright (C) 2000 Free Software Foundation. +Copyright (C) 2000, 2001 Free Software Foundation, Inc. This file is part of the MPFR Library. The MPFR Library is free software; you can redistribute it and/or modify -it under the terms of the GNU Library General Public License as published by -the Free Software Foundation; either version 2 of the License, or (at your +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. The MPFR 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 Library General Public +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 Library General Public License +You should have received a copy of the GNU Lesser General Public License along with the MPFR Library; see the file COPYING.LIB. If not, write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */ @@ -26,38 +26,30 @@ MA 02111-1307, USA. */ #include "longlong.h" /* set f to the rational q */ -void -#if __STDC__ +int mpfr_set_q (mpfr_ptr f, mpq_srcptr q, mp_rnd_t rnd) -#else -mpfr_set_q (f, q, rnd) - mpfr_ptr f; - mpq_srcptr q; - mp_rnd_t rnd; -#endif { - int sign; mpz_srcptr num, den; - unsigned int prec; - mpfr_t n,d; + mpfr_t n, d; + int inexact; MPFR_CLEAR_FLAGS(f); - num = mpq_numref(q); - sign = mpz_cmp_ui(num, 0); - if (sign==0) { - MPFR_SET_ZERO(f); - return; - } + num = mpq_numref (q); + if (mpz_cmp_ui (num, 0) == 0) + { + MPFR_SET_ZERO(f); + return 0; + } den = mpq_denref(q); - prec = MPFR_PREC(f); - mpfr_init2(n, mpz_sizeinbase(num, 2)); - mpfr_set_z(n, num, GMP_RNDZ); /* result is exact */ - mpfr_init2(d, mpz_sizeinbase(den, 2)); - mpfr_set_z(d, den, GMP_RNDZ); /* result is exact */ - MPFR_PREC(f) = prec; - mpfr_div(f, n, d, rnd); - mpfr_clear(n); mpfr_clear(d); + mpfr_init2 (n, mpz_sizeinbase(num, 2)); + mpfr_set_z (n, num, GMP_RNDZ); /* result is exact */ + mpfr_init2 (d, mpz_sizeinbase(den, 2)); + mpfr_set_z (d, den, GMP_RNDZ); /* result is exact */ + inexact = mpfr_div (f, n, d, rnd); + mpfr_clear (n); + mpfr_clear (d); + MPFR_RET(inexact); } |