| /* Software floating-point emulation. |
| Copyright (C) 1997,1998,1999 Free Software Foundation, Inc. |
| This file is part of the GNU C Library. |
| Contributed by Richard Henderson (rth@cygnus.com), |
| Jakub Jelinek (jj@ultra.linux.cz), |
| David S. Miller (davem@redhat.com) and |
| Peter Maydell (pmaydell@chiark.greenend.org.uk). |
| |
| The GNU C Library is free software; you can redistribute it and/or |
| modify it under the terms of the GNU Library General Public License as |
| published by the Free Software Foundation; either version 2 of the |
| License, or (at your option) any later version. |
| |
| The GNU C Library is distributed in the hope that it will be useful, |
| but WITHOUT ANY WARRANTY; without even the implied warranty of |
| MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
| Library General Public License for more details. |
| |
| You should have received a copy of the GNU Library General Public |
| License along with the GNU C Library; see the file COPYING.LIB. If |
| not, write to the Free Software Foundation, Inc., |
| 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */ |
| |
| #ifndef __MATH_EMU_SOFT_FP_H__ |
| #define __MATH_EMU_SOFT_FP_H__ |
| |
| #include <asm/sfp-machine.h> |
| |
| /* Allow sfp-machine to have its own byte order definitions. */ |
| #ifndef __BYTE_ORDER |
| #include <endian.h> |
| #endif |
| |
| #define _FP_WORKBITS 3 |
| #define _FP_WORK_LSB ((_FP_W_TYPE)1 << 3) |
| #define _FP_WORK_ROUND ((_FP_W_TYPE)1 << 2) |
| #define _FP_WORK_GUARD ((_FP_W_TYPE)1 << 1) |
| #define _FP_WORK_STICKY ((_FP_W_TYPE)1 << 0) |
| |
| #ifndef FP_RND_NEAREST |
| # define FP_RND_NEAREST 0 |
| # define FP_RND_ZERO 1 |
| # define FP_RND_PINF 2 |
| # define FP_RND_MINF 3 |
| #ifndef FP_ROUNDMODE |
| # define FP_ROUNDMODE FP_RND_NEAREST |
| #endif |
| #endif |
| |
| /* By default don't care about exceptions. */ |
| #ifndef FP_EX_INVALID |
| #define FP_EX_INVALID 0 |
| #endif |
| #ifndef FP_EX_OVERFLOW |
| #define FP_EX_OVERFLOW 0 |
| #endif |
| #ifndef FP_EX_UNDERFLOW |
| #define FP_EX_UNDERFLOW |
| #endif |
| #ifndef FP_EX_DIVZERO |
| #define FP_EX_DIVZERO 0 |
| #endif |
| #ifndef FP_EX_INEXACT |
| #define FP_EX_INEXACT 0 |
| #endif |
| #ifndef FP_EX_DENORM |
| #define FP_EX_DENORM 0 |
| #endif |
| |
| #ifdef _FP_DECL_EX |
| #define FP_DECL_EX \ |
| int _fex = 0; \ |
| _FP_DECL_EX |
| #else |
| #define FP_DECL_EX int _fex = 0 |
| #endif |
| |
| #ifndef FP_INIT_ROUNDMODE |
| #define FP_INIT_ROUNDMODE do {} while (0) |
| #endif |
| |
| #ifndef FP_HANDLE_EXCEPTIONS |
| #define FP_HANDLE_EXCEPTIONS do {} while (0) |
| #endif |
| |
| /* By default we never flush denormal input operands to signed zero. */ |
| #ifndef FP_DENORM_ZERO |
| #define FP_DENORM_ZERO 0 |
| #endif |
| |
| #ifndef FP_INHIBIT_RESULTS |
| /* By default we write the results always. |
| * sfp-machine may override this and e.g. |
| * check if some exceptions are unmasked |
| * and inhibit it in such a case. |
| */ |
| #define FP_INHIBIT_RESULTS 0 |
| #endif |
| |
| #ifndef FP_TRAPPING_EXCEPTIONS |
| #define FP_TRAPPING_EXCPETIONS 0 |
| #endif |
| |
| #define FP_SET_EXCEPTION(ex) \ |
| _fex |= (ex) |
| |
| #define FP_UNSET_EXCEPTION(ex) \ |
| _fex &= ~(ex) |
| |
| #define FP_CUR_EXCEPTIONS \ |
| (_fex) |
| |
| #define FP_CLEAR_EXCEPTIONS \ |
| _fex = 0 |
| |
| #define _FP_ROUND_NEAREST(wc, X) \ |
| do { \ |
| if ((_FP_FRAC_LOW_##wc(X) & 15) != _FP_WORK_ROUND) \ |
| _FP_FRAC_ADDI_##wc(X, _FP_WORK_ROUND); \ |
| } while (0) |
| |
| #define _FP_ROUND_ZERO(wc, X) 0 |
| |
| #define _FP_ROUND_PINF(wc, X) \ |
| do { \ |
| if (!X##_s && (_FP_FRAC_LOW_##wc(X) & 7)) \ |
| _FP_FRAC_ADDI_##wc(X, _FP_WORK_LSB); \ |
| } while (0) |
| |
| #define _FP_ROUND_MINF(wc, X) \ |
| do { \ |
| if (X##_s && (_FP_FRAC_LOW_##wc(X) & 7)) \ |
| _FP_FRAC_ADDI_##wc(X, _FP_WORK_LSB); \ |
| } while (0) |
| |
| #define _FP_ROUND(wc, X) \ |
| do { \ |
| if (_FP_FRAC_LOW_##wc(X) & 7) \ |
| FP_SET_EXCEPTION(FP_EX_INEXACT); \ |
| switch (FP_ROUNDMODE) \ |
| { \ |
| case FP_RND_NEAREST: \ |
| _FP_ROUND_NEAREST(wc,X); \ |
| break; \ |
| case FP_RND_ZERO: \ |
| _FP_ROUND_ZERO(wc,X); \ |
| break; \ |
| case FP_RND_PINF: \ |
| _FP_ROUND_PINF(wc,X); \ |
| break; \ |
| case FP_RND_MINF: \ |
| _FP_ROUND_MINF(wc,X); \ |
| break; \ |
| } \ |
| } while (0) |
| |
| #define FP_CLS_NORMAL 0 |
| #define FP_CLS_ZERO 1 |
| #define FP_CLS_INF 2 |
| #define FP_CLS_NAN 3 |
| |
| #define _FP_CLS_COMBINE(x,y) (((x) << 2) | (y)) |
| |
| #include <math-emu/op-1.h> |
| #include <math-emu/op-2.h> |
| #include <math-emu/op-4.h> |
| #include <math-emu/op-8.h> |
| #include <math-emu/op-common.h> |
| |
| /* Sigh. Silly things longlong.h needs. */ |
| #define UWtype _FP_W_TYPE |
| #define W_TYPE_SIZE _FP_W_TYPE_SIZE |
| |
| typedef int SItype __attribute__((mode(SI))); |
| typedef int DItype __attribute__((mode(DI))); |
| typedef unsigned int USItype __attribute__((mode(SI))); |
| typedef unsigned int UDItype __attribute__((mode(DI))); |
| #if _FP_W_TYPE_SIZE == 32 |
| typedef unsigned int UHWtype __attribute__((mode(HI))); |
| #elif _FP_W_TYPE_SIZE == 64 |
| typedef USItype UHWtype; |
| #endif |
| |
| #ifndef umul_ppmm |
| #include <stdlib/longlong.h> |
| #endif |
| |
| #endif /* __MATH_EMU_SOFT_FP_H__ */ |