This is the mail archive of the
glibc-bugs@sourceware.org
mailing list for the glibc project.
[Bug math/16338] ldbl-128 logl does not handle subnormals
- From: "cvs-commit at gcc dot gnu.org" <sourceware-bugzilla at sourceware dot org>
- To: glibc-bugs at sourceware dot org
- Date: Wed, 18 Dec 2013 11:40:57 +0000
- Subject: [Bug math/16338] ldbl-128 logl does not handle subnormals
- Auto-submitted: auto-generated
- References: <bug-16338-131 at http dot sourceware dot org/bugzilla/>
http://sourceware.org/bugzilla/show_bug.cgi?id=16338
--- Comment #1 from cvs-commit at gcc dot gnu.org <cvs-commit at gcc dot gnu.org> ---
This is an automated email from the git hooks/post-receive script. It was
generated because a ref change was pushed to the repository containing
the project "GNU C Library master sources".
The branch, master has been updated
via 2dec468fd82c6d968f622c2c519fcb8dc09bddac (commit)
from 4d84e6addd62bdc256627af334ba85489112e008 (commit)
Those revisions listed above that are new to this repository have
not appeared on any other notification email; so we list those
revisions in full, below.
- Log -----------------------------------------------------------------
https://sourceware.org/git/gitweb.cgi?p=glibc.git;h=2dec468fd82c6d968f622c2c519fcb8dc09bddac
commit 2dec468fd82c6d968f622c2c519fcb8dc09bddac
Author: Joseph Myers <joseph@codesourcery.com>
Date: Wed Dec 18 11:38:27 2013 +0000
Fix ldbl-128 logl for subnormals (bug 16338).
This patch fixes bug 16338, ldbl-128 logl not handling subnormals
(with consequent inaccuracy for lgammal as well). The fix is simply
to use __frexpl when determining the exponent, as done already in
log2l and log10l. Given the lack of testing of small arguments to any
of the log* functions, appropriate tests are added for all of them.
Tested x86_64 and x86 and ulps updated accordingly, and spot tests
also run for mips64 to confirm the ldbl-128 fix.
Note that while this fixes lgammal inaccuracy for small positive
arguments, I suspect that there will still be problems with spurious
underflows in that case.
* sysdeps/ieee754/ldbl-128/e_logl.c (__ieee754_logl): Use __frexpl
to determine exponent and adjust argument to have exponent of -1.
* math/auto-libm-test-in: Add more tests of log, log10, log1p and
log2.
* math/auto-libm-test-out: Regenerated.
* sysdeps/x86_64/fpu/libm-test-ulps: Update.
-----------------------------------------------------------------------
Summary of changes:
ChangeLog | 11 +
NEWS | 2 +-
math/auto-libm-test-in | 11 +
math/auto-libm-test-out | 790 +++++++++++++++++++++++++++++++++++++
sysdeps/i386/fpu/libm-test-ulps | 3 +
sysdeps/ieee754/ldbl-128/e_logl.c | 5 +-
sysdeps/x86_64/fpu/libm-test-ulps | 3 +
7 files changed, 821 insertions(+), 4 deletions(-)
--
You are receiving this mail because:
You are on the CC list for the bug.