https://gcc.gnu.org/bugzilla/show_bug.cgi?id=63488
--- Comment #8 from Steve Kargl ---
On Thu, Oct 09, 2014 at 06:27:08PM +, zimmerma+gcc at loria dot fr wrote:
> https://gcc.gnu.org/bugzilla/show_bug.cgi?id=63488
>
> --- Comment #7 from Paul Zimmermann ---
> I agree that near zeroes we can
https://gcc.gnu.org/bugzilla/show_bug.cgi?id=63488
--- Comment #7 from Paul Zimmermann ---
I agree that near zeroes we can expect large errors. However for other
functions
I got only small errors in ulps, maybe I was unlucky. Also the ultimate goal is
to get correct rounding, even near zeroes.
https://gcc.gnu.org/bugzilla/show_bug.cgi?id=63488
--- Comment #6 from Steve Kargl ---
On Thu, Oct 09, 2014 at 05:34:11PM +, kargl at gcc dot gnu.org wrote:
>
> Testing 1 values in a small interval about the lowest
> 10 zeros, the double precision y0() on FreeBSD (which comes
> from fdli
https://gcc.gnu.org/bugzilla/show_bug.cgi?id=63488
kargl at gcc dot gnu.org changed:
What|Removed |Added
CC||kargl at gcc dot gnu.org
--- C
https://gcc.gnu.org/bugzilla/show_bug.cgi?id=63488
--- Comment #4 from Francois-Xavier Coudert ---
(In reply to jos...@codesourcery.com from comment #3)
> (Eventually I think we should provide _Float128 functions directly in
> glibc's libm on x86/x86_64, with the TS 18661-3 names, in which case
https://gcc.gnu.org/bugzilla/show_bug.cgi?id=63488
--- Comment #3 from joseph at codesourcery dot com ---
Note that libquadmath has not been updated from glibc since November 2012.
So, while in the Bessel function case large errors are already known for
all floating-point types in glibc, in ge
https://gcc.gnu.org/bugzilla/show_bug.cgi?id=63488
--- Comment #2 from Francois-Xavier Coudert ---
The code in question is at libquadmath/math/j0q.c, function y0q, in the branch
annotated /* 0 <= x <= 2 */
It has to do with the rational function:
/* Y0(x) = 2/pi * log(x) * J0(x) + R(x^2)
Pe
https://gcc.gnu.org/bugzilla/show_bug.cgi?id=63488
Francois-Xavier Coudert changed:
What|Removed |Added
Status|UNCONFIRMED |NEW
Last reconfirmed|
