dcd506d49e61d77d7eb2867259a894495344739c
[kopensolaris-gnu/glibc.git] / sysdeps / powerpc / fpu / bits / mathinline.h
1 /* Inline math functions for powerpc.
2    Copyright (C) 1995, 1996, 1997, 1998, 1999, 2000, 2004
3    Free Software Foundation, Inc.
4    This file is part of the GNU C Library.
5
6    The GNU C Library is free software; you can redistribute it and/or
7    modify it under the terms of the GNU Lesser General Public
8    License as published by the Free Software Foundation; either
9    version 2.1 of the License, or (at your option) any later version.
10
11    The GNU C Library is distributed in the hope that it will be useful,
12    but WITHOUT ANY WARRANTY; without even the implied warranty of
13    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14    Lesser General Public License for more details.
15
16    You should have received a copy of the GNU Lesser General Public
17    License along with the GNU C Library; if not, write to the Free
18    Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
19    02111-1307 USA.  */
20
21 #if defined __GNUC__ && !defined _SOFT_FLOAT
22
23 #ifdef __USE_ISOC99
24 # if __GNUC_PREREQ (2,96)
25
26 #  define isgreater(x, y) __builtin_isgreater (x, y)
27 #  define isgreaterequal(x, y) __builtin_isgreaterequal (x, y)
28 #  define isless(x, y) __builtin_isless (x, y)
29 #  define islessequal(x, y) __builtin_islessequal (x, y)
30 #  define islessgreater(x, y) __builtin_islessgreater (x, y)
31 #  define isunordered(x, y) __builtin_isunordered (x, y)
32
33 # else
34
35 #  define __unordered_cmp(x, y) \
36   (__extension__                                                              \
37    ({ __typeof__(x) __x = (x); __typeof__(y) __y = (y);                       \
38       unsigned __r;                                                           \
39       __asm__("fcmpu 7,%1,%2 ; mfcr %0" : "=r" (__r) : "f" (__x), "f"(__y)    \
40               : "cr7");  \
41       __r; }))
42
43 #  define isgreater(x, y) (__unordered_cmp (x, y) >> 2 & 1)
44 #  define isgreaterequal(x, y) ((__unordered_cmp (x, y) & 6) != 0)
45 #  define isless(x, y) (__unordered_cmp (x, y) >> 3 & 1)
46 #  define islessequal(x, y) ((__unordered_cmp (x, y) & 0xA) != 0)
47 #  define islessgreater(x, y) ((__unordered_cmp (x, y) & 0xC) != 0)
48 #  define isunordered(x, y) (__unordered_cmp (x, y) & 1)
49
50 # endif /* __GNUC_PREREQ (2,97) */
51 #endif /* __USE_ISOC99 */
52
53 #if !defined __NO_MATH_INLINES && defined __OPTIMIZE__
54
55 #ifdef __cplusplus
56 # define __MATH_INLINE __inline
57 #else
58 # define __MATH_INLINE extern __inline
59 #endif  /* __cplusplus */
60
61 #ifdef __USE_ISOC99
62
63 # ifndef __powerpc64__
64 __MATH_INLINE long int lrint (double __x) __THROW;
65 __MATH_INLINE long int
66 lrint (double __x) __THROW
67 {
68   union {
69     double __d;   
70     int __ll[2];
71   } __u;
72   __asm__ ("fctiw %0,%1" : "=f"(__u.__d) : "f"(__x));
73   return __u.__ll[1];
74 }
75
76 __MATH_INLINE long int lrintf (float __x) __THROW;
77 __MATH_INLINE long int
78 lrintf (float __x) __THROW
79 {
80   union {
81     double __d;
82     int __ll[2];
83   } __u;
84   __asm__ ("fctiw %0,%1" : "=f"(__u.__d) : "f"(__x));
85   return __u.__ll[1];
86 }
87 # endif
88
89 __MATH_INLINE double fdim (double __x, double __y) __THROW;
90 __MATH_INLINE double
91 fdim (double __x, double __y) __THROW
92 {
93   return __x < __y ? 0 : __x - __y;
94 }
95
96 __MATH_INLINE float fdimf (float __x, float __y) __THROW;
97 __MATH_INLINE float
98 fdimf (float __x, float __y) __THROW
99 {
100   return __x < __y ? 0 : __x - __y;
101 }
102
103 #endif /* __USE_ISOC99 */
104 #endif /* !__NO_MATH_INLINES && __OPTIMIZE__ */
105 #endif /* __GNUC__ && !_SOFT_FLOAT */