diff options
author | thevenyp <thevenyp@280ebfd0-de03-0410-8827-d642c229c3f4> | 2007-09-07 14:47:34 +0000 |
---|---|---|
committer | thevenyp <thevenyp@280ebfd0-de03-0410-8827-d642c229c3f4> | 2007-09-07 14:47:34 +0000 |
commit | 7166384e2317f2e513e3fdde4565ad1cde35ee65 (patch) | |
tree | e6eac49d9f2ad6f1593fe729cd8bfa22cd3d30b7 /div_d.c | |
parent | 5d5cf53719304cfde4e836364bab11a173934f83 (diff) | |
download | mpfr-7166384e2317f2e513e3fdde4565ad1cde35ee65.tar.gz |
new arithmetic functions with a double argument
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@4832 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'div_d.c')
-rw-r--r-- | div_d.c | 52 |
1 files changed, 52 insertions, 0 deletions
diff --git a/div_d.c b/div_d.c new file mode 100644 index 000000000..49ce9f85c --- /dev/null +++ b/div_d.c @@ -0,0 +1,52 @@ +/* mpfr_div_d -- divide a multiple precision floating-point number + by a machine double precision float + +Copyright 2007 Free Software Foundation, Inc. +Contributed by the Arenaire and Cacao projects, INRIA. + +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 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 Lesser General Public +License for more details. + +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., 51 Franklin St, Fifth Floor, Boston, +MA 02110-1301, USA. */ + +#include "mpfr-impl.h" + +int +mpfr_div_d (mpfr_ptr a, mpfr_srcptr b, double c, mp_rnd_t rnd_mode) +{ + int inexact; + mpfr_t d; + MPFR_SAVE_EXPO_DECL (expo); + + MPFR_LOG_FUNC (("b[%#R]=%R c%.20g rnd=%d", b, b, c, rnd_mode), + ("a[%#R]=%R", a, a)); + + MPFR_SAVE_EXPO_MARK (expo); + + mpfr_init2 (d, IEEE_DBL_MANT_DIG); + inexact = mpfr_set_d (d, c, rnd_mode); + MPFR_ASSERTN (inexact == 0); + + inexact = mpfr_div (a, b, d, rnd_mode); + if (MPFR_UNLIKELY (MPFR_IS_SINGULAR (a))) + { + MPFR_SAVE_EXPO_UPDATE_FLAGS (expo, __gmpfr_flags); + } + + + mpfr_clear(d); + MPFR_SAVE_EXPO_FREE (expo); + return mpfr_check_range (a, inexact, rnd_mode); +} |