fix threshold constants in j0f, y0f, j1f, y1f
authorSzabolcs Nagy <nsz@port70.net>
Wed, 15 Mar 2017 01:55:49 +0000 (02:55 +0100)
committerRich Felker <dalias@aerifal.cx>
Wed, 15 Mar 2017 04:15:28 +0000 (00:15 -0400)
commit8cba1dc46c8f29261aa441e70bac798f2c2c0f58
tree69ab34ac1a1dfa914a9ffdc3498592a667850309
parentcb525397bb053ea49cf160965477a17b17286eb3
fix threshold constants in j0f, y0f, j1f, y1f

partly following freebsd rev 279491
https://svnweb.freebsd.org/base?view=revision&revision=279491
(musl had some of the fixes before freebsd).

the change should not matter much for j0f, y0f, but it improves
j1f and y1f in [2.5,~3.75] (that is [0x40200000,~0x40700000]).
near roots (e.g. around 3.8317 for j1f) there are still large
ulp errors.

dropped code that tried to raise inexact.
src/math/j0f.c
src/math/j1f.c