https://gcc.gnu.org/bugzilla/show_bug.cgi?id=63488
Francois-Xavier Coudert <fxcoudert at gcc dot gnu.org> changed:
What |Removed |Added
----------------------------------------------------------------------------
Status|UNCONFIRMED |NEW
Last reconfirmed| |2014-10-08
Version|unknown |5.0
Ever confirmed|0 |1
--- Comment #1 from Francois-Xavier Coudert <fxcoudert at gcc dot gnu.org> ---
Confirmed on trunk.
Also visible with the following Fortran source:
real(kind=16), parameter :: x = 0.8935761195587725798762818805462843676_16
real(kind=16) :: y
print *, bessel_y0(x)
y = x ; y = bessel_y0(y)
print *, y
print *, (y - bessel_y0(x)) / bessel_y0(x)
end
which outputs:
-7.44623881999333920107530266264973398E-0007
-7.44623881999333920107530266281566970E-0007
2.22845016536952266229206331898288934E-0029
