SAFE public projects git trees. - safe/jmp/linux-2.6/blob - include/math-emu/soft-fp.h

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