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#define FP_SET_EXCEPTION(ex) \ 101 _fex |= (ex) 102 103#define FP_UNSET_EXCEPTION(ex) \ 104 _fex &= ~(ex) 105 106#define FP_CLEAR_EXCEPTIONS \ 107 _fex = 0 108 109#define _FP_ROUND_NEAREST(wc, X) \ 110do { \ 111 if ((_FP_FRAC_LOW_##wc(X) & 15) != _FP_WORK_ROUND) \ 112 _FP_FRAC_ADDI_##wc(X, _FP_WORK_ROUND); \ 113} while (0) 114 115#define _FP_ROUND_ZERO(wc, X) 0 116 117#define _FP_ROUND_PINF(wc, X) \ 118do { \ 119 if (!X##_s && (_FP_FRAC_LOW_##wc(X) & 7)) \ 120 _FP_FRAC_ADDI_##wc(X, _FP_WORK_LSB); \ 121} while (0) 122 123#define _FP_ROUND_MINF(wc, X) \ 124do { \ 125 if (X##_s && (_FP_FRAC_LOW_##wc(X) & 7)) \ 126 _FP_FRAC_ADDI_##wc(X, _FP_WORK_LSB); \ 127} while (0) 128 129#define _FP_ROUND(wc, X) \ 130do { \ 131 if (_FP_FRAC_LOW_##wc(X) & 7) \ 132 FP_SET_EXCEPTION(FP_EX_INEXACT); \ 133 switch (FP_ROUNDMODE) \ 134 { \ 135 case FP_RND_NEAREST: \ 136 _FP_ROUND_NEAREST(wc,X); \ 137 break; \ 138 case FP_RND_ZERO: \ 139 _FP_ROUND_ZERO(wc,X); \ 140 break; \ 141 case FP_RND_PINF: \ 142 _FP_ROUND_PINF(wc,X); \ 143 break; \ 144 case FP_RND_MINF: \ 145 _FP_ROUND_MINF(wc,X); \ 146 break; \ 147 } \ 148} while (0) 149 150#define FP_CLS_NORMAL 0 151#define FP_CLS_ZERO 1 152#define FP_CLS_INF 2 153#define FP_CLS_NAN 3 154 155#define _FP_CLS_COMBINE(x,y) (((x) << 2) | (y)) 156 157#include <math-emu/op-1.h> 158#include <math-emu/op-2.h> 159#include <math-emu/op-4.h> 160#include <math-emu/op-8.h> 161#include <math-emu/op-common.h> 162 163/* Sigh. Silly things longlong.h needs. */ 164#define UWtype _FP_W_TYPE 165#define W_TYPE_SIZE _FP_W_TYPE_SIZE 166 167typedef int SItype __attribute__((mode(SI))); 168typedef int DItype __attribute__((mode(DI))); 169typedef unsigned int USItype __attribute__((mode(SI))); 170typedef unsigned int UDItype __attribute__((mode(DI))); 171#if _FP_W_TYPE_SIZE == 32 172typedef unsigned int UHWtype __attribute__((mode(HI))); 173#elif _FP_W_TYPE_SIZE == 64 174typedef USItype UHWtype; 175#endif 176 177#ifndef umul_ppmm 178#include <stdlib/longlong.h> 179#endif 180 181#endif /* __MATH_EMU_SOFT_FP_H__ */ 182