ldbl-128: Fix y0 and y1 for -Inf input [BZ #21130]
authorGabriel F. T. Gomes <gftg@linux.vnet.ibm.com>
Fri, 6 Jan 2017 15:05:30 +0000 (13:05 -0200)
committerGabriel F. T. Gomes <gftg@linux.vnet.ibm.com>
Sun, 12 Feb 2017 20:30:38 +0000 (18:30 -0200)
commitb987917e6aa7ffe2fd74f0b6a989438e6edd0727
treeb8b52e949fb772130ffa853ec1f1e0a2944c87ff
parent40b7fbf27310d4f4a899aa2c54740dce5369ec70
ldbl-128: Fix y0 and y1 for -Inf input [BZ #21130]

The Bessel functions of the second type (Yn) are not defined for
negative input and should return NAN with the "invalid" exception
raised, in these cases.  However, current code checks for infinity and
return zero, regardless of the sign.  This error is exposed for long
double when linking with -lieee.  Without this flag, the error is not
exposed, because the wrappers for these functions, which use
__kernel_standard functionality, return the correct value.

Tested for powerpc64le.

[BZ #21130]
* sysdeps/ieee754/ldbl-128/e_j0l.c (__ieee754_y0l): Return NAN
with the "invalid" exception raised when x is -Inf.
* sysdeps/ieee754/ldbl-128/e_j1l.c (__ieee754_y1l): Likewise.
ChangeLog
sysdeps/ieee754/ldbl-128/e_j0l.c
sysdeps/ieee754/ldbl-128/e_j1l.c