/* 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_INVALID_SNAN
  #define FP_EX_INVALID_SNAN	0
  #endif
  /* inf - inf */
  #ifndef FP_EX_INVALID_ISI
  #define FP_EX_INVALID_ISI	0
  #endif
  /* inf / inf */
  #ifndef FP_EX_INVALID_IDI
  #define FP_EX_INVALID_IDI	0
  #endif
  /* 0 / 0 */
  #ifndef FP_EX_INVALID_ZDZ
  #define FP_EX_INVALID_ZDZ	0
  #endif
  /* inf * 0 */
  #ifndef FP_EX_INVALID_IMZ
  #define FP_EX_INVALID_IMZ	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_EXCEPTIONS 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__ */