PowerPC specific math inline functions.
[kopensolaris-gnu/glibc.git] / sysdeps / powerpc / fpu / bits / mathinline.h
1 /* Inline math functions for powerpc.
2    Copyright (C) 1995, 1996, 1997, 1998, 1999 Free Software Foundation, Inc.
3    This file is part of the GNU C Library.
4
5    The GNU C Library is free software; you can redistribute it and/or
6    modify it under the terms of the GNU Library General Public License as
7    published by the Free Software Foundation; either version 2 of the
8    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    Library General Public License for more details.
14
15    You should have received a copy of the GNU Library General Public
16    License along with the GNU C Library; see the file COPYING.LIB.  If not,
17    write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
18    Boston, MA 02111-1307, USA.  */
19
20 #if defined __GNUC__ && !defined _SOFT_FLOAT
21
22 #ifdef __USE_ISOC99
23 # define __unordered_cmp(x, y) \
24   (__extension__                                                              \
25    ({ __typeof__(x) __x = (x); __typeof__(y) __y = (y);                       \
26       unsigned __r;                                                           \
27       __asm__("fcmpu 7,%1,%2 ; mfcr %0" : "=r" (__r) : "f" (__x), "f"(__y)    \
28               : "cr7");  \
29       __r; }))
30
31 # define isgreater(x, y) (__unordered_cmp (x, y) >> 2 & 1)
32 # define isgreaterequal(x, y) ((__unordered_cmp (x, y) & 6) != 0)
33 # define isless(x, y) (__unordered_cmp (x, y) >> 3 & 1)
34 # define islessequal(x, y) ((__unordered_cmp (x, y) & 0xA) != 0)
35 # define islessgreater(x, y) ((__unordered_cmp (x, y) & 0xC) != 0)
36 # define isunordered(x, y) (__unordered_cmp (x, y) & 1)
37 #endif /* __USE_ISOC99 && !_SOFT_FLOAT */
38
39 #if !defined __NO_MATH_INLINES && defined __OPTIMIZE__
40
41 #ifdef __cplusplus
42 # define __MATH_INLINE __inline
43 #else
44 # define __MATH_INLINE extern __inline
45 #endif  /* __cplusplus */
46
47 #ifdef __USE_ISOC99
48 __MATH_INLINE long int lrint (double __x);
49 __MATH_INLINE long int
50 lrint (double __x)
51 {
52   union {
53     double __d;
54     long int __ll[2];
55   } __u;
56   __asm__ ("fctiw %0,%1" : "=f"(__u.__d) : "f"(__x));
57   return __u.__ll[1];
58 }
59
60 __MATH_INLINE long int lrintf (float __x);
61 __MATH_INLINE long int
62 lrintf (float __x)
63 {
64   union {
65     double __d;
66     long int __ll[2];
67   } __u;
68   __asm__ ("fctiw %0,%1" : "=f"(__u.__d) : "f"(__x));
69   return __u.__ll[1];
70 }
71
72 __MATH_INLINE double fdim (double __x, double __y);
73 __MATH_INLINE double
74 fdim (double __x, double __y)
75 {
76   return __x < __y ? 0 : __x - __y;
77 }
78
79 __MATH_INLINE float fdimf (float __x, float __y);
80 __MATH_INLINE float
81 fdimf (float __x, float __y)
82 {
83   return __x < __y ? 0 : __x - __y;
84 }
85
86 #endif /* __USE_ISOC99 */
87 #endif /* !__NO_MATH_INLINES && __OPTIMIZE__ */
88 #endif /* __GNUC__ && !_SOFT_FLOAT */