soft-fp.h
上传用户:jlfgdled
上传日期:2013-04-10
资源大小:33168k
文件大小:5k
源码类别:

Linux/Unix编程

开发平台:

Unix_Linux

  1. /* Software floating-point emulation.
  2.    Copyright (C) 1997,1998,1999 Free Software Foundation, Inc.
  3.    This file is part of the GNU C Library.
  4.    Contributed by Richard Henderson (rth@cygnus.com),
  5.   Jakub Jelinek (jj@ultra.linux.cz),
  6.   David S. Miller (davem@redhat.com) and
  7.   Peter Maydell (pmaydell@chiark.greenend.org.uk).
  8.    The GNU C Library is free software; you can redistribute it and/or
  9.    modify it under the terms of the GNU Library General Public License as
  10.    published by the Free Software Foundation; either version 2 of the
  11.    License, or (at your option) any later version.
  12.    The GNU C Library is distributed in the hope that it will be useful,
  13.    but WITHOUT ANY WARRANTY; without even the implied warranty of
  14.    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  15.    Library General Public License for more details.
  16.    You should have received a copy of the GNU Library General Public
  17.    License along with the GNU C Library; see the file COPYING.LIB.  If
  18.    not, write to the Free Software Foundation, Inc.,
  19.    59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.  */
  20. #ifndef __MATH_EMU_SOFT_FP_H__
  21. #define __MATH_EMU_SOFT_FP_H__
  22. #include <asm/sfp-machine.h>
  23. /* Allow sfp-machine to have its own byte order definitions. */
  24. #ifndef __BYTE_ORDER
  25. #include <endian.h>
  26. #endif
  27. #define _FP_WORKBITS 3
  28. #define _FP_WORK_LSB ((_FP_W_TYPE)1 << 3)
  29. #define _FP_WORK_ROUND ((_FP_W_TYPE)1 << 2)
  30. #define _FP_WORK_GUARD ((_FP_W_TYPE)1 << 1)
  31. #define _FP_WORK_STICKY ((_FP_W_TYPE)1 << 0)
  32. #ifndef FP_RND_NEAREST
  33. # define FP_RND_NEAREST 0
  34. # define FP_RND_ZERO 1
  35. # define FP_RND_PINF 2
  36. # define FP_RND_MINF 3
  37. #ifndef FP_ROUNDMODE
  38. # define FP_ROUNDMODE FP_RND_NEAREST
  39. #endif
  40. #endif
  41. /* By default don't care about exceptions. */
  42. #ifndef FP_EX_INVALID
  43. #define FP_EX_INVALID 0
  44. #endif
  45. #ifndef FP_EX_OVERFLOW
  46. #define FP_EX_OVERFLOW 0
  47. #endif
  48. #ifndef FP_EX_UNDERFLOW
  49. #define FP_EX_UNDERFLOW
  50. #endif
  51. #ifndef FP_EX_DIVZERO
  52. #define FP_EX_DIVZERO 0
  53. #endif
  54. #ifndef FP_EX_INEXACT
  55. #define FP_EX_INEXACT 0
  56. #endif
  57. #ifndef FP_EX_DENORM
  58. #define FP_EX_DENORM 0
  59. #endif
  60. #ifdef _FP_DECL_EX
  61. #define FP_DECL_EX
  62.   int _fex = 0;
  63.   _FP_DECL_EX
  64. #else
  65. #define FP_DECL_EX int _fex = 0
  66. #endif
  67.   
  68. #ifndef FP_INIT_ROUNDMODE
  69. #define FP_INIT_ROUNDMODE do {} while (0)
  70. #endif
  71. #ifndef FP_HANDLE_EXCEPTIONS
  72. #define FP_HANDLE_EXCEPTIONS do {} while (0)
  73. #endif
  74. /* By default we never flush denormal input operands to signed zero. */
  75. #ifndef FP_DENORM_ZERO
  76. #define FP_DENORM_ZERO 0
  77. #endif
  78. #ifndef FP_INHIBIT_RESULTS
  79. /* By default we write the results always.
  80.  * sfp-machine may override this and e.g.
  81.  * check if some exceptions are unmasked
  82.  * and inhibit it in such a case.
  83.  */
  84. #define FP_INHIBIT_RESULTS 0
  85. #endif
  86. #define FP_SET_EXCEPTION(ex)
  87.   _fex |= (ex)
  88.   
  89. #define FP_UNSET_EXCEPTION(ex)
  90.   _fex &= ~(ex)
  91. #define FP_CLEAR_EXCEPTIONS
  92.   _fex = 0
  93. #define _FP_ROUND_NEAREST(wc, X)
  94. do {
  95.     if ((_FP_FRAC_LOW_##wc(X) & 15) != _FP_WORK_ROUND)
  96.       _FP_FRAC_ADDI_##wc(X, _FP_WORK_ROUND);
  97. } while (0)
  98. #define _FP_ROUND_ZERO(wc, X) 0
  99. #define _FP_ROUND_PINF(wc, X)
  100. do {
  101.     if (!X##_s && (_FP_FRAC_LOW_##wc(X) & 7))
  102.       _FP_FRAC_ADDI_##wc(X, _FP_WORK_LSB);
  103. } while (0)
  104. #define _FP_ROUND_MINF(wc, X)
  105. do {
  106.     if (X##_s && (_FP_FRAC_LOW_##wc(X) & 7))
  107.       _FP_FRAC_ADDI_##wc(X, _FP_WORK_LSB);
  108. } while (0)
  109. #define _FP_ROUND(wc, X)
  110. do {
  111. if (_FP_FRAC_LOW_##wc(X) & 7)
  112.   FP_SET_EXCEPTION(FP_EX_INEXACT);
  113. switch (FP_ROUNDMODE)
  114. {
  115.   case FP_RND_NEAREST:
  116.     _FP_ROUND_NEAREST(wc,X);
  117.     break;
  118.   case FP_RND_ZERO:
  119.     _FP_ROUND_ZERO(wc,X);
  120.     break;
  121.   case FP_RND_PINF:
  122.     _FP_ROUND_PINF(wc,X);
  123.     break;
  124.   case FP_RND_MINF:
  125.     _FP_ROUND_MINF(wc,X);
  126.     break;
  127. }
  128. } while (0)
  129. #define FP_CLS_NORMAL 0
  130. #define FP_CLS_ZERO 1
  131. #define FP_CLS_INF 2
  132. #define FP_CLS_NAN 3
  133. #define _FP_CLS_COMBINE(x,y) (((x) << 2) | (y))
  134. #include <math-emu/op-1.h>
  135. #include <math-emu/op-2.h>
  136. #include <math-emu/op-4.h>
  137. #include <math-emu/op-8.h>
  138. #include <math-emu/op-common.h>
  139. /* Sigh.  Silly things longlong.h needs.  */
  140. #define UWtype _FP_W_TYPE
  141. #define W_TYPE_SIZE _FP_W_TYPE_SIZE
  142. typedef int SItype __attribute__((mode(SI)));
  143. typedef int DItype __attribute__((mode(DI)));
  144. typedef unsigned int USItype __attribute__((mode(SI)));
  145. typedef unsigned int UDItype __attribute__((mode(DI)));
  146. #if _FP_W_TYPE_SIZE == 32
  147. typedef unsigned int UHWtype __attribute__((mode(HI)));
  148. #elif _FP_W_TYPE_SIZE == 64
  149. typedef USItype UHWtype;
  150. #endif
  151. #ifndef umul_ppmm
  152. #include <stdlib/longlong.h>
  153. #endif
  154. #endif /* __MATH_EMU_SOFT_FP_H__ */