sysdeps/powerpc/bits/fenvinline.h

   1 /* Inline floating-point environment handling functions for powerpc.
   2    Copyright (C) 1995-2020 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    <https://www.gnu.org/licenses/>.  */
  18
  19 #if defined __GNUC__ && !defined _SOFT_FLOAT && !defined __NO_FPRS__
  20
  21 /* Inline definitions for fegetround.  */
  22 # define __fegetround_ISA300()                                          \
  23   (__extension__  ({                                                    \
  24     union { double __d; unsigned long long __ll; } __u;                 \
  25     __asm__ __volatile__ (                                              \
  26       ".machine push; .machine \"power9\"; mffsl %0; .machine pop"      \
  27       : "=f" (__u.__d));                                                \
  28     __u.__ll & 0x0000000000000003LL;                                    \
  29   }))
  30
  31 # define __fegetround_ISA2()                                            \
  32   (__extension__  ({                                                    \
  33      int __fegetround_result;                                           \
  34      __asm__ __volatile__ ("mcrfs 7,7 ; mfcr %0"                        \
  35                            : "=r"(__fegetround_result) : : "cr7");      \
  36      __fegetround_result & 3;                                           \
  37   }))
  38
  39 # ifdef _ARCH_PWR9
  40 #  define __fegetround() __fegetround_ISA300()
  41 # elif defined __BUILTIN_CPU_SUPPORTS__
  42 #  define __fegetround()                                                \
  43   (__glibc_likely (__builtin_cpu_supports ("arch_3_00"))                \
  44    ? __fegetround_ISA300()                                              \
  45    : __fegetround_ISA2()                                                \
  46   )
  47 # else
  48 #  define __fegetround() __fegetround_ISA2()
  49 # endif
  50
  51 # define fegetround() __fegetround ()
  52
  53 # ifndef __NO_MATH_INLINES
  54
  55 /* Builtins to mtfsb0 and mtfsb1 was introduced on GCC 9.  */
  56 #  if !__GNUC_PREREQ(9, 0)
  57 /* The weird 'i#*X' constraints on the following suppress a gcc
  58    warning when __excepts is not a constant.  Otherwise, they mean the
  59    same as just plain 'i'.  This warning only happens in old GCC
  60    versions (gcc 3 or less).  Otherwise plain 'i' works fine.  */
  61 #   define __MTFSB0(__b) __asm__ __volatile__ ("mtfsb0 %0" : : "i#*X" (__b))
  62 #   define __MTFSB1(__b) __asm__ __volatile__ ("mtfsb1 %0" : : "i#*X" (__b))
  63 #  else
  64 #   define __MTFSB0(__b) __builtin_mtfsb0 (__b)
  65 #   define __MTFSB1(__b) __builtin_mtfsb1 (__b)
  66 #  endif
  67
  68 #  if __GNUC_PREREQ(3, 4)
  69
  70 #include <sys/param.h>
  71
  72 /* Inline definition for feraiseexcept.  */
  73 #   define feraiseexcept(__excepts) \
  74   (__extension__  ({                                                          \
  75     int __e = __excepts;                                                      \
  76     int __ret = 0;                                                            \
  77     if (__builtin_constant_p (__e)                                            \
  78         && __builtin_popcount (__e) == 1                                      \
  79         && __e != FE_INVALID)                                                 \
  80       {                                                                       \
  81         __MTFSB1 ((__builtin_clz (__e)));                                     \
  82       }                                                                       \
  83     else                                                                      \
  84       __ret = feraiseexcept (__e);                                            \
  85     __ret;                                                                    \
  86   }))
  87
  88 /* Inline definition for feclearexcept.  */
  89 #   define feclearexcept(__excepts) \
  90   (__extension__  ({                                                          \
  91     int __e = __excepts;                                                      \
  92     int __ret = 0;                                                            \
  93     if (__builtin_constant_p (__e)                                            \
  94         && __builtin_popcount (__e) == 1                                      \
  95         && __e != FE_INVALID)                                                 \
  96       {                                                                       \
  97         __MTFSB0 ((__builtin_clz (__e)));                                     \
  98       }                                                                       \
  99     else                                                                      \
 100       __ret = feclearexcept (__e);                                            \
 101     __ret;                                                                    \
 102   }))
 103
 104 #  endif /* __GNUC_PREREQ(3, 4).  */
 105
 106 # endif /* !__NO_MATH_INLINES.  */
 107
 108 #endif /* __GNUC__ && !_SOFT_FLOAT && !__NO_FPRS__ */