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