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