nsz Git - musl/blob - src/math/atanhl.c

   1 #include "libm.h"
   2
   3 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
   4 long double atanhl(long double x)
   5 {
   6         return atanh(x);
   7 }
   8 #elif LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384
   9 /* atanh(x) = log((1+x)/(1-x))/2 = log1p(2x/(1-x))/2 ~= x + x^3/3 + o(x^5) */
  10 long double atanhl(long double x)
  11 {
  12         union {
  13                 long double f;
  14                 struct{uint64_t m; uint16_t se; uint16_t pad;} i;
  15         } u = {.f = x};
  16         unsigned e = u.i.se & 0x7fff;
  17         unsigned s = u.i.se >> 15;
  18
  19         /* |x| */
  20         u.i.se = e;
  21         x = u.f;
  22
  23         if (e < 0x3fff - 1) {
  24                 /* |x| < 0.5, up to 1.7ulp error */
  25                 x = 0.5*log1pl(2*x + 2*x*x/(1-x));
  26         } else {
  27                 x = 0.5*log1pl(2*x/(1-x));
  28         }
  29         return s ? -x : x;
  30 }
  31 #endif