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