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

   1 /* origin: FreeBSD /usr/src/lib/msun/src/s_asinhf.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
  19 ln2 = 6.9314718246e-01, /* 0x3f317218 */
  20 huge= 1.0000000000e+30;
  21
  22 float asinhf(float x)
  23 {
  24         float t,w;
  25         int32_t hx,ix;
  26
  27         GET_FLOAT_WORD(hx, x);
  28         ix = hx & 0x7fffffff;
  29         if (ix >= 0x7f800000)   /* x is inf or NaN */
  30                 return x+x;
  31         if (ix < 0x31800000) {  /* |x| < 2**-28 */
  32                 /* return x inexact except 0 */
  33                 if (huge+x > 1.0f)
  34                         return x;
  35         }
  36         if (ix > 0x4d800000) {  /* |x| > 2**28 */
  37                 w = logf(fabsf(x)) + ln2;
  38         } else if (ix > 0x40000000) {  /* 2**28 > |x| > 2.0 */
  39                 t = fabsf(x);
  40                 w = logf(2.0f*t + 1.0f/(sqrtf(x*x+1.0f)+t));
  41         } else {                /* 2.0 > |x| > 2**-28 */
  42                 t = x*x;
  43                 w =log1pf(fabsf(x) + t/(1.0f+sqrtf(1.0f+t)));
  44         }
  45         if (hx > 0)
  46                 return w;
  47         return -w;
  48 }