linux/include/math-emu/soft-fp.h
<<
>>
Prefs
   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_INVALID_SNAN
  55#define FP_EX_INVALID_SNAN      0
  56#endif
  57/* inf - inf */
  58#ifndef FP_EX_INVALID_ISI
  59#define FP_EX_INVALID_ISI       0
  60#endif
  61/* inf / inf */
  62#ifndef FP_EX_INVALID_IDI
  63#define FP_EX_INVALID_IDI       0
  64#endif
  65/* 0 / 0 */
  66#ifndef FP_EX_INVALID_ZDZ
  67#define FP_EX_INVALID_ZDZ       0
  68#endif
  69/* inf * 0 */
  70#ifndef FP_EX_INVALID_IMZ
  71#define FP_EX_INVALID_IMZ       0
  72#endif
  73#ifndef FP_EX_OVERFLOW
  74#define FP_EX_OVERFLOW          0
  75#endif
  76#ifndef FP_EX_UNDERFLOW
  77#define FP_EX_UNDERFLOW         
  78#endif
  79#ifndef FP_EX_DIVZERO
  80#define FP_EX_DIVZERO           0
  81#endif
  82#ifndef FP_EX_INEXACT
  83#define FP_EX_INEXACT           0
  84#endif
  85#ifndef FP_EX_DENORM
  86#define FP_EX_DENORM            0
  87#endif
  88
  89#ifdef _FP_DECL_EX
  90#define FP_DECL_EX                                      \
  91  int _fex = 0;                                         \
  92  _FP_DECL_EX
  93#else
  94#define FP_DECL_EX int _fex = 0
  95#endif
  96  
  97#ifndef FP_INIT_ROUNDMODE
  98#define FP_INIT_ROUNDMODE do {} while (0)
  99#endif
 100
 101#ifndef FP_HANDLE_EXCEPTIONS
 102#define FP_HANDLE_EXCEPTIONS do {} while (0)
 103#endif
 104
 105/* By default we never flush denormal input operands to signed zero. */
 106#ifndef FP_DENORM_ZERO
 107#define FP_DENORM_ZERO 0
 108#endif
 109
 110#ifndef FP_INHIBIT_RESULTS
 111/* By default we write the results always.
 112 * sfp-machine may override this and e.g.
 113 * check if some exceptions are unmasked
 114 * and inhibit it in such a case.
 115 */
 116#define FP_INHIBIT_RESULTS 0
 117#endif
 118
 119#ifndef FP_TRAPPING_EXCEPTIONS
 120#define FP_TRAPPING_EXCEPTIONS 0
 121#endif
 122
 123#define FP_SET_EXCEPTION(ex)                            \
 124  _fex |= (ex)
 125  
 126#define FP_UNSET_EXCEPTION(ex)                          \
 127  _fex &= ~(ex)
 128
 129#define FP_CUR_EXCEPTIONS                               \
 130  (_fex)
 131
 132#define FP_CLEAR_EXCEPTIONS                             \
 133  _fex = 0
 134
 135#define _FP_ROUND_NEAREST(wc, X)                        \
 136do {                                                    \
 137    if ((_FP_FRAC_LOW_##wc(X) & 15) != _FP_WORK_ROUND)  \
 138      _FP_FRAC_ADDI_##wc(X, _FP_WORK_ROUND);            \
 139} while (0)
 140
 141#define _FP_ROUND_ZERO(wc, X)           (void)0
 142
 143#define _FP_ROUND_PINF(wc, X)                           \
 144do {                                                    \
 145    if (!X##_s && (_FP_FRAC_LOW_##wc(X) & 7))           \
 146      _FP_FRAC_ADDI_##wc(X, _FP_WORK_LSB);              \
 147} while (0)
 148
 149#define _FP_ROUND_MINF(wc, X)                           \
 150do {                                                    \
 151    if (X##_s && (_FP_FRAC_LOW_##wc(X) & 7))            \
 152      _FP_FRAC_ADDI_##wc(X, _FP_WORK_LSB);              \
 153} while (0)
 154
 155#define _FP_ROUND(wc, X)                        \
 156do {                                            \
 157        if (_FP_FRAC_LOW_##wc(X) & 7)           \
 158          FP_SET_EXCEPTION(FP_EX_INEXACT);      \
 159        switch (FP_ROUNDMODE)                   \
 160        {                                       \
 161          case FP_RND_NEAREST:                  \
 162            _FP_ROUND_NEAREST(wc,X);            \
 163            break;                              \
 164          case FP_RND_ZERO:                     \
 165            _FP_ROUND_ZERO(wc,X);               \
 166            break;                              \
 167          case FP_RND_PINF:                     \
 168            _FP_ROUND_PINF(wc,X);               \
 169            break;                              \
 170          case FP_RND_MINF:                     \
 171            _FP_ROUND_MINF(wc,X);               \
 172            break;                              \
 173        }                                       \
 174} while (0)
 175
 176#define FP_CLS_NORMAL           0
 177#define FP_CLS_ZERO             1
 178#define FP_CLS_INF              2
 179#define FP_CLS_NAN              3
 180
 181#define _FP_CLS_COMBINE(x,y)    (((x) << 2) | (y))
 182
 183#include <math-emu/op-1.h>
 184#include <math-emu/op-2.h>
 185#include <math-emu/op-4.h>
 186#include <math-emu/op-8.h>
 187#include <math-emu/op-common.h>
 188
 189/* Sigh.  Silly things longlong.h needs.  */
 190#define UWtype          _FP_W_TYPE
 191#define W_TYPE_SIZE     _FP_W_TYPE_SIZE
 192
 193typedef int SItype __attribute__((mode(SI)));
 194typedef int DItype __attribute__((mode(DI)));
 195typedef unsigned int USItype __attribute__((mode(SI)));
 196typedef unsigned int UDItype __attribute__((mode(DI)));
 197#if _FP_W_TYPE_SIZE == 32
 198typedef unsigned int UHWtype __attribute__((mode(HI)));
 199#elif _FP_W_TYPE_SIZE == 64
 200typedef USItype UHWtype;
 201#endif
 202
 203#ifndef umul_ppmm
 204#include <stdlib/longlong.h>
 205#endif
 206
 207#endif /* __MATH_EMU_SOFT_FP_H__ */
 208