Fix ldbl-128 erfl spurious underflows (bug 16287).
authorJoseph Myers <joseph@codesourcery.com>
Tue, 24 Jun 2014 20:56:56 +0000 (20:56 +0000)
committerJoseph Myers <joseph@codesourcery.com>
Tue, 24 Jun 2014 20:56:56 +0000 (20:56 +0000)
commite7dd3c8c1db8c6d293abb995e033893b7328ad19
treedcfbd9476c52dbccdedf8d691b37b7442199a756
parent0ea5d373c7984a1da2d0e990c6e78e6d6c24f91b
Fix ldbl-128 erfl spurious underflows (bug 16287).

This patch fixes bug 16287, spurious underflows from ldbl-128 erfl
arising from it calling erfcl for arguments with absolute value at
least 1.0, although for large positive arguments erfcl correctly
underflows but erfl shouldn't.  The fix is simply to avoid calling
erfcl, and just return 1, for arguments above a cut-off large enough
that erfl correctly rounds to-nearest as 1 but not so large that erfcl
underflows.

Tested mips64.  Also tested x86_64 and x86 to confirm the new tests
(taken from the tests of erfc) don't cause any problems there; no ulps
updates needed.

[BZ #16287]
* sysdeps/ieee754/ldbl-128/s_erfl.c (__erfl): Return 1 without
calling __erfcl for arguments at least 16.
* math/auto-libm-test-in: Add more tests of erf.
* math/auto-libm-test-out: Regenerated.
ChangeLog
NEWS
math/auto-libm-test-in
math/auto-libm-test-out
sysdeps/ieee754/ldbl-128/s_erfl.c