[__LIBC_INTERNAL_MATH_INLINES]
[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 __signbitf (float __x) __THROW
37 {
38   __extension__ union { float __f; int __i; } __u = { __f: __x };
39   return __u.__i < 0;
40 }
41
42 __MATH_INLINE int __signbit (double __x) __THROW
43 {
44   __extension__ union { double __d; long __i; } __u = { __d: __x };
45   return __u.__i < 0;
46 }
47
48 #endif /* C99 */
49
50 /* This code is used internally in the GNU libc.  */
51 #ifdef __LIBC_INTERNAL_MATH_INLINES
52
53 #include <sysdep.h>
54 #include <ldsodefs.h>
55 #include <dl-procinfo.h>
56
57 __MATH_INLINE double
58 __ieee754_sqrt (double x)
59 {
60   double res;
61      
62   asm ( "sqdbr %0,%1" : "=f" (res) : "f" (x) );
63   return res;
64 }
65
66 __MATH_INLINE float
67 __ieee754_sqrtf (float x)
68 {
69   float res;
70   
71   asm ( "sqebr %0,%1" : "=f" (res) : "f" (x) );
72   return res;
73 }
74
75 #endif /* __LIBC_INTERNAL_MATH_INLINES */
76
77 #endif /* __NO_MATH_INLINES */