X-Git-Url: http://erislabs.net/gitweb/?a=blobdiff_plain;ds=sidebyside;f=lib%2Fatanl.c;h=88284a67195a329847da18f2f75a57d1a12fd30b;hb=e2811532647e8280f13347d71579621ac441c0cd;hp=eca41774d048d369b812229b501a1a44b066c223;hpb=57fdfd3f8ec62b105c53bcdf6f127c35c7fe7391;p=gnulib.git diff --git a/lib/atanl.c b/lib/atanl.c index eca41774d..88284a671 100644 --- a/lib/atanl.c +++ b/lib/atanl.c @@ -18,9 +18,21 @@ /* Specification. */ #include -/* s_atanl.c +#if HAVE_SAME_LONG_DOUBLE_AS_DOUBLE + +long double +atanl (long double x) +{ + return atan (x); +} + +#else + +/* Code based on glibc/sysdeps/ieee754/ldbl-128/s_atanl.c. */ + +/* s_atanl.c * - * Inverse circular tangent for 128-bit long double precision + * Inverse circular tangent for 128-bit long double precision * (arctangent) * * @@ -62,8 +74,6 @@ * */ -#include "isnanl.h" - /* arctan(k/8), k = 0, ..., 82 */ static const long double atantbl[84] = { 0.0000000000000000000000000000000000000000E0L, @@ -187,9 +197,9 @@ atanl (long double x) { /* Infinity. */ if (sign) - return -atantbl[83]; + return -atantbl[83]; else - return atantbl[83]; + return atantbl[83]; } if (sign) @@ -203,7 +213,7 @@ atanl (long double x) else { /* Index of nearest table element. - Roundoff to integer is asymmetrical to avoid cancellation when t < 0 + Roundoff to integer is asymmetrical to avoid cancellation when t < 0 (cf. fdlibm). */ k = 8.0 * x + 0.25; u = 0.125 * k; @@ -224,3 +234,5 @@ atanl (long double x) else return u; } + +#endif