Subversion Repositories or1k

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

#define FP_SET_EXCEPTION(ex)                            \

  _fex |= (ex)

#define FP_UNSET_EXCEPTION(ex)                          \

  _fex &= ~(ex)

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