qp.c revision 1.3 1 1.3 jmc /* $NetBSD: qp.c,v 1.3 2002/03/27 03:43:39 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.3 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