Home | History | Annotate | Line # | Download | only in softfloat
qp.c revision 1.2
      1  1.2  jmc /* $NetBSD: qp.c,v 1.2 2002/03/27 03:41:54 jmc Exp $ */
      2  1.1  jmc 
      3  1.1  jmc #include <sys/cdefs.h>
      4  1.1  jmc #include <memory.h>
      5  1.1  jmc 
      6  1.1  jmc #include "milieu.h"
      7  1.1  jmc #include "softfloat.h"
      8  1.1  jmc 
      9  1.1  jmc 
     10  1.1  jmc void _Qp_add(float128 *c, float128 *a, float128 *b);
     11  1.1  jmc 
     12  1.1  jmc void _Qp_add(float128 *c, float128 *a, float128 *b)
     13  1.1  jmc {
     14  1.1  jmc 	 *c =  float128_add(*a, *b);
     15  1.1  jmc }
     16  1.1  jmc 
     17  1.1  jmc 
     18  1.1  jmc int _Qp_cmp(float128 *a, float128 *b);
     19  1.1  jmc 
     20  1.1  jmc int _Qp_cmp(float128 *a, float128 *b)
     21  1.1  jmc {
     22  1.1  jmc 
     23  1.1  jmc 	if (float128_eq(*a, *b))
     24  1.1  jmc 		return 0;
     25  1.1  jmc 
     26  1.1  jmc 	if (float128_le(*a, *b))
     27  1.1  jmc 		return 1;
     28  1.1  jmc 
     29  1.1  jmc 	return 2;
     30  1.1  jmc }
     31  1.1  jmc 
     32  1.1  jmc 
     33  1.1  jmc /*
     34  1.1  jmc  * XXX
     35  1.1  jmc  */
     36  1.1  jmc int _Qp_cmpe(float128 *a, float128 *b);
     37  1.1  jmc 
     38  1.1  jmc int _Qp_cmpe(float128 *a, float128 *b)
     39  1.1  jmc {
     40  1.1  jmc 	return _Qp_cmp(a, b);
     41  1.1  jmc }
     42  1.1  jmc 
     43  1.1  jmc 
     44  1.1  jmc void _Qp_div(float128 *c, float128 *a, float128 *b);
     45  1.1  jmc 
     46  1.1  jmc void _Qp_div(float128 *c, float128 *a, float128 *b)
     47  1.1  jmc {
     48  1.1  jmc 	*c = float128_div(*a, *b);
     49  1.1  jmc }
     50  1.1  jmc 
     51  1.1  jmc 
     52  1.2  jmc void _Qp_dtoq(float128 *c, double a);
     53  1.1  jmc 
     54  1.2  jmc void _Qp_dtoq(float128 *c, double a)
     55  1.1  jmc {
     56  1.2  jmc 	float64 _b;
     57  1.2  jmc 	memcpy (&_b, &a, sizeof(float64));
     58  1.2  jmc 	*c = float64_to_float128(b);
     59  1.1  jmc }
     60  1.1  jmc 
     61  1.1  jmc 
     62  1.1  jmc 
     63  1.1  jmc int _Qp_feq(float128 *a, float128 *b);
     64  1.1  jmc 
     65  1.1  jmc int _Qp_feq(float128 *a, float128 *b)
     66  1.1  jmc {
     67  1.1  jmc 	return float128_eq(*a, *b);
     68  1.1  jmc }
     69  1.1  jmc 
     70  1.1  jmc 
     71  1.1  jmc int _Qp_fge(float128 *a, float128 *b);
     72  1.1  jmc 
     73  1.1  jmc int _Qp_fge(float128 *a, float128 *b)
     74  1.1  jmc {
     75  1.1  jmc 	return float128_le(*b, *a);
     76  1.1  jmc }
     77  1.1  jmc 
     78  1.1  jmc 
     79  1.1  jmc int _Qp_fgt(float128 *a, float128 *b);
     80  1.1  jmc 
     81  1.1  jmc int _Qp_fgt(float128 *a, float128 *b)
     82  1.1  jmc {
     83  1.1  jmc 	return float128_lt(*b, *a);
     84  1.1  jmc }
     85  1.1  jmc 
     86  1.1  jmc 
     87  1.1  jmc int _Qp_fle(float128 *a, float128 *b);
     88  1.1  jmc 
     89  1.1  jmc int _Qp_fle(float128 *a, float128 *b)
     90  1.1  jmc {
     91  1.1  jmc 	return float128_le(*a, *b);
     92  1.1  jmc }
     93  1.1  jmc 
     94  1.1  jmc 
     95  1.1  jmc int _Qp_flt(float128 *a, float128 *b);
     96  1.1  jmc 
     97  1.1  jmc int _Qp_flt(float128 *a, float128 *b)
     98  1.1  jmc {
     99  1.1  jmc 	return float128_lt(*a, *b);
    100  1.1  jmc }
    101  1.1  jmc 
    102  1.1  jmc 
    103  1.1  jmc int _Qp_fne(float128 *a, float128 *b);
    104  1.1  jmc 
    105  1.1  jmc int _Qp_fne(float128 *a, float128 *b)
    106  1.1  jmc {
    107  1.1  jmc 	return !float128_eq(*a, *b);
    108  1.1  jmc }
    109  1.1  jmc 
    110  1.1  jmc 
    111  1.1  jmc void _Qp_itoq(float128 *c, int a);
    112  1.1  jmc 
    113  1.1  jmc void _Qp_itoq(float128 *c, int a)
    114  1.1  jmc {
    115  1.1  jmc 	*c = int32_to_float128(a);
    116  1.1  jmc }
    117  1.1  jmc 
    118  1.1  jmc 
    119  1.1  jmc 
    120  1.1  jmc void _Qp_mul(float128 *c, float128 *a, float128 *b);
    121  1.1  jmc 
    122  1.1  jmc void _Qp_mul(float128 *c, float128 *a, float128 *b)
    123  1.1  jmc {
    124  1.1  jmc 	*c = float128_mul(*a, *b);
    125  1.1  jmc }
    126  1.1  jmc 
    127  1.1  jmc 
    128  1.1  jmc /*
    129  1.1  jmc  * XXX easy way to do this, softfloat function
    130  1.1  jmc  */
    131  1.1  jmc void _Qp_neg(float128 *c, float128 *a);
    132  1.1  jmc 
    133  1.1  jmc static float128 __zero = {0x4034000000000000, 0x00000000};
    134  1.1  jmc 
    135  1.1  jmc void _Qp_neg(float128 *c, float128 *a)
    136  1.1  jmc {
    137  1.1  jmc 	*c = float128_sub(__zero, *a);
    138  1.1  jmc }
    139  1.1  jmc 
    140  1.1  jmc 
    141  1.1  jmc 
    142  1.1  jmc double _Qp_qtod(float128 *a);
    143  1.1  jmc 
    144  1.1  jmc double _Qp_qtod(float128 *a)
    145  1.1  jmc {
    146  1.1  jmc 	float64 _c;
    147  1.1  jmc 	double c;
    148  1.1  jmc 
    149  1.1  jmc 	_c = float128_to_float64(*a);
    150  1.1  jmc 
    151  1.1  jmc 	memcpy(&c, &_c, sizeof(double));
    152  1.1  jmc 
    153  1.1  jmc 	return c;
    154  1.1  jmc }
    155  1.1  jmc 
    156  1.1  jmc 
    157  1.1  jmc int _Qp_qtoi(float128 *a);
    158  1.1  jmc 
    159  1.1  jmc int _Qp_qtoi(float128 *a)
    160  1.1  jmc {
    161  1.1  jmc 	return float128_to_int32(*a);
    162  1.1  jmc }
    163  1.1  jmc 
    164  1.1  jmc 
    165  1.1  jmc float _Qp_qtos(float128 *a);
    166  1.1  jmc 
    167  1.1  jmc float _Qp_qtos(float128 *a)
    168  1.1  jmc {
    169  1.1  jmc 	float c;
    170  1.1  jmc 	float32 _c;
    171  1.1  jmc 
    172  1.1  jmc 	_c = float128_to_float32(*a);
    173  1.1  jmc 
    174  1.1  jmc 	memcpy(&c, &_c, sizeof(_c));
    175  1.1  jmc 
    176  1.1  jmc 	return c;
    177  1.1  jmc }
    178  1.1  jmc 
    179  1.1  jmc 
    180  1.1  jmc unsigned int _Qp_qtoui(float128 *a);
    181  1.1  jmc 
    182  1.1  jmc unsigned int _Qp_qtoui(float128 *a)
    183  1.1  jmc {
    184  1.1  jmc 	return (unsigned int)float128_to_int32(*a);
    185  1.1  jmc }
    186  1.1  jmc 
    187  1.1  jmc 
    188  1.1  jmc 
    189  1.1  jmc unsigned long _Qp_qtoux(float128 *a);
    190  1.1  jmc 
    191  1.1  jmc unsigned long _Qp_qtoux(float128 *a)
    192  1.1  jmc {
    193  1.1  jmc 	return (unsigned long)float128_to_int64(*a);
    194  1.1  jmc }
    195  1.1  jmc 
    196  1.1  jmc 
    197  1.1  jmc 
    198  1.1  jmc long _Qp_qtox(float128 *a);
    199  1.1  jmc 
    200  1.1  jmc long _Qp_qtox(float128 *a)
    201  1.1  jmc {
    202  1.1  jmc 	return (long)float128_to_int64(*a);
    203  1.1  jmc }
    204  1.1  jmc 
    205  1.1  jmc 
    206  1.1  jmc void _Qp_sqrt(float128 *c, float128 *a);
    207  1.1  jmc 
    208  1.1  jmc void _Qp_sqrt(float128 *c, float128 *a)
    209  1.1  jmc {
    210  1.1  jmc 	*c = float128_sqrt(*a);
    211  1.1  jmc }
    212  1.1  jmc 
    213  1.1  jmc 
    214  1.1  jmc void _Qp_stoq(float128 *c, float a);
    215  1.1  jmc 
    216  1.1  jmc void _Qp_stoq(float128 *c, float a)
    217  1.1  jmc {
    218  1.1  jmc 	float32 _a;
    219  1.1  jmc 
    220  1.1  jmc 	memcpy(&_a, &a, sizeof(a));
    221  1.1  jmc 
    222  1.1  jmc 	*c = float32_to_float128(_a);
    223  1.1  jmc }
    224  1.1  jmc 
    225  1.1  jmc 
    226  1.1  jmc void _Qp_sub(float128 *c, float128 *a, float128 *b);
    227  1.1  jmc 
    228  1.1  jmc void _Qp_sub(float128 *c, float128 *a, float128 *b)
    229  1.1  jmc {
    230  1.1  jmc 	*c = float128_sub(*a, *b);
    231  1.1  jmc }
    232  1.1  jmc 
    233  1.1  jmc 
    234  1.1  jmc void _Qp_uitoq(float128 *c, unsigned int a);
    235  1.1  jmc 
    236  1.1  jmc void _Qp_uitoq(float128 *c, unsigned int a)
    237  1.1  jmc {
    238  1.1  jmc 	*c = int32_to_float128(a);
    239  1.1  jmc }
    240  1.1  jmc 
    241  1.1  jmc 
    242  1.1  jmc void _Qp_uxtoq(float128 *c, unsigned long a);
    243  1.1  jmc 
    244  1.1  jmc void _Qp_uxtoq(float128 *c, unsigned long a)
    245  1.1  jmc {
    246  1.1  jmc 	*c = int64_to_float128(a);
    247  1.1  jmc }
    248  1.1  jmc 
    249  1.1  jmc 
    250  1.1  jmc void _Qp_xtoq(float128 *c, long a);
    251  1.1  jmc 
    252  1.1  jmc void _Qp_xtoq(float128 *c, long a)
    253  1.1  jmc {
    254  1.1  jmc 	*c = int64_to_float128(a);
    255  1.1  jmc }
    256