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 */