# ChangeLog for dev-libs/mpfr
# Copyright 1999-2007 Gentoo Foundation; Distributed under the GPL v2
-# $Header: /var/cvsroot/gentoo-x86/dev-libs/mpfr/ChangeLog,v 1.46 2007/02/17 11:52:32 eroyf Exp $
+# $Header: /var/cvsroot/gentoo-x86/dev-libs/mpfr/ChangeLog,v 1.47 2007/04/06 14:07:46 dragonheart Exp $
+
+*mpfr-2.2.1_p5 (06 Apr 2007)
+
+ 06 Apr 2007; Daniel Black <dragonheart@gentoo.org> +files/2.2.1/patch02,
+ +files/2.2.1/patch03, +files/2.2.1/patch04, +files/2.2.1/patch05,
+ +mpfr-2.2.1_p5.ebuild:
+ upstream patches
17 Feb 2007; Alexander H. Færøy <eroyf@gentoo.org>
mpfr-2.2.0_p16.ebuild:
MD5 fd276c3936d6fa5adf9fae1cd57c74c4 files/2.2.1/patch01 1826
RMD160 8e5d21a2b43062a31d13bca242f334d42b552bee files/2.2.1/patch01 1826
SHA256 0b19e4b6484e505aa9e518b8942df3dee686d17d17cdbc2c41f7d444f9afd8f8 files/2.2.1/patch01 1826
+AUX 2.2.1/patch02 1687 RMD160 ef09cb3386c4312952e94d5a9717c29d75dcc92d SHA1 e93c31a2baaef5367a710a5e64b82744f5de4b64 SHA256 eb100b8c4990c047bf61fc53ed1006fa1087277eebe484e91827c018e3ad4768
+MD5 03e638937ec7a2af4f703577396e8e52 files/2.2.1/patch02 1687
+RMD160 ef09cb3386c4312952e94d5a9717c29d75dcc92d files/2.2.1/patch02 1687
+SHA256 eb100b8c4990c047bf61fc53ed1006fa1087277eebe484e91827c018e3ad4768 files/2.2.1/patch02 1687
+AUX 2.2.1/patch03 11655 RMD160 866cc9cc1d01e21027dae61cffc585cb9dd12e66 SHA1 4eaa326c4fe43499a8acbdd6fc71c6a408663650 SHA256 2d2dc48ab57184fdb502d02faf1070d0be574882107fe4cd3fd1644688a0e68e
+MD5 fa3398730e7d23fdf37ee4e5c422bf80 files/2.2.1/patch03 11655
+RMD160 866cc9cc1d01e21027dae61cffc585cb9dd12e66 files/2.2.1/patch03 11655
+SHA256 2d2dc48ab57184fdb502d02faf1070d0be574882107fe4cd3fd1644688a0e68e files/2.2.1/patch03 11655
+AUX 2.2.1/patch04 7509 RMD160 3d4c5ab82824f66259841a8c059ce9625a7a4891 SHA1 9a1fec3686cfd6b227cebcb9e24d4ed8549ad5d9 SHA256 f449733f1a72df18421d8f471118fc2c77e408654361bf5cab7cce41c8aaf26a
+MD5 5d44d16448f85fea2acd5bcbd04cabb0 files/2.2.1/patch04 7509
+RMD160 3d4c5ab82824f66259841a8c059ce9625a7a4891 files/2.2.1/patch04 7509
+SHA256 f449733f1a72df18421d8f471118fc2c77e408654361bf5cab7cce41c8aaf26a files/2.2.1/patch04 7509
+AUX 2.2.1/patch05 16992 RMD160 38a4558d6890ca12b613b5b482046b19d1cf3b6b SHA1 4d559340e1006e84dae50e655723e916a26989ad SHA256 304b03f5e99e4d1790107b4e7d0b1d54305b5d06203fc1723758f59d3d36947b
+MD5 2881b39c442466caa3451725b50f5a7a files/2.2.1/patch05 16992
+RMD160 38a4558d6890ca12b613b5b482046b19d1cf3b6b files/2.2.1/patch05 16992
+SHA256 304b03f5e99e4d1790107b4e7d0b1d54305b5d06203fc1723758f59d3d36947b files/2.2.1/patch05 16992
DIST mpfr-2.2.0.tar.bz2 728490 RMD160 5243eb0387392123d9d30ca009b30147a51370a6 SHA1 3894bf68982e6218fd94ee743f05c0554a35536a SHA256 cea002a11707b5f8fb2147686145563e5f313fc62558e8400444bdcf89971e6a
DIST mpfr-2.2.0_p5 27251 RMD160 7e60b2be4691235de88646534b64cd12cdc1a45d SHA1 b2600dead06a074fc4485d1a2f3e4940a1c54134 SHA256 1e253c48800667d295e5ec3407b8e27a86e385706ecb49a7dc51503839752309
DIST mpfr-2.2.1.tar.bz2 787634 RMD160 cb4c24dc356412e3bae30242f426a2a02d146d04 SHA1 d0abd0943fcfc937a816bf67baaff3f0273a19e7 SHA256 92e2fea0eefd3d50b8e762334326657bfc5be9a6bda58cfa724b2a5cba9f65e4
MD5 b3a03005931420dbb0f7a4255922f0f3 mpfr-2.2.1_p1.ebuild 1425
RMD160 8dd8637bd7619c67ea15df13d45e1ec66288d0a0 mpfr-2.2.1_p1.ebuild 1425
SHA256 d1d396760cc0332c0f1a109ebf09d30be4564db22187c350577582140a7901d3 mpfr-2.2.1_p1.ebuild 1425
-MISC ChangeLog 5973 RMD160 863242400fd4e741cc45c001f15cae78fd5ba2fa SHA1 51bc4d05e6ab10e22044e12321475be28e6725b7 SHA256 299dc60a6a03662100aca87b4eb8618ef081923fb148e2b71e0cf2f2b778ce1b
-MD5 a3137909a9b89b57ce11dee14bf6868f ChangeLog 5973
-RMD160 863242400fd4e741cc45c001f15cae78fd5ba2fa ChangeLog 5973
-SHA256 299dc60a6a03662100aca87b4eb8618ef081923fb148e2b71e0cf2f2b778ce1b ChangeLog 5973
+EBUILD mpfr-2.2.1_p5.ebuild 1425 RMD160 b608401989c135e68bf2e922a1cf59ac9e9e8a00 SHA1 8c278029ccd45f2fde826c18bb027235d7e1733d SHA256 dbaca0552d6d20b55db957b79ee21465df64313e252c2236670da9d8981f2055
+MD5 47de7220593c3a4b80b8f26c9e5383f5 mpfr-2.2.1_p5.ebuild 1425
+RMD160 b608401989c135e68bf2e922a1cf59ac9e9e8a00 mpfr-2.2.1_p5.ebuild 1425
+SHA256 dbaca0552d6d20b55db957b79ee21465df64313e252c2236670da9d8981f2055 mpfr-2.2.1_p5.ebuild 1425
+MISC ChangeLog 6197 RMD160 e4a9eafde33ee48e6be09b24e8774572dfaf0906 SHA1 44b219acdf8e8316f8202f6402416dc9923fb426 SHA256 d22c2c7af8402d6bd5326b526c289a698248e2dc3d2a26cc9f752a1db91e377d
+MD5 4878f21941ec5ce4fca9802803719389 ChangeLog 6197
+RMD160 e4a9eafde33ee48e6be09b24e8774572dfaf0906 ChangeLog 6197
+SHA256 d22c2c7af8402d6bd5326b526c289a698248e2dc3d2a26cc9f752a1db91e377d ChangeLog 6197
MISC metadata.xml 162 RMD160 d002486a43522f2116b1d9d59828c484956d66e2 SHA1 d6b4923897f6ae673b4f93646f5b4ba61d5a2c3c SHA256 65a915d44de1f01d4b7f72d313b4192c38374a9835d24988c00c1e73dca5805a
MD5 567094e03359ffc1c95af7356395228d metadata.xml 162
RMD160 d002486a43522f2116b1d9d59828c484956d66e2 metadata.xml 162
MD5 2e4c88b3a3810fc13043d40304378bbd files/digest-mpfr-2.2.1_p1 235
RMD160 dcfe570bc7a569fdce138a0a1e8e4d4776463f3a files/digest-mpfr-2.2.1_p1 235
SHA256 138a6c0286b4a439b3c422770c9b41093baa2d0822a06074ce71f8aef4c3a927 files/digest-mpfr-2.2.1_p1 235
+MD5 2e4c88b3a3810fc13043d40304378bbd files/digest-mpfr-2.2.1_p5 235
+RMD160 dcfe570bc7a569fdce138a0a1e8e4d4776463f3a files/digest-mpfr-2.2.1_p5 235
+SHA256 138a6c0286b4a439b3c422770c9b41093baa2d0822a06074ce71f8aef4c3a927 files/digest-mpfr-2.2.1_p5 235
-----BEGIN PGP SIGNATURE-----
-Version: GnuPG v1.4.6 (GNU/Linux)
+Version: GnuPG v2.0.2 (GNU/Linux)
-iD8DBQFF1u0Y6Xxn2nLwtdQRAngNAJ9TS1gfxXgAZ6JdQMi+ey9tOCRdZwCgvRyL
-lEPS5ruT28iv3xwF6v1MBYE=
-=kh5P
+iD8DBQFGFlRQmdTrptrqvGERAj18AJ4vV8Xyp58V3QdoOOLv6BXUk3KPaQCeJdzs
+7Ol3tRdZInf+N3YKnIitOoQ=
+=zT87
-----END PGP SIGNATURE-----
--- /dev/null
+diff -Naurd mpfr-2.2.1-p1/VERSION mpfr-2.2.1-p2/VERSION
+--- mpfr-2.2.1-p1/VERSION 2007-02-02 00:01:41.000000000 +0000
++++ mpfr-2.2.1-p2/VERSION 2007-02-15 23:37:58.000000000 +0000
+@@ -1 +1 @@
+-2.2.1-p1
++2.2.1-p2
+diff -Naurd mpfr-2.2.1-p1/mpfr.h mpfr-2.2.1-p2/mpfr.h
+--- mpfr-2.2.1-p1/mpfr.h 2007-02-02 00:01:41.000000000 +0000
++++ mpfr-2.2.1-p2/mpfr.h 2007-02-15 23:37:58.000000000 +0000
+@@ -26,7 +26,7 @@
+ #define MPFR_VERSION_MAJOR 2
+ #define MPFR_VERSION_MINOR 2
+ #define MPFR_VERSION_PATCHLEVEL 1
+-#define MPFR_VERSION_STRING "2.2.1-p1"
++#define MPFR_VERSION_STRING "2.2.1-p2"
+
+ /* Macros dealing with MPFR VERSION */
+ #define MPFR_VERSION_NUM(a,b,c) (((a) << 16L) | ((b) << 8) | (c))
+diff -Naurd mpfr-2.2.1-p1/tests/tget_f.c mpfr-2.2.1-p2/tests/tget_f.c
+--- mpfr-2.2.1-p1/tests/tget_f.c 2006-11-29 09:49:47.000000000 +0000
++++ mpfr-2.2.1-p2/tests/tget_f.c 2007-02-15 23:37:53.000000000 +0000
+@@ -153,8 +153,9 @@
+ }
+ if (i <= - (unsigned long) LONG_MIN)
+ {
+- mpfr_set_si (y, - (long) i, GMP_RNDN);
+- if (mpfr_get_f (x, y, GMP_RNDN) || mpf_cmp_si (x, - (long) i))
++ long j = i < - (unsigned long) LONG_MIN ? - (long) i : LONG_MIN;
++ mpfr_set_si (y, j, GMP_RNDN);
++ if (mpfr_get_f (x, y, GMP_RNDN) || mpf_cmp_si (x, j))
+ {
+ printf ("Error: mpfr_get_f(-%lu) fails\n", i);
+ exit (1);
+diff -Naurd mpfr-2.2.1-p1/version.c mpfr-2.2.1-p2/version.c
+--- mpfr-2.2.1-p1/version.c 2007-02-02 00:01:41.000000000 +0000
++++ mpfr-2.2.1-p2/version.c 2007-02-15 23:37:58.000000000 +0000
+@@ -24,5 +24,5 @@
+ const char *
+ mpfr_get_version (void)
+ {
+- return "2.2.1-p1";
++ return "2.2.1-p2";
+ }
--- /dev/null
+diff -Naurd mpfr-2.2.1-p2/VERSION mpfr-2.2.1-p3/VERSION
+--- mpfr-2.2.1-p2/VERSION 2007-02-15 23:37:58.000000000 +0000
++++ mpfr-2.2.1-p3/VERSION 2007-02-15 23:40:49.000000000 +0000
+@@ -1 +1 @@
+-2.2.1-p2
++2.2.1-p3
+diff -Naurd mpfr-2.2.1-p2/mpfr.h mpfr-2.2.1-p3/mpfr.h
+--- mpfr-2.2.1-p2/mpfr.h 2007-02-15 23:37:58.000000000 +0000
++++ mpfr-2.2.1-p3/mpfr.h 2007-02-15 23:40:49.000000000 +0000
+@@ -26,7 +26,7 @@
+ #define MPFR_VERSION_MAJOR 2
+ #define MPFR_VERSION_MINOR 2
+ #define MPFR_VERSION_PATCHLEVEL 1
+-#define MPFR_VERSION_STRING "2.2.1-p2"
++#define MPFR_VERSION_STRING "2.2.1-p3"
+
+ /* Macros dealing with MPFR VERSION */
+ #define MPFR_VERSION_NUM(a,b,c) (((a) << 16L) | ((b) << 8) | (c))
+@@ -626,13 +626,13 @@
+ _p->_mpfr_sign = 1; \
+ _p->_mpfr_exp = __MPFR_EXP_ZERO; \
+ (void) (_r); 0; }) : \
+- mpfr_set_ui (_f,_u,_r))
++ mpfr_set_ui_2exp ((_f), (_u), 0, (_r)))
+ #endif
+ #undef mpfr_set_si
+ #define mpfr_set_si(_f,_s,_r) \
+ (__builtin_constant_p (_s) && (_s) >= 0 ? \
+ mpfr_set_ui ((_f), (_s), (_r)) : \
+- mpfr_set_si ((_f), (_s), (_r)))
++ mpfr_set_si_2exp ((_f), (_s), 0, (_r)))
+ #endif
+ #endif
+
+diff -Naurd mpfr-2.2.1-p2/set_si.c mpfr-2.2.1-p3/set_si.c
+--- mpfr-2.2.1-p2/set_si.c 2006-11-29 09:49:47.000000000 +0000
++++ mpfr-2.2.1-p3/set_si.c 2007-02-15 23:40:22.000000000 +0000
+@@ -19,62 +19,11 @@
+ the Free Software Foundation, Inc., 51 Franklin St, Fifth Floor, Boston,
+ MA 02110-1301, USA. */
+
+-#define MPFR_NEED_LONGLONG_H
+ #include "mpfr-impl.h"
+
++#undef mpfr_set_si
+ int
+-(mpfr_set_si) (mpfr_ptr x, long i, mp_rnd_t rnd_mode)
++mpfr_set_si (mpfr_ptr x, long i, mp_rnd_t rnd_mode)
+ {
+- int inex;
+- mp_size_t xn;
+- unsigned int cnt, nbits;
+- mp_limb_t ai, *xp;
+-
+- MPFR_CLEAR_FLAGS(x);
+- if (i == 0)
+- {
+- MPFR_SET_ZERO(x);
+- MPFR_SET_POS(x);
+- MPFR_RET(0);
+- }
+-
+- xn = (MPFR_PREC(x)-1)/BITS_PER_MP_LIMB;
+- ai = SAFE_ABS(unsigned long, i);
+- MPFR_ASSERTN(SAFE_ABS(unsigned long, i) == ai);
+- count_leading_zeros(cnt, ai);
+-
+- xp = MPFR_MANT(x);
+- xp[xn] = ai << cnt;
+- /* don't forget to put zero in lower limbs */
+- MPN_ZERO(xp, xn);
+- /* set sign */
+- if (i < 0)
+- MPFR_SET_NEG(x);
+- else
+- MPFR_SET_POS(x);
+-
+- nbits = BITS_PER_MP_LIMB - cnt;
+- MPFR_EXP (x) = nbits; /* may be out-of-range, check range below */
+- inex = mpfr_check_range(x, 0, rnd_mode);
+- if (inex)
+- return inex; /* underflow or overflow */
+-
+- /* round if MPFR_PREC(x) smaller than length of i */
+- if (MPFR_UNLIKELY(MPFR_PREC(x) < nbits))
+- {
+- int carry;
+-
+- carry = mpfr_round_raw(xp+xn, xp+xn, nbits, (i < 0), MPFR_PREC(x),
+- rnd_mode, &inex);
+- if (MPFR_UNLIKELY(carry))
+- {
+- /* nbits is the current exponent */
+- if (MPFR_UNLIKELY((mp_exp_t) nbits == __gmpfr_emax))
+- return mpfr_overflow(x, rnd_mode, (i < 0 ? -1 : 1));
+- MPFR_SET_EXP (x, nbits + 1);
+- xp[xn] = MPFR_LIMB_HIGHBIT;
+- }
+- }
+-
+- MPFR_RET(inex);
++ return mpfr_set_si_2exp (x, i, 0, rnd_mode);
+ }
+diff -Naurd mpfr-2.2.1-p2/set_si_2exp.c mpfr-2.2.1-p3/set_si_2exp.c
+--- mpfr-2.2.1-p2/set_si_2exp.c 2006-11-29 09:49:47.000000000 +0000
++++ mpfr-2.2.1-p3/set_si_2exp.c 2007-02-15 23:40:22.000000000 +0000
+@@ -20,18 +20,54 @@
+ the Free Software Foundation, Inc., 51 Franklin St, Fifth Floor, Boston,
+ MA 02110-1301, USA. */
+
++#define MPFR_NEED_LONGLONG_H
+ #include "mpfr-impl.h"
+
+ int
+ mpfr_set_si_2exp (mpfr_ptr x, long i, mp_exp_t e, mp_rnd_t rnd_mode)
+ {
+- int res;
+- MPFR_SAVE_EXPO_DECL (expo);
++ if (i == 0)
++ {
++ MPFR_SET_ZERO (x);
++ MPFR_SET_POS (x);
++ MPFR_RET (0);
++ }
++ else
++ {
++ mp_size_t xn;
++ unsigned int cnt, nbits;
++ mp_limb_t ai, *xp;
++ int inex = 0;
+
+- MPFR_SAVE_EXPO_MARK (expo);
+- res = mpfr_set_si (x, i, rnd_mode);
+- mpfr_mul_2si (x, x, e, rnd_mode); /* Should be exact */
+- MPFR_SAVE_EXPO_FREE (expo);
+- res = mpfr_check_range(x, res, rnd_mode);
+- return res;
++ /* FIXME: support int limbs (e.g. 16-bit limbs on 16-bit proc) */
++ ai = SAFE_ABS (unsigned long, i);
++ MPFR_ASSERTN (SAFE_ABS (unsigned long, i) == ai);
++
++ /* Position of the highest limb */
++ xn = (MPFR_PREC (x) - 1) / BITS_PER_MP_LIMB;
++ count_leading_zeros (cnt, ai);
++ MPFR_ASSERTD (cnt < BITS_PER_MP_LIMB); /* OK since i != 0 */
++
++ xp = MPFR_MANT(x);
++ xp[xn] = ai << cnt;
++ /* Zero the xn lower limbs. */
++ MPN_ZERO(xp, xn);
++ MPFR_SET_SIGN (x, i < 0 ? MPFR_SIGN_NEG : MPFR_SIGN_POS);
++
++ nbits = BITS_PER_MP_LIMB - cnt;
++ e += nbits; /* exponent _before_ the rounding */
++
++ /* round if MPFR_PREC(x) smaller than length of i */
++ if (MPFR_UNLIKELY (MPFR_PREC (x) < nbits) &&
++ MPFR_UNLIKELY (mpfr_round_raw (xp + xn, xp + xn, nbits, i < 0,
++ MPFR_PREC (x), rnd_mode, &inex)))
++ {
++ e++;
++ xp[xn] = MPFR_LIMB_HIGHBIT;
++ }
++
++ MPFR_CLEAR_FLAGS (x);
++ MPFR_EXP (x) = e;
++ return mpfr_check_range (x, inex, rnd_mode);
++ }
+ }
+diff -Naurd mpfr-2.2.1-p2/set_ui.c mpfr-2.2.1-p3/set_ui.c
+--- mpfr-2.2.1-p2/set_ui.c 2006-11-29 09:49:47.000000000 +0000
++++ mpfr-2.2.1-p3/set_ui.c 2007-02-15 23:40:21.000000000 +0000
+@@ -19,60 +19,11 @@
+ the Free Software Foundation, Inc., 51 Franklin St, Fifth Floor, Boston,
+ MA 02110-1301, USA. */
+
+-#define MPFR_NEED_LONGLONG_H
+ #include "mpfr-impl.h"
+
+ #undef mpfr_set_ui
+ int
+ mpfr_set_ui (mpfr_ptr x, unsigned long i, mp_rnd_t rnd_mode)
+ {
+- int inex;
+-
+- MPFR_SET_POS(x);
+- if (i == 0)
+- {
+- MPFR_SET_ZERO(x);
+- MPFR_RET(0);
+- }
+- else
+- {
+- mp_size_t xn;
+- unsigned int cnt, nbits;
+- mp_limb_t *xp;
+-
+- MPFR_CLEAR_FLAGS(x);
+-
+- xn = (MPFR_PREC(x)-1)/BITS_PER_MP_LIMB;
+- MPFR_ASSERTD(i == (mp_limb_t) i);
+- count_leading_zeros(cnt, (mp_limb_t) i);
+-
+- xp = MPFR_MANT(x);
+- xp[xn] = ((mp_limb_t) i) << cnt;
+- /* don't forget to put zero in lower limbs */
+- MPN_ZERO(xp, xn);
+-
+- nbits = BITS_PER_MP_LIMB - cnt;
+- MPFR_EXP (x) = nbits; /* may be out-of-range, check range below */
+- inex = mpfr_check_range(x, 0, rnd_mode);
+- if (inex)
+- return inex; /* underflow or overflow */
+-
+- /* round if MPFR_PREC(x) smaller than length of i */
+- if (MPFR_UNLIKELY( MPFR_PREC(x) < nbits) )
+- {
+- int carry;
+- carry = mpfr_round_raw(xp+xn, xp+xn, nbits, 0, MPFR_PREC(x),
+- rnd_mode, &inex);
+- if (MPFR_UNLIKELY(carry))
+- {
+- /* nbits is the current exponent */
+- if (MPFR_UNLIKELY((mp_exp_t) nbits == __gmpfr_emax))
+- return mpfr_overflow(x, rnd_mode, 1);
+-
+- MPFR_SET_EXP (x, nbits + 1);
+- xp[xn] = MPFR_LIMB_HIGHBIT;
+- }
+- }
+- }
+- MPFR_RET(inex);
++ return mpfr_set_ui_2exp (x, i, 0, rnd_mode);
+ }
+diff -Naurd mpfr-2.2.1-p2/set_ui_2exp.c mpfr-2.2.1-p3/set_ui_2exp.c
+--- mpfr-2.2.1-p2/set_ui_2exp.c 2006-11-29 09:49:47.000000000 +0000
++++ mpfr-2.2.1-p3/set_ui_2exp.c 2007-02-15 23:40:22.000000000 +0000
+@@ -20,18 +20,53 @@
+ the Free Software Foundation, Inc., 51 Franklin St, Fifth Floor, Boston,
+ MA 02110-1301, USA. */
+
++#define MPFR_NEED_LONGLONG_H
+ #include "mpfr-impl.h"
+
+ int
+ mpfr_set_ui_2exp (mpfr_ptr x, unsigned long i, mp_exp_t e, mp_rnd_t rnd_mode)
+ {
+- int res;
+- MPFR_SAVE_EXPO_DECL (expo);
++ MPFR_SET_POS (x);
+
+- MPFR_SAVE_EXPO_MARK (expo);
+- res = mpfr_set_ui (x, i, rnd_mode);
+- mpfr_mul_2si (x, x, e, rnd_mode); /* Should be exact */
+- MPFR_SAVE_EXPO_FREE (expo);
+- res = mpfr_check_range(x, res, rnd_mode);
+- return res;
++ if (i == 0)
++ {
++ MPFR_SET_ZERO (x);
++ MPFR_RET (0);
++ }
++ else
++ {
++ mp_size_t xn;
++ unsigned int cnt, nbits;
++ mp_limb_t *xp;
++ int inex = 0;
++
++ /* FIXME: support int limbs (e.g. 16-bit limbs on 16-bit proc) */
++ MPFR_ASSERTD (i == (mp_limb_t) i);
++
++ /* Position of the highest limb */
++ xn = (MPFR_PREC (x) - 1) / BITS_PER_MP_LIMB;
++ count_leading_zeros (cnt, (mp_limb_t) i);
++ MPFR_ASSERTD (cnt < BITS_PER_MP_LIMB); /* OK since i != 0 */
++
++ xp = MPFR_MANT(x);
++ xp[xn] = ((mp_limb_t) i) << cnt;
++ /* Zero the xn lower limbs. */
++ MPN_ZERO(xp, xn);
++
++ nbits = BITS_PER_MP_LIMB - cnt;
++ e += nbits; /* exponent _before_ the rounding */
++
++ /* round if MPFR_PREC(x) smaller than length of i */
++ if (MPFR_UNLIKELY (MPFR_PREC (x) < nbits) &&
++ MPFR_UNLIKELY (mpfr_round_raw (xp + xn, xp + xn, nbits, 0,
++ MPFR_PREC (x), rnd_mode, &inex)))
++ {
++ e++;
++ xp[xn] = MPFR_LIMB_HIGHBIT;
++ }
++
++ MPFR_CLEAR_FLAGS (x);
++ MPFR_EXP (x) = e;
++ return mpfr_check_range (x, inex, rnd_mode);
++ }
+ }
+diff -Naurd mpfr-2.2.1-p2/tests/tset_si.c mpfr-2.2.1-p3/tests/tset_si.c
+--- mpfr-2.2.1-p2/tests/tset_si.c 2006-11-29 09:49:47.000000000 +0000
++++ mpfr-2.2.1-p3/tests/tset_si.c 2007-02-15 23:40:21.000000000 +0000
+@@ -69,6 +69,27 @@
+ if (mpfr_cmp_ui (x, 1<<13) || res <= 0)
+ ERROR ("Prec 2 + ui_2exp");
+
++ mpfr_clear_flags ();
++ mpfr_set_ui_2exp (x, 17, MPFR_EMAX_MAX, GMP_RNDN);
++ if (!mpfr_inf_p (x) || MPFR_IS_NEG (x))
++ ERROR ("mpfr_set_ui_2exp and overflow (bad result)");
++ if (!mpfr_overflow_p ())
++ ERROR ("mpfr_set_ui_2exp and overflow (overflow flag not set)");
++
++ mpfr_clear_flags ();
++ mpfr_set_si_2exp (x, 17, MPFR_EMAX_MAX, GMP_RNDN);
++ if (!mpfr_inf_p (x) || MPFR_IS_NEG (x))
++ ERROR ("mpfr_set_si_2exp (pos) and overflow (bad result)");
++ if (!mpfr_overflow_p ())
++ ERROR ("mpfr_set_si_2exp (pos) and overflow (overflow flag not set)");
++
++ mpfr_clear_flags ();
++ mpfr_set_si_2exp (x, -17, MPFR_EMAX_MAX, GMP_RNDN);
++ if (!mpfr_inf_p (x) || MPFR_IS_POS (x))
++ ERROR ("mpfr_set_si_2exp (neg) and overflow (bad result)");
++ if (!mpfr_overflow_p ())
++ ERROR ("mpfr_set_si_2exp (neg) and overflow (overflow flag not set)");
++
+ mpfr_clear (x);
+ }
+
+@@ -112,7 +133,8 @@
+ unsigned long zl, dl;
+ int inex;
+ int r;
+- mp_exp_t emax;
++ mp_exp_t emin, emax;
++ int flag;
+
+ tests_start_mpfr ();
+
+@@ -350,6 +372,43 @@
+ exit (1);
+ }
+
++ emin = mpfr_get_emin ();
++ mpfr_set_prec (x, 2);
++
++ mpfr_set_emin (4);
++ mpfr_clear_flags ();
++ mpfr_set_ui (x, 7, GMP_RNDU);
++ flag = mpfr_underflow_p ();
++ mpfr_set_emin (emin);
++ if (mpfr_cmp_ui (x, 8) != 0)
++ {
++ printf ("Error for mpfr_set_ui (x, 7, GMP_RNDU), prec = 2, emin = 4\n");
++ exit (1);
++ }
++ if (flag)
++ {
++ printf ("mpfr_set_ui (x, 7, GMP_RNDU) should not underflow "
++ "with prec = 2, emin = 4\n");
++ exit (1);
++ }
++
++ mpfr_set_emin (4);
++ mpfr_clear_flags ();
++ mpfr_set_si (x, -7, GMP_RNDD);
++ flag = mpfr_underflow_p ();
++ mpfr_set_emin (emin);
++ if (mpfr_cmp_si (x, -8) != 0)
++ {
++ printf ("Error for mpfr_set_si (x, -7, GMP_RNDD), prec = 2, emin = 4\n");
++ exit (1);
++ }
++ if (flag)
++ {
++ printf ("mpfr_set_si (x, -7, GMP_RNDD) should not underflow "
++ "with prec = 2, emin = 4\n");
++ exit (1);
++ }
++
+ mpfr_clear (x);
+
+ test_2exp ();
+diff -Naurd mpfr-2.2.1-p2/version.c mpfr-2.2.1-p3/version.c
+--- mpfr-2.2.1-p2/version.c 2007-02-15 23:37:58.000000000 +0000
++++ mpfr-2.2.1-p3/version.c 2007-02-15 23:40:49.000000000 +0000
+@@ -24,5 +24,5 @@
+ const char *
+ mpfr_get_version (void)
+ {
+- return "2.2.1-p2";
++ return "2.2.1-p3";
+ }
--- /dev/null
+diff -Naurd mpfr-2.2.1-p3/VERSION mpfr-2.2.1-p4/VERSION
+--- mpfr-2.2.1-p3/VERSION 2007-02-15 23:40:49.000000000 +0000
++++ mpfr-2.2.1-p4/VERSION 2007-02-15 23:43:29.000000000 +0000
+@@ -1 +1 @@
+-2.2.1-p3
++2.2.1-p4
+diff -Naurd mpfr-2.2.1-p3/mpfr.h mpfr-2.2.1-p4/mpfr.h
+--- mpfr-2.2.1-p3/mpfr.h 2007-02-15 23:40:49.000000000 +0000
++++ mpfr-2.2.1-p4/mpfr.h 2007-02-15 23:43:29.000000000 +0000
+@@ -26,7 +26,7 @@
+ #define MPFR_VERSION_MAJOR 2
+ #define MPFR_VERSION_MINOR 2
+ #define MPFR_VERSION_PATCHLEVEL 1
+-#define MPFR_VERSION_STRING "2.2.1-p3"
++#define MPFR_VERSION_STRING "2.2.1-p4"
+
+ /* Macros dealing with MPFR VERSION */
+ #define MPFR_VERSION_NUM(a,b,c) (((a) << 16L) | ((b) << 8) | (c))
+diff -Naurd mpfr-2.2.1-p3/pow_si.c mpfr-2.2.1-p4/pow_si.c
+--- mpfr-2.2.1-p3/pow_si.c 2006-11-29 09:49:47.000000000 +0000
++++ mpfr-2.2.1-p4/pow_si.c 2007-02-15 23:43:17.000000000 +0000
+@@ -67,24 +67,20 @@
+ /* detect exact powers: x^(-n) is exact iff x is a power of 2 */
+ if (mpfr_cmp_si_2exp (x, MPFR_SIGN(x), MPFR_EXP(x) - 1) == 0)
+ {
+- mp_exp_t expx = MPFR_EXP (x); /* warning: x and y may be the same
+- variable */
+- mpfr_set_si (y, (n % 2) ? MPFR_INT_SIGN(x) : 1, rnd);
+- expx --;
++ mp_exp_t expx = MPFR_EXP (x) - 1, expy;
+ MPFR_ASSERTD (n < 0);
+- /* Warning n*expx may overflow!
+- Some systems abort with LONG_MIN / 1 or LONG_MIN/-1*/
+- if (n != -1 && expx > 0 && -expx < MPFR_EXP_MIN / (-n))
+- MPFR_EXP (y) = MPFR_EMIN_MIN - 1; /* Underflow */
+- else if (n != -1 && expx < 0 && -expx > MPFR_EXP_MAX / (-n))
+- MPFR_EXP (y) = MPFR_EMAX_MAX + 1; /* Overflow */
+- else
+- MPFR_EXP (y) += n * expx;
+- return mpfr_check_range (y, 0, rnd);
++ /* Warning: n * expx may overflow!
++ Some systems (apparently alpha-freebsd) abort with
++ LONG_MIN / 1, and LONG_MIN / -1 is undefined. */
++ expy =
++ n != -1 && expx > 0 && expx > (__gmpfr_emin - 1) / n ?
++ MPFR_EMIN_MIN - 2 /* Underflow */ :
++ n != -1 && expx < 0 && expx < (__gmpfr_emax - 1) / n ?
++ MPFR_EMAX_MAX /* Overflow */ : n * expx;
++ return mpfr_set_si_2exp (y, n % 2 ? MPFR_INT_SIGN (x) : 1,
++ expy, rnd);
+ }
+
+- n = -n;
+-
+ /* General case */
+ {
+ /* Declaration of the intermediary variable */
+@@ -94,9 +90,12 @@
+ mp_prec_t Nt; /* working precision */
+ mp_exp_t err; /* error */
+ int inexact;
++ unsigned long abs_n;
+ MPFR_SAVE_EXPO_DECL (expo);
+ MPFR_ZIV_DECL (loop);
+
++ abs_n = - (unsigned long) n;
++
+ /* compute the precision of intermediary variable */
+ /* the optimal number of bits : see algorithms.tex */
+ Nt = Ny + 3 + MPFR_INT_CEIL_LOG2 (Ny);
+@@ -109,17 +108,17 @@
+ MPFR_ZIV_INIT (loop, Nt);
+ for (;;)
+ {
+- /* compute 1/(x^n) n>0*/
+- mpfr_pow_ui (t, x, (unsigned long int) n, GMP_RNDN);
++ /* compute 1/(x^n), with n > 0 */
++ mpfr_pow_ui (t, x, abs_n, GMP_RNDN);
+ mpfr_ui_div (t, 1, t, GMP_RNDN);
+- /* FIXME: old code improved, but I think this is still incorrect. */
++ /* FIXME: old code improved, but I think this is still incorrect. */
+ if (MPFR_UNLIKELY (MPFR_IS_ZERO (t)))
+ {
+ MPFR_ZIV_FREE (loop);
+ mpfr_clear (t);
+ MPFR_SAVE_EXPO_FREE (expo);
+ return mpfr_underflow (y, rnd == GMP_RNDN ? GMP_RNDZ : rnd,
+- (unsigned) n & 1 ? MPFR_SIGN (x) :
++ abs_n & 1 ? MPFR_SIGN (x) :
+ MPFR_SIGN_POS);
+ }
+ if (MPFR_UNLIKELY (MPFR_IS_INF (t)))
+@@ -127,8 +126,7 @@
+ MPFR_ZIV_FREE (loop);
+ mpfr_clear (t);
+ MPFR_SAVE_EXPO_FREE (expo);
+- return mpfr_overflow (y, rnd,
+- (unsigned) n & 1 ? MPFR_SIGN (x) :
++ return mpfr_overflow (y, rnd, abs_n & 1 ? MPFR_SIGN (x) :
+ MPFR_SIGN_POS);
+ }
+ /* error estimate -- see pow function in algorithms.ps */
+diff -Naurd mpfr-2.2.1-p3/tests/tpow.c mpfr-2.2.1-p4/tests/tpow.c
+--- mpfr-2.2.1-p3/tests/tpow.c 2006-11-29 09:49:47.000000000 +0000
++++ mpfr-2.2.1-p4/tests/tpow.c 2007-02-15 23:43:17.000000000 +0000
+@@ -23,6 +23,7 @@
+ #include <stdlib.h>
+ #include <float.h>
+ #include <math.h>
++#include <limits.h>
+
+ #include "mpfr-test.h"
+
+@@ -231,6 +232,55 @@
+ mpfr_pow_si (x, x, -2, GMP_RNDN);
+ MPFR_ASSERTN(mpfr_inf_p (x) && mpfr_sgn (x) > 0);
+
++ mpfr_set_si (x, 2, GMP_RNDN);
++ mpfr_pow_si (x, x, LONG_MAX, GMP_RNDN); /* 2^LONG_MAX */
++ if (LONG_MAX > mpfr_get_emax () - 1) /* LONG_MAX + 1 > emax */
++ {
++ MPFR_ASSERTN (mpfr_inf_p (x));
++ }
++ else
++ {
++ MPFR_ASSERTN (mpfr_cmp_si_2exp (x, 1, LONG_MAX));
++ }
++
++ mpfr_set_si (x, 2, GMP_RNDN);
++ mpfr_pow_si (x, x, LONG_MIN, GMP_RNDN); /* 2^LONG_MIN */
++ if (LONG_MIN + 1 < mpfr_get_emin ())
++ {
++ MPFR_ASSERTN (mpfr_zero_p (x));
++ }
++ else
++ {
++ MPFR_ASSERTN (mpfr_cmp_si_2exp (x, 1, LONG_MIN));
++ }
++
++ mpfr_set_si (x, 2, GMP_RNDN);
++ mpfr_pow_si (x, x, LONG_MIN + 1, GMP_RNDN); /* 2^(LONG_MIN+1) */
++ if (mpfr_nan_p (x))
++ {
++ printf ("Error in pow_si(2, LONG_MIN+1): got NaN\n");
++ exit (1);
++ }
++ if (LONG_MIN + 2 < mpfr_get_emin ())
++ {
++ MPFR_ASSERTN (mpfr_zero_p (x));
++ }
++ else
++ {
++ MPFR_ASSERTN (mpfr_cmp_si_2exp (x, 1, LONG_MIN + 1));
++ }
++
++ mpfr_set_si_2exp (x, 1, -1, GMP_RNDN); /* 0.5 */
++ mpfr_pow_si (x, x, LONG_MIN, GMP_RNDN); /* 2^(-LONG_MIN) */
++ if (LONG_MIN < 1 - mpfr_get_emax ()) /* 1 - LONG_MIN > emax */
++ {
++ MPFR_ASSERTN (mpfr_inf_p (x));
++ }
++ else
++ {
++ MPFR_ASSERTN (mpfr_cmp_si_2exp (x, 2, - (LONG_MIN + 1)));
++ }
++
+ mpfr_clear (x);
+ }
+
+@@ -270,6 +320,35 @@
+ }
+
+ static void
++pow_si_long_min (void)
++{
++ mpfr_t x, y, z;
++ int inex;
++
++ mpfr_inits2 (sizeof(long) * CHAR_BIT + 32, x, y, z, (void *) 0);
++ mpfr_set_si_2exp (x, 3, -1, GMP_RNDN); /* 1.5 */
++
++ inex = mpfr_set_si (y, LONG_MIN, GMP_RNDN);
++ MPFR_ASSERTN (inex == 0);
++ mpfr_nextbelow (y);
++ mpfr_pow (y, x, y, GMP_RNDD);
++
++ inex = mpfr_set_si (z, LONG_MIN, GMP_RNDN);
++ MPFR_ASSERTN (inex == 0);
++ mpfr_nextabove (z);
++ mpfr_pow (z, x, z, GMP_RNDU);
++
++ mpfr_pow_si (x, x, LONG_MIN, GMP_RNDN); /* 1.5^LONG_MIN */
++ if (mpfr_cmp (x, y) < 0 || mpfr_cmp (x, z) > 0)
++ {
++ printf ("Error in pow_si_long_min\n");
++ exit (1);
++ }
++
++ mpfr_clears (x, y, z, (void *) 0);
++}
++
++static void
+ check_inexact (mp_prec_t p)
+ {
+ mpfr_t x, y, z, t;
+@@ -751,6 +830,7 @@
+ check_pow_ui ();
+ check_pow_si ();
+ check_special_pow_si ();
++ pow_si_long_min ();
+ for (p = 2; p < 100; p++)
+ check_inexact (p);
+ underflows ();
+diff -Naurd mpfr-2.2.1-p3/version.c mpfr-2.2.1-p4/version.c
+--- mpfr-2.2.1-p3/version.c 2007-02-15 23:40:49.000000000 +0000
++++ mpfr-2.2.1-p4/version.c 2007-02-15 23:43:29.000000000 +0000
+@@ -24,5 +24,5 @@
+ const char *
+ mpfr_get_version (void)
+ {
+- return "2.2.1-p3";
++ return "2.2.1-p4";
+ }
--- /dev/null
+diff -Naurd mpfr-2.2.1-p4/VERSION mpfr-2.2.1-p5/VERSION
+--- mpfr-2.2.1-p4/VERSION 2007-02-15 23:43:29.000000000 +0000
++++ mpfr-2.2.1-p5/VERSION 2007-02-15 23:46:49.000000000 +0000
+@@ -1 +1 @@
+-2.2.1-p4
++2.2.1-p5
+diff -Naurd mpfr-2.2.1-p4/mpfr.h mpfr-2.2.1-p5/mpfr.h
+--- mpfr-2.2.1-p4/mpfr.h 2007-02-15 23:43:29.000000000 +0000
++++ mpfr-2.2.1-p5/mpfr.h 2007-02-15 23:46:49.000000000 +0000
+@@ -26,7 +26,7 @@
+ #define MPFR_VERSION_MAJOR 2
+ #define MPFR_VERSION_MINOR 2
+ #define MPFR_VERSION_PATCHLEVEL 1
+-#define MPFR_VERSION_STRING "2.2.1-p4"
++#define MPFR_VERSION_STRING "2.2.1-p5"
+
+ /* Macros dealing with MPFR VERSION */
+ #define MPFR_VERSION_NUM(a,b,c) (((a) << 16L) | ((b) << 8) | (c))
+diff -Naurd mpfr-2.2.1-p4/mpfr.info mpfr-2.2.1-p5/mpfr.info
+--- mpfr-2.2.1-p4/mpfr.info 2006-11-29 09:51:26.000000000 +0000
++++ mpfr-2.2.1-p5/mpfr.info 2007-02-15 23:46:21.000000000 +0000
+@@ -1,10 +1,10 @@
+ This is mpfr.info, produced by makeinfo version 4.8 from mpfr.texi.
+
+ This manual documents how to install and use the Multiple Precision
+-Floating-Point Reliable Library, version 2.2.1.
++Floating-Point Reliable Library, version 2.2.2-p5.
+
+ Copyright 1991, 1993, 1994, 1995, 1996, 1997, 1998, 1999, 2000,
+-2001, 2002, 2003, 2004, 2005, 2006 Free Software Foundation, Inc.
++2001, 2002, 2003, 2004, 2005, 2006, 2007 Free Software Foundation, Inc.
+
+ Permission is granted to copy, distribute and/or modify this
+ document under the terms of the GNU Free Documentation License, Version
+@@ -26,10 +26,10 @@
+ ****
+
+ This manual documents how to install and use the Multiple Precision
+-Floating-Point Reliable Library, version 2.2.1.
++Floating-Point Reliable Library, version 2.2.2-p5.
+
+ Copyright 1991, 1993, 1994, 1995, 1996, 1997, 1998, 1999, 2000,
+-2001, 2002, 2003, 2004, 2005, 2006 Free Software Foundation, Inc.
++2001, 2002, 2003, 2004, 2005, 2006, 2007 Free Software Foundation, Inc.
+
+ Permission is granted to copy, distribute and/or modify this
+ document under the terms of the GNU Free Documentation License, Version
+@@ -173,7 +173,7 @@
+ This will copy the files `mpfr.h' and `mpf2mpfr.h' to the directory
+ `/usr/local/include', the file `libmpfr.a' to the directory
+ `/usr/local/lib', and the file `mpfr.info' to the directory
+- `/usr/local/info' (or if you passed the `--prefix' option to
++ `/usr/local/share/info' (or if you passed the `--prefix' option to
+ `configure', using the prefix directory given as argument to
+ `--prefix' instead of `/usr/local').
+
+@@ -396,6 +396,53 @@
+ should return an overflow or an underflow if `1' is not representable
+ in the current exponent range.
+
++4.5 Exceptions
++==============
++
++MPFR supports 5 exception types:
++
++ * Underflow: An underflow occurs when the exact result of a function
++ is a non-zero real number and the result obtained after the
++ rounding, assuming an unbounded exponent range (for the rounding),
++ has an exponent smaller than the minimum exponent of the current
++ range. In the round-to-nearest mode, the halfway case is rounded
++ toward zero.
++
++ Note: This is not the single definition of the underflow. MPFR
++ chooses to consider the underflow after rounding. The underflow
++ before rounding can also be defined. For instance, consider a
++ function that has the exact result 7 multiplied by two to the power
++ E-4, where E is the smallest exponent (for a mantissa between 1/2
++ and 1) in the current range, with a 2-bit target precision and
++ rounding towards plus infinity. The exact result has the exponent
++ E-1. With the underflow before rounding, such a function call
++ would yield an underflow, as E-1 is outside the current exponent
++ range. However, MPFR first considers the rounded result assuming
++ an unbounded exponent range. The exact result cannot be
++ represented exactly in precision 2, and here, it is rounded to 0.5
++ times 2 to E, which is representable in the current exponent
++ range. As a consequence, this will not yield an underflow in MPFR.
++
++ * Overflow: An overflow occurs when the exact result of a function
++ is a non-zero real number and the result obtained after the
++ rounding, assuming an unbounded exponent range (for the rounding),
++ has an exponent larger than the maximum exponent of the current
++ range. In the round-to-nearest mode, the result is infinite.
++
++ * Invalid (NaN): An invalid (or NaN) exception occurs when the
++ result of a function is a NaN.
++
++ * Inexact: An inexact exception occurs when the result of a function
++ cannot be represented exactly and must be rounded.
++
++ * Range error: A range exception occurs when a function that does
++ not return a MPFR number (such as comparisons and conversions to
++ an integer) has an invalid result.
++
++
++ MPFR has a global flag for each exception, which can be cleared, set
++or tested by functions described in *Note Exceptions::.
++
+ \1f
+ File: mpfr.info, Node: MPFR Interface, Next: Contributors, Prev: MPFR Basics, Up: Top
+
+@@ -950,6 +997,12 @@
+ Set ROP to the absolute value of OP, rounded in the direction RND.
+ Just changes the sign if ROP and OP are the same variable.
+
++ -- Function: int mpfr_dim (mpfr_t ROP, mpfr_t OP1, mpfr_t OP2,
++ mp_rnd_t RND)
++ Set ROP to the positive difference of OP1 and OP2, i.e., OP1 - OP2
++ rounded in the direction RND if OP1 > OP2, and +0 otherwise.
++ Returns NaN when OP1 or OP2 is NaN.
++
+ -- Function: int mpfr_mul_2ui (mpfr_t ROP, mpfr_t OP1, unsigned long
+ int OP2, mp_rnd_t RND)
+ -- Function: int mpfr_mul_2si (mpfr_t ROP, mpfr_t OP1, long int OP2,
+@@ -1873,9 +1926,12 @@
+ Meunier helped in the design of the `mpfr_erf' code.
+
+ The development of the MPFR library would not have been possible
+-without the continuous support of LORIA, INRIA and INRIA Lorraine. The
+-development of MPFR was also supported by a grant (202F0659 00 MPN 121)
+-from the Conseil Régional de Lorraine in 2002.
++without the continuous support of INRIA, and of the LORIA and LIP
++laboratories. In particular the main authors were or are members of the
++PolKA, Spaces, Cacao project-teams at LORIA (Nancy, France) and of the
++Arenaire project-team at LIP (Lyon, France). The development of MPFR
++was also supported by a grant (202F0659 00 MPN 121) from the Conseil
++Régional de Lorraine in 2002.
+
+ \1f
+ File: mpfr.info, Node: References, Next: GNU Free Documentation License, Prev: Contributors, Up: Top
+@@ -2387,8 +2443,8 @@
+ * Precision <1>: MPFR Interface. (line 20)
+ * Precision: MPFR Basics. (line 26)
+ * Reporting bugs: Reporting Bugs. (line 6)
+-* Rounding Modes <1>: Rounding Modes. (line 6)
+-* Rounding Modes: MPFR Basics. (line 32)
++* Rounding Modes <1>: MPFR Basics. (line 32)
++* Rounding Modes: Rounding Modes. (line 6)
+ * Special functions: Special Functions. (line 6)
+
+ \1f
+@@ -2477,14 +2533,16 @@
+ * mpfr_custom_init_set: Custom Interface. (line 47)
+ * mpfr_custom_move: Custom Interface. (line 84)
+ * MPFR_DECL_INIT: Advanced Functions. (line 10)
++* mpfr_dim: Basic Arithmetic Functions.
++ (line 155)
+ * mpfr_div: Basic Arithmetic Functions.
+ (line 58)
+ * mpfr_div_2exp: Compatibility with MPF.
+ (line 47)
+ * mpfr_div_2si: Basic Arithmetic Functions.
+- (line 164)
++ (line 170)
+ * mpfr_div_2ui: Basic Arithmetic Functions.
+- (line 162)
++ (line 168)
+ * mpfr_div_q: Basic Arithmetic Functions.
+ (line 70)
+ * mpfr_div_si: Basic Arithmetic Functions.
+@@ -2625,9 +2683,9 @@
+ * mpfr_mul_2exp: Compatibility with MPF.
+ (line 45)
+ * mpfr_mul_2si: Basic Arithmetic Functions.
+- (line 157)
++ (line 163)
+ * mpfr_mul_2ui: Basic Arithmetic Functions.
+- (line 155)
++ (line 161)
+ * mpfr_mul_q: Basic Arithmetic Functions.
+ (line 48)
+ * mpfr_mul_si: Basic Arithmetic Functions.
+@@ -2808,33 +2866,33 @@
+
+ \1f
+ Tag Table:
+-Node: Top\7f949
+-Node: Copying\7f2260
+-Node: Introduction to MPFR\7f3977
+-Node: Installing MPFR\7f5675
+-Node: Reporting Bugs\7f8343
+-Node: MPFR Basics\7f9895
+-Node: MPFR Interface\7f16074
+-Node: Initialization Functions\7f18238
+-Node: Assignment Functions\7f21660
+-Node: Combined Initialization and Assignment Functions\7f28550
+-Node: Conversion Functions\7f29832
+-Node: Basic Arithmetic Functions\7f35481
+-Node: Comparison Functions\7f42745
+-Node: Special Functions\7f46215
+-Node: Input and Output Functions\7f55364
+-Node: Integer Related Functions\7f57300
+-Node: Miscellaneous Functions\7f60294
+-Node: Rounding Modes\7f64675
+-Node: Exceptions\7f66024
+-Node: Advanced Functions\7f71338
+-Node: Compatibility with MPF\7f73943
+-Node: Custom Interface\7f76212
+-Node: Internals\7f80255
+-Node: Contributors\7f82275
+-Node: References\7f84024
+-Node: GNU Free Documentation License\7f84769
+-Node: Concept Index\7f107212
+-Node: Function Index\7f111074
++Node: Top\7f958
++Node: Copying\7f2278
++Node: Introduction to MPFR\7f3995
++Node: Installing MPFR\7f5693
++Node: Reporting Bugs\7f8367
++Node: MPFR Basics\7f9919
++Node: MPFR Interface\7f18384
++Node: Initialization Functions\7f20548
++Node: Assignment Functions\7f23970
++Node: Combined Initialization and Assignment Functions\7f30860
++Node: Conversion Functions\7f32142
++Node: Basic Arithmetic Functions\7f37791
++Node: Comparison Functions\7f45323
++Node: Special Functions\7f48793
++Node: Input and Output Functions\7f57942
++Node: Integer Related Functions\7f59878
++Node: Miscellaneous Functions\7f62872
++Node: Rounding Modes\7f67253
++Node: Exceptions\7f68602
++Node: Advanced Functions\7f73916
++Node: Compatibility with MPF\7f76521
++Node: Custom Interface\7f78790
++Node: Internals\7f82833
++Node: Contributors\7f84853
++Node: References\7f86790
++Node: GNU Free Documentation License\7f87535
++Node: Concept Index\7f109978
++Node: Function Index\7f113840
+ \1f
+ End Tag Table
+diff -Naurd mpfr-2.2.1-p4/mpfr.texi mpfr-2.2.1-p5/mpfr.texi
+--- mpfr-2.2.1-p4/mpfr.texi 2006-11-29 09:49:47.000000000 +0000
++++ mpfr-2.2.1-p5/mpfr.texi 2007-02-15 23:46:09.000000000 +0000
+@@ -2,8 +2,8 @@
+ @c %**start of header
+ @setfilename mpfr.info
+ @documentencoding ISO-8859-1
+-@set VERSION 2.2.1
+-@set UPDATED-MONTH November 2006
++@set VERSION 2.2.2-p5
++@set UPDATED-MONTH February 2007
+ @settitle MPFR @value{VERSION}
+ @synindex tp fn
+ @iftex
+@@ -15,7 +15,7 @@
+ This manual documents how to install and use the Multiple Precision
+ Floating-Point Reliable Library, version @value{VERSION}.
+
+-Copyright 1991, 1993, 1994, 1995, 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006 Free Software Foundation, Inc.
++Copyright 1991, 1993, 1994, 1995, 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007 Free Software Foundation, Inc.
+
+ Permission is granted to copy, distribute and/or modify this document under
+ the terms of the GNU Free Documentation License, Version 1.1 or any later
+@@ -312,7 +312,7 @@
+ This will copy the files @file{mpfr.h} and @file{mpf2mpfr.h} to the directory
+ @file{/usr/local/include}, the file @file{libmpfr.a} to the directory
+ @file{/usr/local/lib}, and the file @file{mpfr.info} to the directory
+-@file{/usr/local/info} (or if you passed the @samp{--prefix} option to
++@file{/usr/local/share/info} (or if you passed the @samp{--prefix} option to
+ @file{configure}, using the prefix directory given as argument to
+ @samp{--prefix} instead of @file{/usr/local}).
+ @end enumerate
+@@ -559,6 +559,61 @@
+ for special cases (like @code{acos(0)}) should return an overflow or
+ an underflow if @code{1} is not representable in the current exponent range.
+
++@section Exceptions
++
++MPFR supports 5 exception types:
++
++@itemize @bullet
++
++@item Underflow:
++An underflow occurs when the exact result of a function is a non-zero
++real number and the result obtained after the rounding, assuming an
++unbounded exponent range (for the rounding), has an exponent smaller
++than the minimum exponent of the current range. In the round-to-nearest
++mode, the halfway case is rounded toward zero.
++
++Note: This is not the single definition of the underflow. MPFR chooses
++to consider the underflow after rounding. The underflow before rounding
++can also be defined. For instance, consider a function that has the
++exact result @m{7 \times 2^{e-4}, 7 multiplied by two to the power
++@var{e}@minus{}4}, where @var{e} is the smallest exponent (for a
++mantissa between 1/2 and 1) in the current
++range, with a 2-bit target precision and rounding towards plus infinity.
++The exact result has the exponent @var{e}@minus{}1. With the underflow
++before rounding, such a function call would yield an underflow, as
++@var{e}@minus{}1 is outside the current exponent range. However, MPFR
++first considers the rounded result assuming an unbounded exponent range.
++The exact result cannot be represented exactly in precision 2, and here,
++it is rounded to @m{0.5 @times 2^e, 0.5 times 2 to @var{e}}, which is
++representable in the current exponent range. As a consequence, this will
++not yield an underflow in MPFR.
++
++@item Overflow:
++An overflow occurs when the exact result of a function is a non-zero
++real number and the result obtained after the rounding, assuming an
++unbounded exponent range (for the rounding), has an exponent larger
++than the maximum exponent of the current range. In the round-to-nearest
++mode, the result is infinite.
++
++@item Invalid (NaN):
++An invalid (or NaN) exception occurs when the result of a function is
++a NaN.
++@c NaN is defined above. So, we don't say anything more.
++
++@item Inexact:
++An inexact exception occurs when the result of a function cannot be
++represented exactly and must be rounded.
++
++@item Range error:
++A range exception occurs when a function that does not return a MPFR
++number (such as comparisons and conversions to an integer) has an
++invalid result.
++
++@end itemize
++
++MPFR has a global flag for each exception, which can be cleared, set
++or tested by functions described in @ref{Exceptions}.
++
+ @node MPFR Interface, Contributors, MPFR Basics, Top
+ @comment node-name, next, previous, up
+ @chapter MPFR Interface
+@@ -1130,6 +1185,13 @@
+ Just changes the sign if @var{rop} and @var{op} are the same variable.
+ @end deftypefun
+
++@deftypefun int mpfr_dim (mpfr_t @var{rop}, mpfr_t @var{op1}, mpfr_t @var{op2}, mp_rnd_t @var{rnd})
++Set @var{rop} to the positive difference of @var{op1} and @var{op2}, i.e.,
++@math{@var{op1} - @var{op2}} rounded in the direction @var{rnd}
++if @math{@var{op1} > @var{op2}}, and +0 otherwise.
++Returns NaN when @var{op1} or @var{op2} is NaN.
++@end deftypefun
++
+ @deftypefun int mpfr_mul_2ui (mpfr_t @var{rop}, mpfr_t @var{op1}, unsigned long int @var{op2}, mp_rnd_t @var{rnd})
+ @deftypefunx int mpfr_mul_2si (mpfr_t @var{rop}, mpfr_t @var{op1}, long int @var{op2}, mp_rnd_t @var{rnd})
+ Set @var{rop} to @m{@var{op1} \times 2^{op2}, @var{op1} times 2 raised
+@@ -2196,7 +2258,10 @@
+ Ludovic Meunier helped in the design of the @code{mpfr_erf} code.
+
+ The development of the MPFR library would not have been possible without the
+-continuous support of LORIA, INRIA and INRIA Lorraine.
++continuous support of INRIA, and of the LORIA and LIP laboratories.
++In particular the main authors were or are members of the
++PolKA, Spaces, Cacao project-teams at LORIA (Nancy, France)
++and of the Arenaire project-team at LIP (Lyon, France).
+ The development of MPFR was also supported by a grant
+ (202F0659 00 MPN 121) from the Conseil R@'egional de Lorraine in 2002.
+
+diff -Naurd mpfr-2.2.1-p4/version.c mpfr-2.2.1-p5/version.c
+--- mpfr-2.2.1-p4/version.c 2007-02-15 23:43:29.000000000 +0000
++++ mpfr-2.2.1-p5/version.c 2007-02-15 23:46:49.000000000 +0000
+@@ -24,5 +24,5 @@
+ const char *
+ mpfr_get_version (void)
+ {
+- return "2.2.1-p4";
++ return "2.2.1-p5";
+ }
--- /dev/null
+MD5 40bf06f8081461d8db7d6f4ad5b9f6bd mpfr-2.2.1.tar.bz2 787634
+RMD160 cb4c24dc356412e3bae30242f426a2a02d146d04 mpfr-2.2.1.tar.bz2 787634
+SHA256 92e2fea0eefd3d50b8e762334326657bfc5be9a6bda58cfa724b2a5cba9f65e4 mpfr-2.2.1.tar.bz2 787634
--- /dev/null
+# Copyright 1999-2007 Gentoo Foundation
+# Distributed under the terms of the GNU General Public License v2
+# $Header: /var/cvsroot/gentoo-x86/dev-libs/mpfr/mpfr-2.2.1_p5.ebuild,v 1.1 2007/04/06 14:07:46 dragonheart Exp $
+
+# NOTE: we cannot depend on autotools here starting with gcc-4.3.x
+inherit eutils
+
+MY_PV=${PV/_p*}
+MY_P=${PN}-${MY_PV}
+PLEVEL=${PV/*p}
+DESCRIPTION="library for multiple-precision floating-point computations with exact rounding"
+HOMEPAGE="http://www.mpfr.org/"
+SRC_URI="http://www.mpfr.org/mpfr-current/${MY_P}.tar.bz2"
+
+LICENSE="LGPL-2.1"
+SLOT="0"
+KEYWORDS="~alpha ~amd64 ~arm ~hppa ~ia64 ~mips ~ppc ~ppc64 ~s390 ~sh ~sparc ~sparc-fbsd ~x86 ~x86-fbsd"
+IUSE=""
+
+DEPEND=">=dev-libs/gmp-4.1.4-r2"
+
+S=${WORKDIR}/${MY_P}
+
+src_unpack() {
+ unpack "${MY_P}.tar.bz2"
+ cd "${S}"
+ [[ ${PLEVEL} == ${PV} ]] && return 0
+ for ((i=1; i<=PLEVEL; ++i)) ; do
+ patch=patch$(printf '%02d' ${i})
+ if [[ -f ${FILESDIR}/${MY_PV}/${patch} ]] ; then
+ epatch "${FILESDIR}"/${MY_PV}/${patch}
+ elif [[ -f ${DISTDIR}/${PN}-${MY_PV}_p${i} ]] ; then
+ epatch "${DISTDIR}"/${PN}-${MY_PV}_p${i}
+ else
+ ewarn "${DISTDIR}/${PN}-${MY_PV}_p${i}"
+ die "patch ${i} missing - please report to bugs.gentoo.org"
+ fi
+ done
+}
+
+src_compile() {
+ econf \
+ --enable-shared \
+ --enable-static \
+ || die
+ emake || die
+}
+
+src_install() {
+ emake install DESTDIR="${D}" || die
+ dodoc AUTHORS BUGS ChangeLog NEWS README TODO
+ dohtml *.html
+}
projects / gentoo.git / commitdiff