1 /* @(#)e_atanh.c 5.1 93/09/24 */
3 * ====================================================
4 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
6 * Developed at SunPro, a Sun Microsystems, Inc. business.
7 * Permission to use, copy, modify, and distribute this
8 * software is freely granted, provided that this notice
10 * ====================================================
15 * 1.Reduced x to positive by atanh(-x) = -atanh(x)
18 * atanh(x) = --- * log(1 + -------) = 0.5 * log1p(2 * --------)
22 * atanh(x) = 0.5*log1p(2x+2x*x/(1-x))
25 * atanh(x) is NaN if |x| > 1 with signal;
26 * atanh(NaN) is that NaN with no signal;
27 * atanh(+-1) is +-INF with signal.
34 #include "math_private.h"
36 static const double one
= 1.0, huge
= 1e300
;
37 static const double zero
= 0.0;
45 EXTRACT_WORDS(hx
,lx
,x
);
47 if ((ix
|((lx
|(-lx
))>>31))>0x3ff00000) /* |x|>1 */
51 if(ix
<0x3e300000&&(huge
+x
)>zero
) return x
; /* x<2**-28 */
53 if(ix
<0x3fe00000) { /* x < 0.5 */
55 t
= 0.5*log1p(t
+t
*x
/(one
-x
));
57 t
= 0.5*log1p((x
+x
)/(one
-x
));
58 if(hx
>=0) return t
; else return -t
;
61 #if LDBL_MANT_DIG == DBL_MANT_DIG
62 __strong_alias(atanhl
, atanh
);
63 #endif /* LDBL_MANT_DIG == DBL_MANT_DIG */