Home | History | Annotate | Line # | Download | only in softfloat
qp.c revision 1.1
      1  1.1  jmc /* $NetBSD: qp.c,v 1.1 2002/02/05 07:53:05 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.1  jmc void _Qp_dtoq(float128 *c, float64 *a);
     53  1.1  jmc 
     54  1.1  jmc void _Qp_dtoq(float128 *c, float64 *a)
     55  1.1  jmc {
     56  1.1  jmc 	*c = float64_to_float128(*a);
     57  1.1  jmc }
     58  1.1  jmc 
     59  1.1  jmc 
     60  1.1  jmc 
     61  1.1  jmc int _Qp_feq(float128 *a, float128 *b);
     62  1.1  jmc 
     63  1.1  jmc int _Qp_feq(float128 *a, float128 *b)
     64  1.1  jmc {
     65  1.1  jmc 	return float128_eq(*a, *b);
     66  1.1  jmc }
     67  1.1  jmc 
     68  1.1  jmc 
     69  1.1  jmc int _Qp_fge(float128 *a, float128 *b);
     70  1.1  jmc 
     71  1.1  jmc int _Qp_fge(float128 *a, float128 *b)
     72  1.1  jmc {
     73  1.1  jmc 	return float128_le(*b, *a);
     74  1.1  jmc }
     75  1.1  jmc 
     76  1.1  jmc 
     77  1.1  jmc int _Qp_fgt(float128 *a, float128 *b);
     78  1.1  jmc 
     79  1.1  jmc int _Qp_fgt(float128 *a, float128 *b)
     80  1.1  jmc {
     81  1.1  jmc 	return float128_lt(*b, *a);
     82  1.1  jmc }
     83  1.1  jmc 
     84  1.1  jmc 
     85  1.1  jmc int _Qp_fle(float128 *a, float128 *b);
     86  1.1  jmc 
     87  1.1  jmc int _Qp_fle(float128 *a, float128 *b)
     88  1.1  jmc {
     89  1.1  jmc 	return float128_le(*a, *b);
     90  1.1  jmc }
     91  1.1  jmc 
     92  1.1  jmc 
     93  1.1  jmc int _Qp_flt(float128 *a, float128 *b);
     94  1.1  jmc 
     95  1.1  jmc int _Qp_flt(float128 *a, float128 *b)
     96  1.1  jmc {
     97  1.1  jmc 	return float128_lt(*a, *b);
     98  1.1  jmc }
     99  1.1  jmc 
    100  1.1  jmc 
    101  1.1  jmc int _Qp_fne(float128 *a, float128 *b);
    102  1.1  jmc 
    103  1.1  jmc int _Qp_fne(float128 *a, float128 *b)
    104  1.1  jmc {
    105  1.1  jmc 	return !float128_eq(*a, *b);
    106  1.1  jmc }
    107  1.1  jmc 
    108  1.1  jmc 
    109  1.1  jmc void _Qp_itoq(float128 *c, int a);
    110  1.1  jmc 
    111  1.1  jmc void _Qp_itoq(float128 *c, int a)
    112  1.1  jmc {
    113  1.1  jmc 	*c = int32_to_float128(a);
    114  1.1  jmc }
    115  1.1  jmc 
    116  1.1  jmc 
    117  1.1  jmc 
    118  1.1  jmc void _Qp_mul(float128 *c, float128 *a, float128 *b);
    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 	*c = float128_mul(*a, *b);
    123  1.1  jmc }
    124  1.1  jmc 
    125  1.1  jmc 
    126  1.1  jmc /*
    127  1.1  jmc  * XXX easy way to do this, softfloat function
    128  1.1  jmc  */
    129  1.1  jmc void _Qp_neg(float128 *c, float128 *a);
    130  1.1  jmc 
    131  1.1  jmc static float128 __zero = {0x4034000000000000, 0x00000000};
    132  1.1  jmc 
    133  1.1  jmc void _Qp_neg(float128 *c, float128 *a)
    134  1.1  jmc {
    135  1.1  jmc 	*c = float128_sub(__zero, *a);
    136  1.1  jmc }
    137  1.1  jmc 
    138  1.1  jmc 
    139  1.1  jmc 
    140  1.1  jmc double _Qp_qtod(float128 *a);
    141  1.1  jmc 
    142  1.1  jmc double _Qp_qtod(float128 *a)
    143  1.1  jmc {
    144  1.1  jmc 	float64 _c;
    145  1.1  jmc 	double c;
    146  1.1  jmc 
    147  1.1  jmc 	_c = float128_to_float64(*a);
    148  1.1  jmc 
    149  1.1  jmc 	memcpy(&c, &_c, sizeof(double));
    150  1.1  jmc 
    151  1.1  jmc 	return c;
    152  1.1  jmc }
    153  1.1  jmc 
    154  1.1  jmc 
    155  1.1  jmc int _Qp_qtoi(float128 *a);
    156  1.1  jmc 
    157  1.1  jmc int _Qp_qtoi(float128 *a)
    158  1.1  jmc {
    159  1.1  jmc 	return float128_to_int32(*a);
    160  1.1  jmc }
    161  1.1  jmc 
    162  1.1  jmc 
    163  1.1  jmc float _Qp_qtos(float128 *a);
    164  1.1  jmc 
    165  1.1  jmc float _Qp_qtos(float128 *a)
    166  1.1  jmc {
    167  1.1  jmc 	float c;
    168  1.1  jmc 	float32 _c;
    169  1.1  jmc 
    170  1.1  jmc 	_c = float128_to_float32(*a);
    171  1.1  jmc 
    172  1.1  jmc 	memcpy(&c, &_c, sizeof(_c));
    173  1.1  jmc 
    174  1.1  jmc 	return c;
    175  1.1  jmc }
    176  1.1  jmc 
    177  1.1  jmc 
    178  1.1  jmc unsigned int _Qp_qtoui(float128 *a);
    179  1.1  jmc 
    180  1.1  jmc unsigned int _Qp_qtoui(float128 *a)
    181  1.1  jmc {
    182  1.1  jmc 	return (unsigned int)float128_to_int32(*a);
    183  1.1  jmc }
    184  1.1  jmc 
    185  1.1  jmc 
    186  1.1  jmc 
    187  1.1  jmc unsigned long _Qp_qtoux(float128 *a);
    188  1.1  jmc 
    189  1.1  jmc unsigned long _Qp_qtoux(float128 *a)
    190  1.1  jmc {
    191  1.1  jmc 	return (unsigned long)float128_to_int64(*a);
    192  1.1  jmc }
    193  1.1  jmc 
    194  1.1  jmc 
    195  1.1  jmc 
    196  1.1  jmc long _Qp_qtox(float128 *a);
    197  1.1  jmc 
    198  1.1  jmc long _Qp_qtox(float128 *a)
    199  1.1  jmc {
    200  1.1  jmc 	return (long)float128_to_int64(*a);
    201  1.1  jmc }
    202  1.1  jmc 
    203  1.1  jmc 
    204  1.1  jmc void _Qp_sqrt(float128 *c, float128 *a);
    205  1.1  jmc 
    206  1.1  jmc void _Qp_sqrt(float128 *c, float128 *a)
    207  1.1  jmc {
    208  1.1  jmc 	*c = float128_sqrt(*a);
    209  1.1  jmc }
    210  1.1  jmc 
    211  1.1  jmc 
    212  1.1  jmc void _Qp_stoq(float128 *c, float a);
    213  1.1  jmc 
    214  1.1  jmc void _Qp_stoq(float128 *c, float a)
    215  1.1  jmc {
    216  1.1  jmc 	float32 _a;
    217  1.1  jmc 
    218  1.1  jmc 	memcpy(&_a, &a, sizeof(a));
    219  1.1  jmc 
    220  1.1  jmc 	*c = float32_to_float128(_a);
    221  1.1  jmc }
    222  1.1  jmc 
    223  1.1  jmc 
    224  1.1  jmc void _Qp_sub(float128 *c, float128 *a, float128 *b);
    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 	*c = float128_sub(*a, *b);
    229  1.1  jmc }
    230  1.1  jmc 
    231  1.1  jmc 
    232  1.1  jmc void _Qp_uitoq(float128 *c, unsigned int a);
    233  1.1  jmc 
    234  1.1  jmc void _Qp_uitoq(float128 *c, unsigned int a)
    235  1.1  jmc {
    236  1.1  jmc 	*c = int32_to_float128(a);
    237  1.1  jmc }
    238  1.1  jmc 
    239  1.1  jmc 
    240  1.1  jmc void _Qp_uxtoq(float128 *c, unsigned long a);
    241  1.1  jmc 
    242  1.1  jmc void _Qp_uxtoq(float128 *c, unsigned long a)
    243  1.1  jmc {
    244  1.1  jmc 	*c = int64_to_float128(a);
    245  1.1  jmc }
    246  1.1  jmc 
    247  1.1  jmc 
    248  1.1  jmc void _Qp_xtoq(float128 *c, long a);
    249  1.1  jmc 
    250  1.1  jmc void _Qp_xtoq(float128 *c, long a)
    251  1.1  jmc {
    252  1.1  jmc 	*c = int64_to_float128(a);
    253  1.1  jmc }
    254