[PATCH v2] Fix y0 and y1 exception handling for zero input [BZ #21134]

[Joseph Myers]

On Tue, 14 Feb 2017, Gabriel F. T. Gomes wrote:

> 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.
> 
> 2017-02-12  Gabriel F. T. Gomes  <gftg@linux.vnet.ibm.com>
> 
> 	[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.

OK.

-- 
Joseph S. Myers
joseph@codesourcery.com