mathinline.h 3.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108
  1. /* Inline math functions for powerpc.
  2. Copyright (C) 1995, 1996, 1997, 1998, 1999, 2000, 2004, 2006
  3. Free Software Foundation, Inc.
  4. This file is part of the GNU C Library.
  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. The GNU C Library is distributed in the hope that it will be useful,
  10. but WITHOUT ANY WARRANTY; without even the implied warranty of
  11. MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
  12. Lesser General Public License for more details.
  13. You should have received a copy of the GNU Lesser General Public
  14. License along with the GNU C Library; if not, see
  15. <http://www.gnu.org/licenses/>. */
  16. #include <features.h>
  17. #ifndef _MATH_H
  18. # error "Never use <bits/mathinline.h> directly; include <math.h> instead."
  19. #endif
  20. #ifdef __cplusplus
  21. # define __MATH_INLINE __inline
  22. #else
  23. # define __MATH_INLINE __extern_inline
  24. #endif /* __cplusplus */
  25. #if defined __GNUC__ && !defined _SOFT_FLOAT
  26. #ifdef __USE_ISOC99
  27. # if !__GNUC_PREREQ (2,97)
  28. # define __unordered_cmp(x, y) \
  29. (__extension__ \
  30. ({ __typeof__(x) __x = (x); __typeof__(y) __y = (y); \
  31. unsigned __r; \
  32. __asm__("fcmpu 7,%1,%2 ; mfcr %0" : "=r" (__r) : "f" (__x), "f"(__y) \
  33. : "cr7"); \
  34. __r; }))
  35. # undef isgreater
  36. # undef isgreaterequal
  37. # undef isless
  38. # undef islessequal
  39. # undef islessgreater
  40. # undef isunordered
  41. # define isgreater(x, y) (__unordered_cmp (x, y) >> 2 & 1)
  42. # define isgreaterequal(x, y) ((__unordered_cmp (x, y) & 6) != 0)
  43. # define isless(x, y) (__unordered_cmp (x, y) >> 3 & 1)
  44. # define islessequal(x, y) ((__unordered_cmp (x, y) & 0xA) != 0)
  45. # define islessgreater(x, y) ((__unordered_cmp (x, y) & 0xC) != 0)
  46. # define isunordered(x, y) (__unordered_cmp (x, y) & 1)
  47. # endif /* __GNUC_PREREQ (2,97) */
  48. #endif /* __USE_ISOC99 */
  49. #if !defined __NO_MATH_INLINES && defined __OPTIMIZE__
  50. #ifdef __USE_ISOC99
  51. __MATH_INLINE long int lrint (double __x) __THROW;
  52. __MATH_INLINE long int
  53. __NTH (lrint (double __x))
  54. {
  55. union {
  56. double __d;
  57. int __ll[2];
  58. } __u;
  59. __asm__ ("fctiw %0,%1" : "=f"(__u.__d) : "f"(__x));
  60. return __u.__ll[1];
  61. }
  62. __MATH_INLINE long int lrintf (float __x) __THROW;
  63. __MATH_INLINE long int
  64. __NTH (lrintf (float __x))
  65. {
  66. union {
  67. double __d;
  68. int __ll[2];
  69. } __u;
  70. __asm__ ("fctiw %0,%1" : "=f"(__u.__d) : "f"(__x));
  71. return __u.__ll[1];
  72. }
  73. __MATH_INLINE double fdim (double __x, double __y) __THROW;
  74. __MATH_INLINE double
  75. __NTH (fdim (double __x, double __y))
  76. {
  77. return __x <= __y ? 0 : __x - __y;
  78. }
  79. __MATH_INLINE float fdimf (float __x, float __y) __THROW;
  80. __MATH_INLINE float
  81. __NTH (fdimf (float __x, float __y))
  82. {
  83. return __x <= __y ? 0 : __x - __y;
  84. }
  85. #endif /* __USE_ISOC99 */
  86. #endif /* !__NO_MATH_INLINES && __OPTIMIZE__ */
  87. #endif /* __GNUC__ && !_SOFT_FLOAT */