1 /* Inline math functions for powerpc.
2 Copyright (C) 1995-2013 Free Software Foundation, Inc.
3 This file is part of the GNU C Library.
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.
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.
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 <http://www.gnu.org/licenses/>. */
20 # error "Never use <bits/mathinline.h> directly; include <math.h> instead."
23 #ifndef __extern_inline
24 # define __MATH_INLINE __inline
26 # define __MATH_INLINE __extern_inline
27 #endif /* __cplusplus */
29 #if defined __GNUC__ && !defined _SOFT_FLOAT
32 # if !__GNUC_PREREQ (2,97)
33 # define __unordered_cmp(x, y) \
35 ({ __typeof__(x) __x = (x); __typeof__(y) __y = (y); \
37 __asm__("fcmpu 7,%1,%2 ; mfcr %0" : "=r" (__r) : "f" (__x), "f"(__y) \
42 # undef isgreaterequal
48 # define isgreater(x, y) (__unordered_cmp (x, y) >> 2 & 1)
49 # define isgreaterequal(x, y) ((__unordered_cmp (x, y) & 6) != 0)
50 # define isless(x, y) (__unordered_cmp (x, y) >> 3 & 1)
51 # define islessequal(x, y) ((__unordered_cmp (x, y) & 0xA) != 0)
52 # define islessgreater(x, y) ((__unordered_cmp (x, y) & 0xC) != 0)
53 # define isunordered(x, y) (__unordered_cmp (x, y) & 1)
55 # endif /* __GNUC_PREREQ (2,97) */
57 /* The gcc, version 2.7 or below, has problems with all this inlining
58 code. So disable it for this version of the compiler. */
59 # if __GNUC_PREREQ (2, 8)
60 /* Test for negative number. Used in the signbit() macro. */
62 __NTH (__signbitf (float __x
))
64 __extension__
union { float __f
; int __i
; } __u
= { __f
: __x
};
68 __NTH (__signbit (double __x
))
70 __extension__
union { double __d
; int __i
[2]; } __u
= { __d
: __x
};
71 return __u
.__i
[0] < 0;
73 # ifdef __LONG_DOUBLE_128__
75 __NTH (__signbitl (long double __x
))
77 __extension__
union { long double __d
; int __i
[4]; } __u
= { __d
: __x
};
78 return __u
.__i
[0] < 0;
82 #endif /* __USE_ISOC99 */
84 #if !defined __NO_MATH_INLINES && defined __OPTIMIZE__
88 # ifndef __powerpc64__
89 __MATH_INLINE
long int lrint (double __x
) __THROW
;
90 __MATH_INLINE
long int
91 __NTH (lrint (double __x
))
97 __asm__ ("fctiw %0,%1" : "=f"(__u
.__d
) : "f"(__x
));
101 __MATH_INLINE
long int lrintf (float __x
) __THROW
;
102 __MATH_INLINE
long int
103 __NTH (lrintf (float __x
))
109 __asm__ ("fctiw %0,%1" : "=f"(__u
.__d
) : "f"(__x
));
114 __MATH_INLINE
double fdim (double __x
, double __y
) __THROW
;
116 __NTH (fdim (double __x
, double __y
))
118 return __x
<= __y
? 0 : __x
- __y
;
121 __MATH_INLINE
float fdimf (float __x
, float __y
) __THROW
;
123 __NTH (fdimf (float __x
, float __y
))
125 return __x
<= __y
? 0 : __x
- __y
;
128 #endif /* __USE_ISOC99 */
129 #endif /* !__NO_MATH_INLINES && __OPTIMIZE__ */
130 #endif /* __GNUC__ && !_SOFT_FLOAT */