307 lines
12 KiB
C
307 lines
12 KiB
C
#ifndef _FLOAT_
|
|
#define _FLOAT_
|
|
|
|
#ifdef __cplusplus
|
|
extern "C" {
|
|
#endif
|
|
|
|
#include <ansidecl.h>
|
|
|
|
|
|
#ifndef _CRTIMP
|
|
#define _CRTIMP
|
|
#endif
|
|
|
|
extern u8 _Nan[8];
|
|
extern u8 _Inf[8];
|
|
|
|
//#define NAN ((unsigned __int64)0xfff8000000000000)
|
|
|
|
#define NAN *((double*)(&_Nan[0])) // instead of using __int64, store a NAN somwhere in memory (misc_vars.c)
|
|
#define INF *((double*)(&_Inf[0])) // Constant INF
|
|
|
|
#define DBL_DIG 15 // # of decimal digits of precision
|
|
#define DBL_EPSILON 2.2204460492503131e-016 // smallest such that 1.0+DBL_EPSILON != 1.0
|
|
#define DBL_MANT_DIG 53 // # of bits in mantissa
|
|
#define DBL_MAX 1.7976931348623158e+308 // max value
|
|
#define DBL_MAX_10_EXP 308 // max decimal exponent
|
|
#define DBL_MAX_EXP 1024 // max binary exponent
|
|
#define DBL_MIN 2.2250738585072014e-308 // min positive value
|
|
#define DBL_MIN_10_EXP (-307) // min decimal exponent
|
|
#define DBL_MIN_EXP (-1021) // min binary exponent
|
|
#define _DBL_RADIX 2 // exponent radix
|
|
#define _DBL_ROUNDS 1 // addition rounding: near
|
|
|
|
#define FLT_DIG 6 // # of decimal digits of precision
|
|
#define FLT_EPSILON 1.192092896e-07F // smallest such that 1.0+FLT_EPSILON != 1.0
|
|
#define FLT_GUARD 0
|
|
#define FLT_MANT_DIG 24 // # of bits in mantissa
|
|
#define FLT_MAX 3.402823466e+38F // max value
|
|
#define FLT_MAX_10_EXP 38 // max decimal exponent
|
|
#define FLT_MAX_EXP 128 // max binary exponent
|
|
#define FLT_MIN 1.175494351e-38F // min positive value
|
|
#define FLT_MIN_10_EXP (-37) // min decimal exponent
|
|
#define FLT_MIN_EXP (-125) // min binary exponent
|
|
#define FLT_NORMALIZE 0
|
|
#define FLT_RADIX 2 // exponent radix
|
|
#define FLT_ROUNDS 1 // addition rounding: near
|
|
|
|
#ifndef _M_M68K
|
|
#define LDBL_DIG DBL_DIG // # of decimal digits of precision
|
|
#define LDBL_EPSILON DBL_EPSILON // smallest such that 1.0+LDBL_EPSILON != 1.0
|
|
#define LDBL_MANT_DIG DBL_MANT_DIG // # of bits in mantissa
|
|
#define LDBL_MAX DBL_MAX // max value
|
|
#define LDBL_MAX_10_EXP DBL_MAX_10_EXP // max decimal exponent
|
|
#define LDBL_MAX_EXP DBL_MAX_EXP // max binary exponent
|
|
#define LDBL_MIN DBL_MIN // min positive value
|
|
#define LDBL_MIN_10_EXP DBL_MIN_10_EXP // min decimal exponent
|
|
#define LDBL_MIN_EXP DBL_MIN_EXP // min binary exponent
|
|
#define _LDBL_RADIX DBL_RADIX // exponent radix
|
|
#define _LDBL_ROUNDS DBL_ROUNDS // addition rounding: near
|
|
#else
|
|
#define LDBL_DIG 18 // # of decimal digits of precision
|
|
#define LDBL_EPSILON 1.08420217248550443412e-019L // smallest such that 1.0+LDBL_EPSILON != 1.0
|
|
#define LDBL_MANT_DIG 64 // # of bits in mantissa
|
|
#define LDBL_MAX 1.189731495357231765e+4932L // max value
|
|
#define LDBL_MAX_10_EXP 4932 // max decimal exponent
|
|
#define LDBL_MAX_EXP 16384 // max binary exponent
|
|
#define LDBL_MIN 3.3621031431120935063e-4932L // min positive value
|
|
#define LDBL_MIN_10_EXP (-4931) // min decimal exponent
|
|
#define LDBL_MIN_EXP (-16381) // min binary exponent
|
|
#define _LDBL_RADIX 2 // exponent radix
|
|
#define _LDBL_ROUNDS 1 // addition rounding: near
|
|
#endif
|
|
|
|
// Function prototypes
|
|
|
|
_CRTIMP unsigned int __cdecl _clearfp(void);
|
|
_CRTIMP unsigned int __cdecl _controlfp(unsigned int,unsigned int);
|
|
_CRTIMP unsigned int __cdecl _statusfp(void);
|
|
_CRTIMP void __cdecl _fpreset(void);
|
|
|
|
#ifndef _MAC
|
|
#define _clear87 _clearfp
|
|
#define _status87 _statusfp
|
|
#endif // _MAC
|
|
|
|
//
|
|
// Abstract User Control Word Mask and bit definitions
|
|
//
|
|
|
|
#ifndef _MAC
|
|
#define _MCW_EM 0x0008001f // interrupt Exception Masks
|
|
#else
|
|
#define _MCW_EM 0x0000001f // interrupt Exception Masks
|
|
#endif
|
|
#define _EM_INEXACT 0x00000001 // inexact (precision)
|
|
#define _EM_UNDERFLOW 0x00000002 // underflow
|
|
#define _EM_OVERFLOW 0x00000004 // overflow
|
|
#define _EM_ZERODIVIDE 0x00000008 // zero divide
|
|
#define _EM_INVALID 0x00000010 // invalid
|
|
|
|
#define _MCW_RC 0x00000300 // Rounding Control
|
|
#define _RC_NEAR 0x00000000 // near
|
|
#define _RC_DOWN 0x00000100 // down
|
|
#define _RC_UP 0x00000200 // up
|
|
#define _RC_CHOP 0x00000300 // chop
|
|
|
|
//
|
|
// Abstract User Status Word bit definitions
|
|
|
|
|
|
#define _SW_INEXACT 0x00000001 // inexact (precision)
|
|
#define _SW_UNDERFLOW 0x00000002 // underflow
|
|
#define _SW_OVERFLOW 0x00000004 // overflow
|
|
#define _SW_ZERODIVIDE 0x00000008 // zero divide
|
|
#define _SW_INVALID 0x00000010 // invalid
|
|
|
|
|
|
//
|
|
// i386 specific definitions
|
|
|
|
#define _MCW_PC 0x00030000 // Precision Control
|
|
#ifdef _M_MPPC
|
|
//
|
|
* PowerMac specific definitions(no precision control)
|
|
|
|
#define _PC_64 0x00000000 // 64 bits
|
|
#define _PC_53 0x00000000 // 53 bits
|
|
#define _PC_24 0x00000000 // 24 bits
|
|
#else
|
|
#define _PC_64 0x00000000 // 64 bits
|
|
#define _PC_53 0x00010000 // 53 bits
|
|
#define _PC_24 0x00020000 // 24 bits
|
|
#endif
|
|
|
|
#define _MCW_IC 0x00040000 // Infinity Control
|
|
#define _IC_AFFINE 0x00040000 // affine
|
|
#define _IC_PROJECTIVE 0x00000000 // projective
|
|
|
|
#define _EM_DENORMAL 0x00080000 // denormal exception mask (_control87 only)
|
|
|
|
#define _SW_DENORMAL 0x00080000 // denormal status bit
|
|
|
|
|
|
_CRTIMP unsigned int __cdecl _control87(unsigned int,unsigned int);
|
|
|
|
|
|
//
|
|
// RISC specific definitions
|
|
|
|
|
|
#define _MCW_DN 0x03000000 // Denormal Control
|
|
#define _DN_SAVE 0x00000000 // save denormal results and operands
|
|
#define _DN_FLUSH 0x01000000 // flush denormal results and operands to zero
|
|
#define _DN_FLUSH_OPERANDS_SAVE_RESULTS 0x02000000 // flush operands to zero and save results
|
|
#define _DN_SAVE_OPERANDS_FLUSH_RESULTS 0x03000000 // save operands and flush results to zero
|
|
|
|
|
|
// initial Control Word value
|
|
|
|
#if defined(_M_IX86)
|
|
|
|
#define _CW_DEFAULT ( _RC_NEAR + _PC_53 + _EM_INVALID + _EM_ZERODIVIDE + _EM_OVERFLOW + _EM_UNDERFLOW + _EM_INEXACT + _EM_DENORMAL)
|
|
|
|
#elif defined(_MAC)
|
|
|
|
#define _CW_DEFAULT ( _RC_NEAR + _PC_64 + _EM_INVALID + _EM_ZERODIVIDE + _EM_OVERFLOW + _EM_UNDERFLOW + _EM_INEXACT )
|
|
|
|
#elif defined(_M_MRX000) || defined (_M_ALPHA) || defined(_M_PPC)
|
|
|
|
#define _CW_DEFAULT ( _RC_NEAR + _DN_FLUSH + _EM_INVALID + _EM_ZERODIVIDE + _EM_OVERFLOW + _EM_UNDERFLOW + _EM_INEXACT )
|
|
|
|
#endif
|
|
|
|
// Global variable holding floating point error code
|
|
|
|
#if defined(_MT) || defined(_DLL)
|
|
_CRTIMP extern int * __cdecl __fpecode(void);
|
|
#define _fpecode (*__fpecode())
|
|
#else // ndef _MT && ndef _DLL
|
|
extern int _fpecode;
|
|
#endif // _MT || _DLL
|
|
|
|
// invalid subconditions (_SW_INVALID also set)
|
|
|
|
#define _SW_UNEMULATED 0x0040 // unemulated instruction
|
|
#define _SW_SQRTNEG 0x0080 // square root of a neg number
|
|
#define _SW_STACKOVERFLOW 0x0200 // FP stack overflow
|
|
#define _SW_STACKUNDERFLOW 0x0400 // FP stack underflow
|
|
|
|
// Floating point error signals and return codes
|
|
|
|
#define _FPE_INVALID 0x81
|
|
#define _FPE_DENORMAL 0x82
|
|
#define _FPE_ZERODIVIDE 0x83
|
|
#define _FPE_OVERFLOW 0x84
|
|
#define _FPE_UNDERFLOW 0x85
|
|
#define _FPE_INEXACT 0x86
|
|
|
|
#define _FPE_UNEMULATED 0x87
|
|
#define _FPE_SQRTNEG 0x88
|
|
#define _FPE_STACKOVERFLOW 0x8a
|
|
#define _FPE_STACKUNDERFLOW 0x8b
|
|
|
|
#define _FPE_EXPLICITGEN 0x8c // raise( SIGFPE );
|
|
|
|
|
|
// IEEE recommended functions
|
|
|
|
_CRTIMP double __cdecl _copysign (double, double);
|
|
_CRTIMP double __cdecl _chgsign (double);
|
|
_CRTIMP double __cdecl _scalb(double, long);
|
|
_CRTIMP double __cdecl _logb(double);
|
|
_CRTIMP double __cdecl _nextafter(double, double);
|
|
_CRTIMP int __cdecl _finite(double);
|
|
_CRTIMP int __cdecl _isnan(double);
|
|
_CRTIMP int __cdecl _fpclass(double);
|
|
|
|
#define _FPCLASS_SNAN 0x0001 // signaling NaN
|
|
#define _FPCLASS_QNAN 0x0002 // quiet NaN
|
|
#define _FPCLASS_NINF 0x0004 // negative infinity
|
|
#define _FPCLASS_NN 0x0008 // negative normal
|
|
#define _FPCLASS_ND 0x0010 // negative denormal
|
|
#define _FPCLASS_NZ 0x0020 // -0
|
|
#define _FPCLASS_PZ 0x0040 // +0
|
|
#define _FPCLASS_PD 0x0080 // positive denormal
|
|
#define _FPCLASS_PN 0x0100 // positive normal
|
|
#define _FPCLASS_PINF 0x0200 // positive infinity
|
|
|
|
|
|
#if !__STDC__
|
|
|
|
// Non-ANSI names for compatibility
|
|
|
|
#ifndef _MAC
|
|
#define clear87 _clear87
|
|
#define status87 _status87
|
|
#define control87 _control87
|
|
#endif // _MAC
|
|
|
|
_CRTIMP void __cdecl fpreset(void);
|
|
|
|
#define DBL_RADIX _DBL_RADIX
|
|
#define DBL_ROUNDS _DBL_ROUNDS
|
|
|
|
#define LDBL_RADIX _LDBL_RADIX
|
|
#define LDBL_ROUNDS _LDBL_ROUNDS
|
|
|
|
#define MCW_EM _MCW_EM
|
|
#define EM_INVALID _EM_INVALID
|
|
#define EM_DENORMAL _EM_DENORMAL
|
|
#define EM_ZERODIVIDE _EM_ZERODIVIDE
|
|
#define EM_OVERFLOW _EM_OVERFLOW
|
|
#define EM_UNDERFLOW _EM_UNDERFLOW
|
|
#define EM_INEXACT _EM_INEXACT
|
|
|
|
#define MCW_IC _MCW_IC
|
|
#define IC_AFFINE _IC_AFFINE
|
|
#define IC_PROJECTIVE _IC_PROJECTIVE
|
|
|
|
#define MCW_RC _MCW_RC
|
|
#define RC_CHOP _RC_CHOP
|
|
#define RC_UP _RC_UP
|
|
#define RC_DOWN _RC_DOWN
|
|
#define RC_NEAR _RC_NEAR
|
|
|
|
#define MCW_PC _MCW_PC
|
|
#define PC_24 _PC_24
|
|
#define PC_53 _PC_53
|
|
#define PC_64 _PC_64
|
|
|
|
#define CW_DEFAULT _CW_DEFAULT
|
|
|
|
#define SW_INVALID _SW_INVALID
|
|
#define SW_DENORMAL _SW_DENORMAL
|
|
#define SW_ZERODIVIDE _SW_ZERODIVIDE
|
|
#define SW_OVERFLOW _SW_OVERFLOW
|
|
#define SW_UNDERFLOW _SW_UNDERFLOW
|
|
#define SW_INEXACT _SW_INEXACT
|
|
|
|
#define SW_UNEMULATED _SW_UNEMULATED
|
|
#define SW_SQRTNEG _SW_SQRTNEG
|
|
#define SW_STACKOVERFLOW _SW_STACKOVERFLOW
|
|
#define SW_STACKUNDERFLOW _SW_STACKUNDERFLOW
|
|
|
|
#define FPE_INVALID _FPE_INVALID
|
|
#define FPE_DENORMAL _FPE_DENORMAL
|
|
#define FPE_ZERODIVIDE _FPE_ZERODIVIDE
|
|
#define FPE_OVERFLOW _FPE_OVERFLOW
|
|
#define FPE_UNDERFLOW _FPE_UNDERFLOW
|
|
#define FPE_INEXACT _FPE_INEXACT
|
|
|
|
#define FPE_UNEMULATED _FPE_UNEMULATED
|
|
#define FPE_SQRTNEG _FPE_SQRTNEG
|
|
#define FPE_STACKOVERFLOW _FPE_STACKOVERFLOW
|
|
#define FPE_STACKUNDERFLOW _FPE_STACKUNDERFLOW
|
|
|
|
#define FPE_EXPLICITGEN _FPE_EXPLICITGEN
|
|
|
|
#endif
|
|
|
|
#ifdef __cplusplus
|
|
}
|
|
#endif
|
|
|
|
#endif
|