1*840b91ccSNathan Whitehorn /* $NetBSD: softfloat.h,v 1.6 2002/05/12 13:12:46 bjh21 Exp $ */ 2*840b91ccSNathan Whitehorn /* $FreeBSD$ */ 3*840b91ccSNathan Whitehorn 4*840b91ccSNathan Whitehorn /* This is a derivative work. */ 5*840b91ccSNathan Whitehorn 6*840b91ccSNathan Whitehorn /* 7*840b91ccSNathan Whitehorn =============================================================================== 8*840b91ccSNathan Whitehorn 9*840b91ccSNathan Whitehorn This C header file is part of the SoftFloat IEC/IEEE Floating-point 10*840b91ccSNathan Whitehorn Arithmetic Package, Release 2a. 11*840b91ccSNathan Whitehorn 12*840b91ccSNathan Whitehorn Written by John R. Hauser. This work was made possible in part by the 13*840b91ccSNathan Whitehorn International Computer Science Institute, located at Suite 600, 1947 Center 14*840b91ccSNathan Whitehorn Street, Berkeley, California 94704. Funding was partially provided by the 15*840b91ccSNathan Whitehorn National Science Foundation under grant MIP-9311980. The original version 16*840b91ccSNathan Whitehorn of this code was written as part of a project to build a fixed-point vector 17*840b91ccSNathan Whitehorn processor in collaboration with the University of California at Berkeley, 18*840b91ccSNathan Whitehorn overseen by Profs. Nelson Morgan and John Wawrzynek. More information 19*840b91ccSNathan Whitehorn is available through the Web page `http://HTTP.CS.Berkeley.EDU/~jhauser/ 20*840b91ccSNathan Whitehorn arithmetic/SoftFloat.html'. 21*840b91ccSNathan Whitehorn 22*840b91ccSNathan Whitehorn THIS SOFTWARE IS DISTRIBUTED AS IS, FOR FREE. Although reasonable effort 23*840b91ccSNathan Whitehorn has been made to avoid it, THIS SOFTWARE MAY CONTAIN FAULTS THAT WILL AT 24*840b91ccSNathan Whitehorn TIMES RESULT IN INCORRECT BEHAVIOR. USE OF THIS SOFTWARE IS RESTRICTED TO 25*840b91ccSNathan Whitehorn PERSONS AND ORGANIZATIONS WHO CAN AND WILL TAKE FULL RESPONSIBILITY FOR ANY 26*840b91ccSNathan Whitehorn AND ALL LOSSES, COSTS, OR OTHER PROBLEMS ARISING FROM ITS USE. 27*840b91ccSNathan Whitehorn 28*840b91ccSNathan Whitehorn Derivative works are acceptable, even for commercial purposes, so long as 29*840b91ccSNathan Whitehorn (1) they include prominent notice that the work is derivative, and (2) they 30*840b91ccSNathan Whitehorn include prominent notice akin to these four paragraphs for those parts of 31*840b91ccSNathan Whitehorn this code that are retained. 32*840b91ccSNathan Whitehorn 33*840b91ccSNathan Whitehorn =============================================================================== 34*840b91ccSNathan Whitehorn */ 35*840b91ccSNathan Whitehorn 36*840b91ccSNathan Whitehorn /* 37*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 38*840b91ccSNathan Whitehorn The macro `FLOATX80' must be defined to enable the extended double-precision 39*840b91ccSNathan Whitehorn floating-point format `floatx80'. If this macro is not defined, the 40*840b91ccSNathan Whitehorn `floatx80' type will not be defined, and none of the functions that either 41*840b91ccSNathan Whitehorn input or output the `floatx80' type will be defined. The same applies to 42*840b91ccSNathan Whitehorn the `FLOAT128' macro and the quadruple-precision format `float128'. 43*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 44*840b91ccSNathan Whitehorn */ 45*840b91ccSNathan Whitehorn /* #define FLOATX80 */ 46*840b91ccSNathan Whitehorn /* #define FLOAT128 */ 47*840b91ccSNathan Whitehorn 48*840b91ccSNathan Whitehorn #include <machine/ieeefp.h> 49*840b91ccSNathan Whitehorn 50*840b91ccSNathan Whitehorn /* 51*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 52*840b91ccSNathan Whitehorn Software IEC/IEEE floating-point types. 53*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 54*840b91ccSNathan Whitehorn */ 55*840b91ccSNathan Whitehorn typedef unsigned int float32; 56*840b91ccSNathan Whitehorn typedef unsigned long long float64; 57*840b91ccSNathan Whitehorn #ifdef FLOATX80 58*840b91ccSNathan Whitehorn typedef struct { 59*840b91ccSNathan Whitehorn unsigned short high; 60*840b91ccSNathan Whitehorn unsigned long long low; 61*840b91ccSNathan Whitehorn } floatx80; 62*840b91ccSNathan Whitehorn #endif 63*840b91ccSNathan Whitehorn #ifdef FLOAT128 64*840b91ccSNathan Whitehorn typedef struct { 65*840b91ccSNathan Whitehorn unsigned long long high, low; 66*840b91ccSNathan Whitehorn } float128; 67*840b91ccSNathan Whitehorn #endif 68*840b91ccSNathan Whitehorn 69*840b91ccSNathan Whitehorn /* 70*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 71*840b91ccSNathan Whitehorn Software IEC/IEEE floating-point underflow tininess-detection mode. 72*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 73*840b91ccSNathan Whitehorn */ 74*840b91ccSNathan Whitehorn #ifndef SOFTFLOAT_FOR_GCC 75*840b91ccSNathan Whitehorn extern int8 float_detect_tininess; 76*840b91ccSNathan Whitehorn #endif 77*840b91ccSNathan Whitehorn enum { 78*840b91ccSNathan Whitehorn float_tininess_after_rounding = 0, 79*840b91ccSNathan Whitehorn float_tininess_before_rounding = 1 80*840b91ccSNathan Whitehorn }; 81*840b91ccSNathan Whitehorn 82*840b91ccSNathan Whitehorn /* 83*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 84*840b91ccSNathan Whitehorn Software IEC/IEEE floating-point rounding mode. 85*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 86*840b91ccSNathan Whitehorn */ 87*840b91ccSNathan Whitehorn extern fp_rnd_t float_rounding_mode; 88*840b91ccSNathan Whitehorn enum { 89*840b91ccSNathan Whitehorn float_round_nearest_even = FP_RN, 90*840b91ccSNathan Whitehorn float_round_to_zero = FP_RZ, 91*840b91ccSNathan Whitehorn float_round_down = FP_RM, 92*840b91ccSNathan Whitehorn float_round_up = FP_RP 93*840b91ccSNathan Whitehorn }; 94*840b91ccSNathan Whitehorn 95*840b91ccSNathan Whitehorn /* 96*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 97*840b91ccSNathan Whitehorn Software IEC/IEEE floating-point exception flags. 98*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 99*840b91ccSNathan Whitehorn */ 100*840b91ccSNathan Whitehorn typedef fp_except_t fp_except; 101*840b91ccSNathan Whitehorn 102*840b91ccSNathan Whitehorn extern fp_except float_exception_flags; 103*840b91ccSNathan Whitehorn extern fp_except float_exception_mask; 104*840b91ccSNathan Whitehorn enum { 105*840b91ccSNathan Whitehorn float_flag_inexact = FP_X_IMP, 106*840b91ccSNathan Whitehorn float_flag_underflow = FP_X_UFL, 107*840b91ccSNathan Whitehorn float_flag_overflow = FP_X_OFL, 108*840b91ccSNathan Whitehorn float_flag_divbyzero = FP_X_DZ, 109*840b91ccSNathan Whitehorn float_flag_invalid = FP_X_INV 110*840b91ccSNathan Whitehorn }; 111*840b91ccSNathan Whitehorn 112*840b91ccSNathan Whitehorn /* 113*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 114*840b91ccSNathan Whitehorn Routine to raise any or all of the software IEC/IEEE floating-point 115*840b91ccSNathan Whitehorn exception flags. 116*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 117*840b91ccSNathan Whitehorn */ 118*840b91ccSNathan Whitehorn void float_raise( fp_except ); 119*840b91ccSNathan Whitehorn 120*840b91ccSNathan Whitehorn /* 121*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 122*840b91ccSNathan Whitehorn Software IEC/IEEE integer-to-floating-point conversion routines. 123*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 124*840b91ccSNathan Whitehorn */ 125*840b91ccSNathan Whitehorn float32 int32_to_float32( int ); 126*840b91ccSNathan Whitehorn float64 int32_to_float64( int ); 127*840b91ccSNathan Whitehorn #ifdef FLOATX80 128*840b91ccSNathan Whitehorn floatx80 int32_to_floatx80( int ); 129*840b91ccSNathan Whitehorn #endif 130*840b91ccSNathan Whitehorn #ifdef FLOAT128 131*840b91ccSNathan Whitehorn float128 int32_to_float128( int ); 132*840b91ccSNathan Whitehorn #endif 133*840b91ccSNathan Whitehorn float32 int64_to_float32( long long ); 134*840b91ccSNathan Whitehorn float64 int64_to_float64( long long ); 135*840b91ccSNathan Whitehorn #ifdef FLOATX80 136*840b91ccSNathan Whitehorn floatx80 int64_to_floatx80( long long ); 137*840b91ccSNathan Whitehorn #endif 138*840b91ccSNathan Whitehorn #ifdef FLOAT128 139*840b91ccSNathan Whitehorn float128 int64_to_float128( long long ); 140*840b91ccSNathan Whitehorn #endif 141*840b91ccSNathan Whitehorn 142*840b91ccSNathan Whitehorn /* 143*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 144*840b91ccSNathan Whitehorn Software IEC/IEEE single-precision conversion routines. 145*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 146*840b91ccSNathan Whitehorn */ 147*840b91ccSNathan Whitehorn int float32_to_int32( float32 ); 148*840b91ccSNathan Whitehorn int float32_to_int32_round_to_zero( float32 ); 149*840b91ccSNathan Whitehorn unsigned int float32_to_uint32_round_to_zero( float32 ); 150*840b91ccSNathan Whitehorn long long float32_to_int64( float32 ); 151*840b91ccSNathan Whitehorn long long float32_to_int64_round_to_zero( float32 ); 152*840b91ccSNathan Whitehorn float64 float32_to_float64( float32 ); 153*840b91ccSNathan Whitehorn #ifdef FLOATX80 154*840b91ccSNathan Whitehorn floatx80 float32_to_floatx80( float32 ); 155*840b91ccSNathan Whitehorn #endif 156*840b91ccSNathan Whitehorn #ifdef FLOAT128 157*840b91ccSNathan Whitehorn float128 float32_to_float128( float32 ); 158*840b91ccSNathan Whitehorn #endif 159*840b91ccSNathan Whitehorn 160*840b91ccSNathan Whitehorn /* 161*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 162*840b91ccSNathan Whitehorn Software IEC/IEEE single-precision operations. 163*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 164*840b91ccSNathan Whitehorn */ 165*840b91ccSNathan Whitehorn float32 float32_round_to_int( float32 ); 166*840b91ccSNathan Whitehorn float32 float32_add( float32, float32 ); 167*840b91ccSNathan Whitehorn float32 float32_sub( float32, float32 ); 168*840b91ccSNathan Whitehorn float32 float32_mul( float32, float32 ); 169*840b91ccSNathan Whitehorn float32 float32_div( float32, float32 ); 170*840b91ccSNathan Whitehorn float32 float32_rem( float32, float32 ); 171*840b91ccSNathan Whitehorn float32 float32_sqrt( float32 ); 172*840b91ccSNathan Whitehorn int float32_eq( float32, float32 ); 173*840b91ccSNathan Whitehorn int float32_le( float32, float32 ); 174*840b91ccSNathan Whitehorn int float32_lt( float32, float32 ); 175*840b91ccSNathan Whitehorn int float32_eq_signaling( float32, float32 ); 176*840b91ccSNathan Whitehorn int float32_le_quiet( float32, float32 ); 177*840b91ccSNathan Whitehorn int float32_lt_quiet( float32, float32 ); 178*840b91ccSNathan Whitehorn #ifndef SOFTFLOAT_FOR_GCC 179*840b91ccSNathan Whitehorn int float32_is_signaling_nan( float32 ); 180*840b91ccSNathan Whitehorn #endif 181*840b91ccSNathan Whitehorn 182*840b91ccSNathan Whitehorn /* 183*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 184*840b91ccSNathan Whitehorn Software IEC/IEEE double-precision conversion routines. 185*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 186*840b91ccSNathan Whitehorn */ 187*840b91ccSNathan Whitehorn int float64_to_int32( float64 ); 188*840b91ccSNathan Whitehorn int float64_to_int32_round_to_zero( float64 ); 189*840b91ccSNathan Whitehorn unsigned int float64_to_uint32_round_to_zero( float64 ); 190*840b91ccSNathan Whitehorn long long float64_to_int64( float64 ); 191*840b91ccSNathan Whitehorn long long float64_to_int64_round_to_zero( float64 ); 192*840b91ccSNathan Whitehorn float32 float64_to_float32( float64 ); 193*840b91ccSNathan Whitehorn #ifdef FLOATX80 194*840b91ccSNathan Whitehorn floatx80 float64_to_floatx80( float64 ); 195*840b91ccSNathan Whitehorn #endif 196*840b91ccSNathan Whitehorn #ifdef FLOAT128 197*840b91ccSNathan Whitehorn float128 float64_to_float128( float64 ); 198*840b91ccSNathan Whitehorn #endif 199*840b91ccSNathan Whitehorn 200*840b91ccSNathan Whitehorn /* 201*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 202*840b91ccSNathan Whitehorn Software IEC/IEEE double-precision operations. 203*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 204*840b91ccSNathan Whitehorn */ 205*840b91ccSNathan Whitehorn float64 float64_round_to_int( float64 ); 206*840b91ccSNathan Whitehorn float64 float64_add( float64, float64 ); 207*840b91ccSNathan Whitehorn float64 float64_sub( float64, float64 ); 208*840b91ccSNathan Whitehorn float64 float64_mul( float64, float64 ); 209*840b91ccSNathan Whitehorn float64 float64_div( float64, float64 ); 210*840b91ccSNathan Whitehorn float64 float64_rem( float64, float64 ); 211*840b91ccSNathan Whitehorn float64 float64_sqrt( float64 ); 212*840b91ccSNathan Whitehorn int float64_eq( float64, float64 ); 213*840b91ccSNathan Whitehorn int float64_le( float64, float64 ); 214*840b91ccSNathan Whitehorn int float64_lt( float64, float64 ); 215*840b91ccSNathan Whitehorn int float64_eq_signaling( float64, float64 ); 216*840b91ccSNathan Whitehorn int float64_le_quiet( float64, float64 ); 217*840b91ccSNathan Whitehorn int float64_lt_quiet( float64, float64 ); 218*840b91ccSNathan Whitehorn #ifndef SOFTFLOAT_FOR_GCC 219*840b91ccSNathan Whitehorn int float64_is_signaling_nan( float64 ); 220*840b91ccSNathan Whitehorn #endif 221*840b91ccSNathan Whitehorn 222*840b91ccSNathan Whitehorn #ifdef FLOATX80 223*840b91ccSNathan Whitehorn 224*840b91ccSNathan Whitehorn /* 225*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 226*840b91ccSNathan Whitehorn Software IEC/IEEE extended double-precision conversion routines. 227*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 228*840b91ccSNathan Whitehorn */ 229*840b91ccSNathan Whitehorn int floatx80_to_int32( floatx80 ); 230*840b91ccSNathan Whitehorn int floatx80_to_int32_round_to_zero( floatx80 ); 231*840b91ccSNathan Whitehorn long long floatx80_to_int64( floatx80 ); 232*840b91ccSNathan Whitehorn long long floatx80_to_int64_round_to_zero( floatx80 ); 233*840b91ccSNathan Whitehorn float32 floatx80_to_float32( floatx80 ); 234*840b91ccSNathan Whitehorn float64 floatx80_to_float64( floatx80 ); 235*840b91ccSNathan Whitehorn #ifdef FLOAT128 236*840b91ccSNathan Whitehorn float128 floatx80_to_float128( floatx80 ); 237*840b91ccSNathan Whitehorn #endif 238*840b91ccSNathan Whitehorn 239*840b91ccSNathan Whitehorn /* 240*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 241*840b91ccSNathan Whitehorn Software IEC/IEEE extended double-precision rounding precision. Valid 242*840b91ccSNathan Whitehorn values are 32, 64, and 80. 243*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 244*840b91ccSNathan Whitehorn */ 245*840b91ccSNathan Whitehorn extern int floatx80_rounding_precision; 246*840b91ccSNathan Whitehorn 247*840b91ccSNathan Whitehorn /* 248*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 249*840b91ccSNathan Whitehorn Software IEC/IEEE extended double-precision operations. 250*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 251*840b91ccSNathan Whitehorn */ 252*840b91ccSNathan Whitehorn floatx80 floatx80_round_to_int( floatx80 ); 253*840b91ccSNathan Whitehorn floatx80 floatx80_add( floatx80, floatx80 ); 254*840b91ccSNathan Whitehorn floatx80 floatx80_sub( floatx80, floatx80 ); 255*840b91ccSNathan Whitehorn floatx80 floatx80_mul( floatx80, floatx80 ); 256*840b91ccSNathan Whitehorn floatx80 floatx80_div( floatx80, floatx80 ); 257*840b91ccSNathan Whitehorn floatx80 floatx80_rem( floatx80, floatx80 ); 258*840b91ccSNathan Whitehorn floatx80 floatx80_sqrt( floatx80 ); 259*840b91ccSNathan Whitehorn int floatx80_eq( floatx80, floatx80 ); 260*840b91ccSNathan Whitehorn int floatx80_le( floatx80, floatx80 ); 261*840b91ccSNathan Whitehorn int floatx80_lt( floatx80, floatx80 ); 262*840b91ccSNathan Whitehorn int floatx80_eq_signaling( floatx80, floatx80 ); 263*840b91ccSNathan Whitehorn int floatx80_le_quiet( floatx80, floatx80 ); 264*840b91ccSNathan Whitehorn int floatx80_lt_quiet( floatx80, floatx80 ); 265*840b91ccSNathan Whitehorn int floatx80_is_signaling_nan( floatx80 ); 266*840b91ccSNathan Whitehorn 267*840b91ccSNathan Whitehorn #endif 268*840b91ccSNathan Whitehorn 269*840b91ccSNathan Whitehorn #ifdef FLOAT128 270*840b91ccSNathan Whitehorn 271*840b91ccSNathan Whitehorn /* 272*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 273*840b91ccSNathan Whitehorn Software IEC/IEEE quadruple-precision conversion routines. 274*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 275*840b91ccSNathan Whitehorn */ 276*840b91ccSNathan Whitehorn int float128_to_int32( float128 ); 277*840b91ccSNathan Whitehorn int float128_to_int32_round_to_zero( float128 ); 278*840b91ccSNathan Whitehorn long long float128_to_int64( float128 ); 279*840b91ccSNathan Whitehorn long long float128_to_int64_round_to_zero( float128 ); 280*840b91ccSNathan Whitehorn float32 float128_to_float32( float128 ); 281*840b91ccSNathan Whitehorn float64 float128_to_float64( float128 ); 282*840b91ccSNathan Whitehorn #ifdef FLOATX80 283*840b91ccSNathan Whitehorn floatx80 float128_to_floatx80( float128 ); 284*840b91ccSNathan Whitehorn #endif 285*840b91ccSNathan Whitehorn 286*840b91ccSNathan Whitehorn /* 287*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 288*840b91ccSNathan Whitehorn Software IEC/IEEE quadruple-precision operations. 289*840b91ccSNathan Whitehorn ------------------------------------------------------------------------------- 290*840b91ccSNathan Whitehorn */ 291*840b91ccSNathan Whitehorn float128 float128_round_to_int( float128 ); 292*840b91ccSNathan Whitehorn float128 float128_add( float128, float128 ); 293*840b91ccSNathan Whitehorn float128 float128_sub( float128, float128 ); 294*840b91ccSNathan Whitehorn float128 float128_mul( float128, float128 ); 295*840b91ccSNathan Whitehorn float128 float128_div( float128, float128 ); 296*840b91ccSNathan Whitehorn float128 float128_rem( float128, float128 ); 297*840b91ccSNathan Whitehorn float128 float128_sqrt( float128 ); 298*840b91ccSNathan Whitehorn int float128_eq( float128, float128 ); 299*840b91ccSNathan Whitehorn int float128_le( float128, float128 ); 300*840b91ccSNathan Whitehorn int float128_lt( float128, float128 ); 301*840b91ccSNathan Whitehorn int float128_eq_signaling( float128, float128 ); 302*840b91ccSNathan Whitehorn int float128_le_quiet( float128, float128 ); 303*840b91ccSNathan Whitehorn int float128_lt_quiet( float128, float128 ); 304*840b91ccSNathan Whitehorn int float128_is_signaling_nan( float128 ); 305*840b91ccSNathan Whitehorn 306*840b91ccSNathan Whitehorn #endif 307*840b91ccSNathan Whitehorn 308