Update ENTERI() macro
Bruce Evans
brde at optusnet.com.au
Thu Mar 7 04:09:14 UTC 2019
On Wed, 6 Mar 2019, Steve Kargl wrote:
> On Thu, Mar 07, 2019 at 06:30:42AM +1100, Bruce Evans wrote:
>> ...
>> I now see that you implemented 2 more versions of __ldexpl_cexpl() by
>> cloning the old double precision version. Apparently the includes
>> are to unpolluted for the compiler to see the multiple versions :-).
>>
>> Using the version in k_expl.h almost forces inlining of expl()'s kernel
>> and its large tables, just like for hyperbolic functions. This wastes
>> a lot of space, especially for duplicating the tables. It is only a
>> small optimization for time. It is done for the hyperbolic functions
>> to get this optimization, and for __ldexpl_cexpl() just for convenience.
>
> The version in k_expl.h has 2 bugs. You note the first (cos instead
> of cosl). The second is
>
> In file included from /data/kargl/trunk/math/libm/msun/ld80/s_cexpl.c:43:
> /data/kargl/trunk/math/libm/msun/ld80/k_expl.h:288:22: error: magnitude of
> floating-point constant too large for type 'double'; maximum is
> 1.7976931348623157E+308 [-Werror,-Wliteral-range]
> exp_x = (lo + hi) * 0x1p16382;
> ^
It is missing an L suffix. Clearly not tested.
>>> [errors for cexpl()]
>>
>> Check a few denormals, Infs and NaNs.
>
> Exceptional cases give
>
> % ./testl -e
> cexpl(1, inf) = (nan,nan) expecting (nan,nan)
> cexpl(1,-inf) = (nan,nan) expecting (nan,nan)
> cexpl(nan,-inf) = (nan,nan) expecting (nan,nan)
> cexpl(nan, inf) = (nan,nan) expecting (nan,nan)
> cexpl(inf, inf) = (inf,nan) expecting (inf,nan)
> cexpl(inf,-inf) = (inf,nan) expecting (inf,nan)
> cexpl(inf,nan) = (inf,nan) expecting (inf,nan)
> cexpl(-inf,nan) = (0,0) expecting (0,0)
> cexpl(-inf,-inf) = (0,0) expecting (0,0)
> cexpl(-inf, inf) = (0,0) expecting (0,0)
But does it preserve NaN bits as much as possible, and produce the same
NaN bits as double precision after conversion to double precision? The
report spells NaN and Inf using C99's fuzzy bad names. Even the C99
specification spells NaN and Inf correctly in Appendix F.
The complex hyperbolic functions have better comments, so are almost
an easier place to start.
Bruce
More information about the freebsd-numerics
mailing list