[PATCH] PPC fixes for long double

[Steven Munroe]

This patch resolves the ulps failures in cprojl, ctanhl, ctanl, ceill,
floorl, and the Bessel functions (j0, j1, jn, y0, y1, yn).

The double double runtime does not propigate -0.0 correctly in some
case. So I have to add fixup code to some complex functions. Similarly
ceill, floorl needed explicite handling for the correct sign of 0.0;

Also found that the k_cosl, k_sinl, k_sincosl did not compute the index
into the __sincosl_table correcly for double double. The result is bad
values for larger radian inputs and cause the Bessel function failures.
Hacked around it by converting the high 32-bit (integer form) of the
high double to ieee extended and used the index calc from the ldbl-128
(ieee extended) version.


-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: ppc-ldsjm-20060130.txt
URL: <http://sourceware.org/pipermail/libc-alpha/attachments/20060131/fc0c82cb/attachment.txt>