]> git.ipfire.org Git - thirdparty/glibc.git/commit - sysdeps/ieee754/ldbl-128/e_j0l.c
Fix y0 and y1 exception handling for zero input [BZ #21134]
authorGabriel F. T. Gomes <gftg@linux.vnet.ibm.com>
Mon, 13 Feb 2017 00:36:27 +0000 (22:36 -0200)
committerGabriel F. T. Gomes <gftg@linux.vnet.ibm.com>
Wed, 15 Feb 2017 12:30:59 +0000 (10:30 -0200)
commit4918e5f4cd63290bb0b2c614f52092ca6a779126
tree1fde6c355b58f36b89fee9dddecc60a7fb29216e
parent10303eb74bfe33d46ef167d2ea31c746ea1cd6ad
Fix y0 and y1 exception handling for zero input [BZ #21134]

The Bessel functions of the second type (Yn) should raise the "divide
by zero" exception when input is zero (both positive and negative).
Current code gives the right output, but fails to set the exception.
This error is exposed for float, double, and 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, set the exception as expected.

Tested for powerpc64le.

[BZ #21134]
* sysdeps/ieee754/dbl-64/e_j0.c (__ieee754_y0): Raise the
"divide by zero" exception when the input is zero.
* sysdeps/ieee754/dbl-64/e_j1.c (__ieee754_y1): Likewise.
* sysdeps/ieee754/flt-32/e_j0f.c (__ieee754_y0f): Likewise.
* sysdeps/ieee754/flt-32/e_j1f.c (__ieee754_y1f): Likewise.
* sysdeps/ieee754/ldbl-128/e_j0l.c (__ieee754_y0l): Likewise.
* sysdeps/ieee754/ldbl-128/e_j1l.c (__ieee754_y1l): Likewise.
ChangeLog
sysdeps/ieee754/dbl-64/e_j0.c
sysdeps/ieee754/dbl-64/e_j1.c
sysdeps/ieee754/flt-32/e_j0f.c
sysdeps/ieee754/flt-32/e_j1f.c
sysdeps/ieee754/ldbl-128/e_j0l.c
sysdeps/ieee754/ldbl-128/e_j1l.c