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

   1 /* origin: FreeBSD /usr/src/lib/msun/src/e_atanhf.c */
   2 /*
   3  * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
   4  */
   5 /*
   6  * ====================================================
   7  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
   8  *
   9  * Developed at SunPro, a Sun Microsystems, Inc. business.
  10  * Permission to use, copy, modify, and distribute this
  11  * software is freely granted, provided that this notice
  12  * is preserved.
  13  * ====================================================
  14  */
  15
  16 #include "libm.h"
  17
  18 static const float huge = 1e30;
  19
  20 float atanhf(float x)
  21 {
  22         float t;
  23         int32_t hx,ix;
  24
  25         GET_FLOAT_WORD(hx, x);
  26         ix = hx & 0x7fffffff;
  27         if (ix > 0x3f800000)                   /* |x| > 1 */
  28                 return (x-x)/(x-x);
  29         if (ix == 0x3f800000)
  30                 return x/0.0f;
  31         if (ix < 0x31800000 && huge+x > 0.0f)  /* x < 2**-28 */
  32                 return x;
  33         SET_FLOAT_WORD(x, ix);
  34         if (ix < 0x3f000000) {                 /* x < 0.5 */
  35                 t = x+x;
  36                 t = 0.5f*log1pf(t + t*x/(1.0f-x));
  37         } else
  38                 t = 0.5f*log1pf((x+x)/(1.0f-x));
  39         if (hx >= 0)
  40                 return t;
  41         return -t;
  42 }