[Bug libquadmath/63488] large errors in y0q

[fxcoudert at gcc dot gnu.org]

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