Remove unnecessary includes.
[kopensolaris-gnu/glibc.git] / sysdeps / s390 / fpu / bits / mathinline.h
1 /* Inline math functions for s390.
2    Copyright (C) 2004 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 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, write to the Free
17    Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
18    02111-1307 USA.  */
19
20 #ifndef _MATH_H
21 # error "Never use <bits/mathinline.h> directly; include <math.h> instead."
22 #endif
23
24 #ifdef __cplusplus
25 # define __MATH_INLINE __inline
26 #else
27 # define __MATH_INLINE extern __inline
28 #endif
29
30 #if (!defined __NO_MATH_INLINES || defined __LIBC_INTERNAL_MATH_INLINES) \
31     && defined __OPTIMIZE__
32
33 #ifdef __USE_ISOC99
34
35 /* Test for negative number.  Used in the signbit() macro.  */
36 __MATH_INLINE int
37 __NTH (__signbitf (float __x))
38 {
39   __extension__ union { float __f; int __i; } __u = { __f: __x };
40   return __u.__i < 0;
41 }
42
43 __MATH_INLINE int
44 __NTH (__signbit (double __x))
45 {
46   __extension__ union { double __d; long __i; } __u = { __d: __x };
47   return __u.__i < 0;
48 }
49
50 #endif /* C99 */
51
52 /* This code is used internally in the GNU libc.  */
53 #ifdef __LIBC_INTERNAL_MATH_INLINES
54
55 __MATH_INLINE double
56 __NTH (__ieee754_sqrt (double x))
57 {
58   double res;
59
60   asm ( "sqdbr %0,%1" : "=f" (res) : "f" (x) );
61   return res;
62 }
63
64 __MATH_INLINE float
65 __NTH (__ieee754_sqrtf (float x))
66 {
67   float res;
68
69   asm ( "sqebr %0,%1" : "=f" (res) : "f" (x) );
70   return res;
71 }
72
73 #endif /* __LIBC_INTERNAL_MATH_INLINES */
74
75 #endif /* __NO_MATH_INLINES */