]> git.ipfire.org Git - thirdparty/glibc.git/blame - sysdeps/ieee754/ldbl-96/w_expl.c
Use non-signaling floating-point comparisons in math functions.
[thirdparty/glibc.git] / sysdeps / ieee754 / ldbl-96 / w_expl.c
CommitLineData
92221550 1/* Copyright (C) 2011, 2012 Free Software Foundation, Inc.
bcf01e6d
UD
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
59ba27a6
PE
16 License along with the GNU C Library; if not, see
17 <http://www.gnu.org/licenses/>. */
bcf01e6d
UD
18
19#include <math.h>
20#include <math_private.h>
21
ee188d55 22static const long double
ee188d55
RM
23o_threshold= 1.135652340629414394949193107797076489134e4,
24 /* 0x400C, 0xB17217F7, 0xD1CF79AC */
25u_threshold= -1.140019167866942050398521670162263001513e4;
26 /* 0x400C, 0xB220C447, 0x69C201E8 */
27
bcf01e6d
UD
28
29/* wrapper expl */
30long double
31__expl (long double x)
ee188d55 32{
92221550 33 if (__builtin_expect (isgreater (x, o_threshold), 0))
bcf01e6d
UD
34 {
35 if (_LIB_VERSION != _IEEE_)
36 return __kernel_standard (x, x, 206);
37 }
92221550 38 else if (__builtin_expect (isless (x, u_threshold), 0))
bcf01e6d
UD
39 {
40 if (_LIB_VERSION != _IEEE_)
41 return __kernel_standard (x, x, 207);
42 }
43
44 return __ieee754_expl (x);
ee188d55 45}
b7805d0b 46hidden_def (__expl)
ee188d55 47weak_alias (__expl, expl)