1 /* s_asinhl.c -- long double version of s_asinh.c.
2 * Conversion to long double by Ulrich Drepper,
3 * Cygnus Support, drepper@cygnus.com.
7 * ====================================================
8 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
10 * Developed at SunPro, a Sun Microsystems, Inc. business.
11 * Permission to use, copy, modify, and distribute this
12 * software is freely granted, provided that this notice
14 * ====================================================
17 #if defined(LIBM_SCCS) && !defined(lint)
18 static char rcsid
[] = "$NetBSD: $";
24 * asinhl(x) = signl(x) * logl [ |x| + sqrtl(x*x+1) ]
26 * asinhl(x) := x if 1+x*x=1,
27 * := signl(x)*(logl(x)+ln2)) for large |x|, else
28 * := signl(x)*logl(2|x|+1/(|x|+sqrtl(x*x+1))) if|x|>2, else
29 * := signl(x)*log1pl(|x| + x^2/(1 + sqrtl(1+x^2)))
33 #include "math_private.h"
36 static const long double
41 ln2
= 6.931471805599453094172321214581765681e-1L,
46 __asinhl (long double x
)
55 ieee854_long_double_shape_type u
;
59 ix
= sign
& 0x7fffffff;
61 return x
+ x
; /* x is inf or NaN */
65 return x
; /* return x inexact except 0 */
70 w
= __ieee754_logl (u
.value
) + ln2
;
72 else if (ix
>0x40000000)
73 { /* 2^ 54 > |x| > 2.0 */
75 w
= __ieee754_logl (2.0 * t
+ one
/ (__ieee754_sqrtl (x
* x
+ one
) + t
));
78 { /* 2.0 > |x| > 2 ^ -56 */
80 w
= __log1pl (u
.value
+ t
/ (one
+ __ieee754_sqrtl (one
+ t
)));
82 if (sign
& 0x80000000)
87 weak_alias (__asinhl
, asinhl
)