sysdeps/s390/fpu/bits/mathinline.h

   1 /* Inline math functions for s390.
   2    Copyright (C) 2004-2014 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, see
  17    <http://www.gnu.org/licenses/>.  */
  18
  19 #ifndef _MATH_H
  20 # error "Never use <bits/mathinline.h> directly; include <math.h> instead."
  21 #endif
  22
  23 #ifndef __extern_inline
  24 # define __MATH_INLINE __inline
  25 #else
  26 # define __MATH_INLINE __extern_inline
  27 #endif
  28
  29 #if (!defined __NO_MATH_INLINES || defined __LIBC_INTERNAL_MATH_INLINES) \
  30     && defined __OPTIMIZE__
  31
  32 #ifdef __USE_ISOC99
  33
  34 /* Test for negative number.  Used in the signbit() macro.  */
  35 __MATH_INLINE int
  36 __NTH (__signbitf (float __x))
  37 {
  38   __extension__ union { float __f; int __i; } __u = { __f: __x };
  39   return __u.__i < 0;
  40 }
  41
  42 __MATH_INLINE int
  43 __NTH (__signbit (double __x))
  44 {
  45   __extension__ union { double __d; long __i; } __u = { __d: __x };
  46   return __u.__i < 0;
  47 }
  48
  49 # ifndef __NO_LONG_DOUBLE_MATH
  50 __MATH_INLINE int
  51 __NTH (__signbitl (long double __x))
  52 {
  53   __extension__ union { long double __l; int __i[4]; } __u = { __l: __x };
  54   return __u.__i[0] < 0;
  55 }
  56 # else
  57 __MATH_INLINE int
  58 __NTH (__signbitl (long double __x))
  59 {
  60   return __signbit ((double) __x);
  61 }
  62 # endif
  63
  64 #endif /* C99 */
  65
  66 /* This code is used internally in the GNU libc.  */
  67 #ifdef __LIBC_INTERNAL_MATH_INLINES
  68
  69 __MATH_INLINE double
  70 __NTH (__ieee754_sqrt (double x))
  71 {
  72   double res;
  73
  74   asm ( "sqdbr %0,%1" : "=f" (res) : "f" (x) );
  75   return res;
  76 }
  77
  78 __MATH_INLINE float
  79 __NTH (__ieee754_sqrtf (float x))
  80 {
  81   float res;
  82
  83   asm ( "sqebr %0,%1" : "=f" (res) : "f" (x) );
  84   return res;
  85 }
  86
  87 # if !defined __NO_LONG_DOUBLE_MATH
  88 __MATH_INLINE long double
  89 __NTH (sqrtl (long double __x))
  90 {
  91   long double res;
  92
  93   asm ( "sqxbr %0,%1" : "=f" (res) : "f" (__x) );
  94   return res;
  95 }
  96 # endif /* !__NO_LONG_DOUBLE_MATH */
  97
  98 #endif /* __LIBC_INTERNAL_MATH_INLINES */
  99
 100 #endif /* __NO_MATH_INLINES */