1 /* w_expl.c -- long double version of w_exp.c.
2 */
3
4 /*
5 * ====================================================
6 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
7 *
8 * Developed at SunPro, a Sun Microsystems, Inc. business.
9 * Permission to use, copy, modify, and distribute this
10 * software is freely granted, provided that this notice
11 * is preserved.
12 * ====================================================
13 */
14
15 #if defined(LIBM_SCCS) && !defined(lint)
16 static char rcsid[] = "$NetBSD: $";
17 #endif
18
19 /*
20 * wrapper expl(x)
21 */
22
23 #include <math.h>
24 #include <math_private.h>
25 #include <math-svid-compat.h>
26 #include <libm-alias-ldouble.h>
27
28 #if LIBM_SVID_COMPAT
__expl(long double x)29 long double __expl(long double x) /* wrapper exp */
30 {
31 # ifdef _IEEE_LIBM
32 return __ieee754_expl(x);
33 # else
34 long double z = __ieee754_expl (x);
35 if (__glibc_unlikely (!isfinite (z) || z == 0)
36 && isfinite (x) && _LIB_VERSION != _IEEE_)
37 return __kernel_standard_l (x, x, 206 + !!signbit (x));
38
39 return z;
40 # endif
41 }
42 libm_hidden_def (__expl)
43 libm_alias_ldouble (__exp, exp)
44 #endif
45