]> git.ipfire.org Git - thirdparty/glibc.git/blob - math/w_exp10l_compat.c
Improve the accuracy of tgamma (BZ #26983)
[thirdparty/glibc.git] / math / w_exp10l_compat.c
1 /* Copyright (C) 2011-2021 Free Software Foundation, Inc.
2 This file is part of the GNU C Library.
3 Contributed by Ulrich Drepper <drepper@gmail.com>, 2011.
4
5 The GNU C Library is free software; you can redistribute it and/or
6 modify it under the terms of the GNU Lesser General Public
7 License as published by the Free Software Foundation; either
8 version 2.1 of the License, or (at your option) any later version.
9
10 The GNU C Library is distributed in the hope that it will be useful,
11 but WITHOUT ANY WARRANTY; without even the implied warranty of
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 Lesser General Public License for more details.
14
15 You should have received a copy of the GNU Lesser General Public
16 License along with the GNU C Library; if not, see
17 <https://www.gnu.org/licenses/>. */
18
19
20 /*
21 * wrapper exp10l(x)
22 */
23
24 #include <math.h>
25 #include <math_private.h>
26 #include <math-svid-compat.h>
27 #include <libm-alias-ldouble.h>
28
29 #if LIBM_SVID_COMPAT
30 long double
31 __exp10l (long double x)
32 {
33 long double z = __ieee754_exp10l (x);
34 if (__builtin_expect (!isfinite (z) || z == 0, 0)
35 && isfinite (x) && _LIB_VERSION != _IEEE_)
36 /* exp10l overflow (246) if x > 0, underflow (247) if x < 0. */
37 return __kernel_standard_l (x, x, 246 + !!signbit (x));
38
39 return z;
40 }
41 libm_alias_ldouble (__exp10, exp10)
42 # if SHLIB_COMPAT (libm, GLIBC_2_1, GLIBC_2_27)
43 strong_alias (__exp10l, __pow10l)
44 compat_symbol (libm, __pow10l, pow10l, GLIBC_2_1);
45 # endif
46 #endif