Log message:
Update MPFR to version 3.0.1 patch level 3.
Patch 3 fixes the following bug:
The mpfr_atan function calls mpfr_check_range on the unrounded
result (no longer allocated) instead of the rounded result.
Consequence of this bug: possible incorrect or invalid result and/or
incorrect flags if the unrounded result or rounded result is not in
the current exponent range. Or possibly worse if the memory used by
the unrounded result is modified before the return of the function.
This bug has been present since MPFR 2.1.0 (exceptions were not
supported before).
Patch 2 fixes the following bug:
The mpfr_rec_sqrt function can have undefined behavior due to the
call of mpn_add_1 or mpn_sub_1 (from GMP) with a null size, which is
not allowed. The fourth argument of mpn_add_1 or mpn_sub_1 (i.e. the
value that should have been added or subtracted) is here also 0, so
that the consequences are limited. In practice, with the current GMP
version, this yields a small buffer overflow (1 limb, i.e. 32 bits
or 64 bits), as shown by valgrind, with the possible consequences:
a crash, or memory corruption (very unlikely) if another process or
thread modifies the limb at the same time (as in this context, the
value is first read then written back by GMP). This bug has always
been present since the introduction of the mpfr_rec_sqrt function
in MPFR 2.4.0.
Patch 1 fixes the following bug:
If the exponent range has been reduced such that emax = 1 or 2,
mpfr_asin may give an incorrect result on +1 and -1.
which should probably never occur in practical codes.
|
Log message:
Update MPFR to version 3.0.0.
Contributed by Stathis Kamperis on IRC.
Changes from versions 2.4.* to version 3.0.0:
- The "boudin aux pommes" release.
- MPFR 3.0.0 is binary incompatible with previous versions but (almost)
API compatible. More precisely the obsolete functions mpfr_random
and mpfr_random2 have been removed, the meaning of the return type
of the function mpfr_get_f has changed, and the return type of the
function mpfr_get_z is now int instead of void. In practice, this
should not break any existing code.
- MPFR is now distributed under the GNU Lesser General Public License
version 3 or later (LGPL v3+).
- Rounding modes GMP_RNDx are now MPFR_RNDx (GMP_RNDx kept for
compatibility).
- A new rounding mode (MPFR_RNDA) is available to round away from zero.
- The rounding mode type is now mpfr_rnd_t (as in previous versions,
both mpfr_rnd_t and mp_rnd_t are accepted, but mp_rnd_t may be
removed in the future).
- The precision type is now mpfr_prec_t (as in previous versions, both
mpfr_prec_t and mp_prec_t are accepted, but mp_prec_t may be removed
in the future) and it is now signed (it was unsigned in MPFR 2.*, but
this was not documented). In practice, this change should not affect
existing code that assumed nothing on the precision type.
- MPFR now has its own exponent type mpfr_exp_t, which is currently
the same as GMP's mp_exp_t.
- Functions mpfr_random and mpfr_random2 have been removed.
- mpfr_get_f and mpfr_get_z now return a ternary value.
- mpfr_strtofr now accepts bases from 37 to 62.
- mpfr_custom_get_mantissa was renamed to mpfr_custom_get_significand
(mpfr_custom_get_mantissa is still available via a #define).
- Functions mpfr_get_si, mpfr_get_ui, mpfr_get_sj, mpfr_get_uj,
mpfr_get_z and mpfr_get_z_2exp no longer have cases with undefined
behavior; in these cases, the behavior is now specified, and in
particular, the erange flag is set.
- New functions mpfr_buildopt_tls_p and mpfr_buildopt_decimal_p giving
information about options used at MPFR build time.
- New function mpfr_regular_p.
- New function mpfr_set_zero.
- New function mpfr_digamma.
- New function mpfr_ai (incomplete, experimental).
- New functions mpfr_set_flt and mpfr_get_flt to convert from/to the
float type.
- New function mpfr_urandom.
- New function mpfr_set_z_2exp (companion to mpfr_get_z_2exp, which
was renamed from mpfr_get_z_exp in previous versions).
- Speed improvement for large operands in the trigonometric functions
(mpfr_sin, mpfr_cos, mpfr_tan, mpfr_sin_cos): speedup of about 2.5
for 10^5 digits, of about 5 for 10^6 digits.
- Speed improvement for large operands of the inverse trigonometric
functions (arcsin, arccos, arctan): about 2 for 10^3 digits, up to
2.7 for 10^6 digits.
- Some documentation files are installed in $docdir.
- The detection of a GMP build directory (more precisely, the internal
header files of GMP) was previously done separately from the use of
the --with-gmp-build configure option. This was not consistent with
the documentation and with other parts of the configure script. So,
as of MPFR 3.0.0, the internal header files of GMP are now used if
and only if the --with-gmp-build configure option is given.
- The configure script recognizes some extra "long double" formats
(double big endian, double little endian, double-double big endian).
- MPFR manual: added "API Compatibility" section.
- Test coverage: 97.1% lines of code.
- Bug fixes.
|