sysdeps/s390/fpu/bits/mathinline.h

   1 /* Inline math functions for s390.
   2    Copyright (C) 2004, 2007 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 #ifndef __extern_inline
  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 # ifndef __NO_LONG_DOUBLE_MATH
  51 __MATH_INLINE int
  52 __NTH (__signbitl (long double __x))
  53 {
  54   __extension__ union { long double __l; int __i[4]; } __u = { __l: __x };
  55   return __u.__i[0] < 0;
  56 }
  57 # else
  58 __MATH_INLINE int
  59 __NTH (__signbitl (long double __x))
  60 {
  61   return __signbit ((double) __x);
  62 }
  63 # endif
  64
  65 #endif /* C99 */
  66
  67 /* This code is used internally in the GNU libc.  */
  68 #ifdef __LIBC_INTERNAL_MATH_INLINES
  69
  70 __MATH_INLINE double
  71 __NTH (__ieee754_sqrt (double x))
  72 {
  73   double res;
  74
  75   asm ( "sqdbr %0,%1" : "=f" (res) : "f" (x) );
  76   return res;
  77 }
  78
  79 __MATH_INLINE float
  80 __NTH (__ieee754_sqrtf (float x))
  81 {
  82   float res;
  83
  84   asm ( "sqebr %0,%1" : "=f" (res) : "f" (x) );
  85   return res;
  86 }
  87
  88 # if !defined __NO_LONG_DOUBLE_MATH
  89 __MATH_INLINE long double
  90 __NTH (sqrtl (long double __x))
  91 {
  92   long double res;
  93
  94   asm ( "sqxbr %0,%1" : "=f" (res) : "f" (__x) );
  95   return res;
  96 }
  97 # endif /* !__NO_LONG_DOUBLE_MATH */
  98
  99 #endif /* __LIBC_INTERNAL_MATH_INLINES */
 100
 101 #endif /* __NO_MATH_INLINES */