numeric slightly revised fdlibm 5.3

David G Hough at validlab.com validlabavalidlab.com
Thu May 6 18:28:00 PDT 2004


A further slight revision has been place at
http://www.validlab.com/software/fdlibm53.tar.gz

The code is the same, but a configure script has been contributed by
Nelson Beebe.

The previous changes to pow and tan are described in the README:

What's new in FDLIBM 5.3?

CONFIGURE
        To build FDLIBM, edit the supplied Makefile or create
        a local Makefile by running "sh configure"
        using the supplied configure script contributed by Nelson Beebe

BUGS FIXED

    1. e_pow.c incorrect results when
        x is very close to -1.0 and y is very large, e.g.
        pow(-1.0000000000000002e+00,4.5035996273704970e+15) = 0
        pow(-9.9999999999999978e-01,4.5035996273704970e+15) = 0
        Correct results are close to -e and -1/e.

    2. k_tan.c error was > 1 ulp target for FDLIBM
        5.2: Worst error at least 1.45 ulp at
        tan(1.7765241907548024E+269) = 1.7733884462610958E+16
        5.3: Worst error 0.96 ulp

NOT FIXED YET

    3. Compiler failure on non-standard code
        Statements like
                    *(1+(int*)&t1) = 0;
        are not standard C and cause some optimizing compilers (e.g. GCC)
        to generate bad code under optimization.    These cases
        are to be addressed in the next release.
       




More information about the Numeric-interest mailing list