[J3] C TS 18661

[Steven G. Kargl]

On Wed, Jan 30, 2019 at 01:42:39PM -0700, Keith Bierman via J3 wrote:
> On Wed, Jan 30, 2019 at 1:29 PM Vipul Parekh via J3 <
> j3 at mailman.j3-fortran.org> wrote:
> 
> > On Wed, Jan 30, 2019 at 11:09 AM Steve Lionel via J3
> > <j3 at mailman.j3-fortran.org> wrote:
> > >
> > > ..  A case for the other intrinsics mentioned would
> > > have to be made using different arguments. ..
> > > ..
> > ...e BESSEL functions
> > introduced in Fortran 2008 available some place?
> >
> >
> There's a lot that goes into such things. Getting the fp details right is
> hard (fdlibm is still probably the closest portable library), and as Van
> noted, it took him 17 years to find a bug. Since Oracle seems to have
> destaffed the project (long after my departure) I'm not sure if/when anyone
> will fix *that* bug. No doubt others may emerge....
> 

Actually, Steve (that would be me) found the fdlibm jn.c bug.
FreeBSD's libm is derived from fdlibm with the float functions
coming from Cygnus Solutions port of the fdlibm's double function.
If you look through the FreeBSD source code repository for libm,
you'll find my fingerprints fixing numerous issues in the
float and double version as well as providing a number of the
missing long double versions.

openlibm, maintained by the Juila project [1], is likely now the
most portable, freely, available libm.  openlibm is derived from
FreeBSD's libm.

[1] https://github.com/JuliaMath/openlibm

-- 
Steve