]>
Commit | Line | Data |
---|---|---|
1 | /* | |
2 | * Linux/PA-RISC Project (http://www.parisc-linux.org/) | |
3 | * | |
4 | * Floating-point emulation code | |
5 | * Copyright (C) 2001 Hewlett-Packard (Paul Bame) <bame@debian.org> | |
6 | * | |
7 | * This program is free software; you can redistribute it and/or modify | |
8 | * it under the terms of the GNU General Public License as published by | |
9 | * the Free Software Foundation; either version 2, or (at your option) | |
10 | * any later version. | |
11 | * | |
12 | * This program is distributed in the hope that it will be useful, | |
13 | * but WITHOUT ANY WARRANTY; without even the implied warranty of | |
14 | * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the | |
15 | * GNU General Public License for more details. | |
16 | * | |
17 | * You should have received a copy of the GNU General Public License | |
18 | * along with this program; if not, write to the Free Software | |
19 | * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA | |
20 | */ | |
21 | /* | |
22 | * BEGIN_DESC | |
23 | * | |
24 | * File: | |
25 | * @(#) pa/spmath/fcnvfut.c $Revision: 1.1 $ | |
26 | * | |
27 | * Purpose: | |
28 | * Floating-point to Unsigned Fixed-point Converts with Truncation | |
29 | * | |
30 | * External Interfaces: | |
31 | * dbl_to_dbl_fcnvfut(srcptr,nullptr,dstptr,status) | |
32 | * dbl_to_sgl_fcnvfut(srcptr,nullptr,dstptr,status) | |
33 | * sgl_to_dbl_fcnvfut(srcptr,nullptr,dstptr,status) | |
34 | * sgl_to_sgl_fcnvfut(srcptr,nullptr,dstptr,status) | |
35 | * | |
36 | * Internal Interfaces: | |
37 | * | |
38 | * Theory: | |
39 | * <<please update with a overview of the operation of this file>> | |
40 | * | |
41 | * END_DESC | |
42 | */ | |
43 | ||
44 | ||
45 | #include "float.h" | |
46 | #include "sgl_float.h" | |
47 | #include "dbl_float.h" | |
48 | #include "cnv_float.h" | |
49 | ||
50 | /************************************************************************ | |
51 | * Floating-point to Unsigned Fixed-point Converts with Truncation * | |
52 | ************************************************************************/ | |
53 | ||
54 | /* | |
55 | * Convert single floating-point to single fixed-point format | |
56 | * with truncated result | |
57 | */ | |
58 | /*ARGSUSED*/ | |
59 | int | |
60 | sgl_to_sgl_fcnvfut (sgl_floating_point * srcptr, unsigned int *nullptr, | |
61 | unsigned int *dstptr, unsigned int *status) | |
62 | { | |
63 | register unsigned int src, result; | |
64 | register int src_exponent; | |
65 | ||
66 | src = *srcptr; | |
67 | src_exponent = Sgl_exponent(src) - SGL_BIAS; | |
68 | ||
69 | /* | |
70 | * Test for overflow | |
71 | */ | |
72 | if (src_exponent > SGL_FX_MAX_EXP + 1) { | |
73 | if (Sgl_isone_sign(src)) { | |
74 | result = 0; | |
75 | } else { | |
76 | result = 0xffffffff; | |
77 | } | |
78 | if (Is_invalidtrap_enabled()) { | |
79 | return(INVALIDEXCEPTION); | |
80 | } | |
81 | Set_invalidflag(); | |
82 | *dstptr = result; | |
83 | return(NOEXCEPTION); | |
84 | } | |
85 | /* | |
86 | * Generate result | |
87 | */ | |
88 | if (src_exponent >= 0) { | |
89 | /* | |
90 | * Check sign. | |
91 | * If negative, trap unimplemented. | |
92 | */ | |
93 | if (Sgl_isone_sign(src)) { | |
94 | result = 0; | |
95 | if (Is_invalidtrap_enabled()) { | |
96 | return(INVALIDEXCEPTION); | |
97 | } | |
98 | Set_invalidflag(); | |
99 | *dstptr = result; | |
100 | return(NOEXCEPTION); | |
101 | } | |
102 | Sgl_clear_signexponent_set_hidden(src); | |
103 | Suint_from_sgl_mantissa(src,src_exponent,result); | |
104 | *dstptr = result; | |
105 | ||
106 | /* check for inexact */ | |
107 | if (Sgl_isinexact_to_unsigned(src,src_exponent)) { | |
108 | if (Is_inexacttrap_enabled()) return(INEXACTEXCEPTION); | |
109 | else Set_inexactflag(); | |
110 | } | |
111 | } | |
112 | else { | |
113 | *dstptr = 0; | |
114 | ||
115 | /* check for inexact */ | |
116 | if (Sgl_isnotzero_exponentmantissa(src)) { | |
117 | if (Is_inexacttrap_enabled()) return(INEXACTEXCEPTION); | |
118 | else Set_inexactflag(); | |
119 | } | |
120 | } | |
121 | return(NOEXCEPTION); | |
122 | } | |
123 | ||
124 | /* | |
125 | * Single Floating-point to Double Unsigned Fixed | |
126 | */ | |
127 | /*ARGSUSED*/ | |
128 | int | |
129 | sgl_to_dbl_fcnvfut (sgl_floating_point * srcptr, unsigned int *nullptr, | |
130 | dbl_unsigned * dstptr, unsigned int *status) | |
131 | { | |
132 | register int src_exponent; | |
133 | register unsigned int src, resultp1, resultp2; | |
134 | ||
135 | src = *srcptr; | |
136 | src_exponent = Sgl_exponent(src) - SGL_BIAS; | |
137 | ||
138 | /* | |
139 | * Test for overflow | |
140 | */ | |
141 | if (src_exponent > DBL_FX_MAX_EXP + 1) { | |
142 | if (Sgl_isone_sign(src)) { | |
143 | resultp1 = resultp2 = 0; | |
144 | } else { | |
145 | resultp1 = resultp2 = 0xffffffff; | |
146 | } | |
147 | if (Is_invalidtrap_enabled()) { | |
148 | return(INVALIDEXCEPTION); | |
149 | } | |
150 | Set_invalidflag(); | |
151 | Duint_copytoptr(resultp1,resultp2,dstptr); | |
152 | return(NOEXCEPTION); | |
153 | } | |
154 | /* | |
155 | * Generate result | |
156 | */ | |
157 | if (src_exponent >= 0) { | |
158 | /* | |
159 | * Check sign. | |
160 | * If negative, trap unimplemented. | |
161 | */ | |
162 | if (Sgl_isone_sign(src)) { | |
163 | resultp1 = resultp2 = 0; | |
164 | if (Is_invalidtrap_enabled()) { | |
165 | return(INVALIDEXCEPTION); | |
166 | } | |
167 | Set_invalidflag(); | |
168 | Duint_copytoptr(resultp1,resultp2,dstptr); | |
169 | return(NOEXCEPTION); | |
170 | } | |
171 | Sgl_clear_signexponent_set_hidden(src); | |
172 | Duint_from_sgl_mantissa(src,src_exponent,resultp1,resultp2); | |
173 | Duint_copytoptr(resultp1,resultp2,dstptr); | |
174 | ||
175 | /* check for inexact */ | |
176 | if (Sgl_isinexact_to_unsigned(src,src_exponent)) { | |
177 | if (Is_inexacttrap_enabled()) return(INEXACTEXCEPTION); | |
178 | else Set_inexactflag(); | |
179 | } | |
180 | } | |
181 | else { | |
182 | Duint_setzero(resultp1,resultp2); | |
183 | Duint_copytoptr(resultp1,resultp2,dstptr); | |
184 | ||
185 | /* check for inexact */ | |
186 | if (Sgl_isnotzero_exponentmantissa(src)) { | |
187 | if (Is_inexacttrap_enabled()) return(INEXACTEXCEPTION); | |
188 | else Set_inexactflag(); | |
189 | } | |
190 | } | |
191 | return(NOEXCEPTION); | |
192 | } | |
193 | ||
194 | /* | |
195 | * Double Floating-point to Single Unsigned Fixed | |
196 | */ | |
197 | /*ARGSUSED*/ | |
198 | int | |
199 | dbl_to_sgl_fcnvfut (dbl_floating_point * srcptr, unsigned int *nullptr, | |
200 | unsigned int *dstptr, unsigned int *status) | |
201 | { | |
202 | register unsigned int srcp1, srcp2, result; | |
203 | register int src_exponent; | |
204 | ||
205 | Dbl_copyfromptr(srcptr,srcp1,srcp2); | |
206 | src_exponent = Dbl_exponent(srcp1) - DBL_BIAS; | |
207 | ||
208 | /* | |
209 | * Test for overflow | |
210 | */ | |
211 | if (src_exponent > SGL_FX_MAX_EXP + 1) { | |
212 | if (Dbl_isone_sign(srcp1)) { | |
213 | result = 0; | |
214 | } else { | |
215 | result = 0xffffffff; | |
216 | } | |
217 | if (Is_invalidtrap_enabled()) { | |
218 | return(INVALIDEXCEPTION); | |
219 | } | |
220 | Set_invalidflag(); | |
221 | *dstptr = result; | |
222 | return(NOEXCEPTION); | |
223 | } | |
224 | /* | |
225 | * Generate result | |
226 | */ | |
227 | if (src_exponent >= 0) { | |
228 | /* | |
229 | * Check sign. | |
230 | * If negative, trap unimplemented. | |
231 | */ | |
232 | if (Dbl_isone_sign(srcp1)) { | |
233 | result = 0; | |
234 | if (Is_invalidtrap_enabled()) { | |
235 | return(INVALIDEXCEPTION); | |
236 | } | |
237 | Set_invalidflag(); | |
238 | *dstptr = result; | |
239 | return(NOEXCEPTION); | |
240 | } | |
241 | Dbl_clear_signexponent_set_hidden(srcp1); | |
242 | Suint_from_dbl_mantissa(srcp1,srcp2,src_exponent,result); | |
243 | *dstptr = result; | |
244 | ||
245 | /* check for inexact */ | |
246 | if (Dbl_isinexact_to_unsigned(srcp1,srcp2,src_exponent)) { | |
247 | if (Is_inexacttrap_enabled()) return(INEXACTEXCEPTION); | |
248 | else Set_inexactflag(); | |
249 | } | |
250 | } | |
251 | else { | |
252 | *dstptr = 0; | |
253 | ||
254 | /* check for inexact */ | |
255 | if (Dbl_isnotzero_exponentmantissa(srcp1,srcp2)) { | |
256 | if (Is_inexacttrap_enabled()) return(INEXACTEXCEPTION); | |
257 | else Set_inexactflag(); | |
258 | } | |
259 | } | |
260 | return(NOEXCEPTION); | |
261 | } | |
262 | ||
263 | /* | |
264 | * Double Floating-point to Double Unsigned Fixed | |
265 | */ | |
266 | /*ARGSUSED*/ | |
267 | int | |
268 | dbl_to_dbl_fcnvfut (dbl_floating_point * srcptr, unsigned int *nullptr, | |
269 | dbl_unsigned * dstptr, unsigned int *status) | |
270 | { | |
271 | register int src_exponent; | |
272 | register unsigned int srcp1, srcp2, resultp1, resultp2; | |
273 | ||
274 | Dbl_copyfromptr(srcptr,srcp1,srcp2); | |
275 | src_exponent = Dbl_exponent(srcp1) - DBL_BIAS; | |
276 | ||
277 | /* | |
278 | * Test for overflow | |
279 | */ | |
280 | if (src_exponent > DBL_FX_MAX_EXP + 1) { | |
281 | if (Dbl_isone_sign(srcp1)) { | |
282 | resultp1 = resultp2 = 0; | |
283 | } else { | |
284 | resultp1 = resultp2 = 0xffffffff; | |
285 | } | |
286 | if (Is_invalidtrap_enabled()) { | |
287 | return(INVALIDEXCEPTION); | |
288 | } | |
289 | Set_invalidflag(); | |
290 | Duint_copytoptr(resultp1,resultp2,dstptr); | |
291 | return(NOEXCEPTION); | |
292 | } | |
293 | /* | |
294 | * Generate result | |
295 | */ | |
296 | if (src_exponent >= 0) { | |
297 | /* | |
298 | * Check sign. | |
299 | * If negative, trap unimplemented. | |
300 | */ | |
301 | if (Dbl_isone_sign(srcp1)) { | |
302 | resultp1 = resultp2 = 0; | |
303 | if (Is_invalidtrap_enabled()) { | |
304 | return(INVALIDEXCEPTION); | |
305 | } | |
306 | Set_invalidflag(); | |
307 | Duint_copytoptr(resultp1,resultp2,dstptr); | |
308 | return(NOEXCEPTION); | |
309 | } | |
310 | Dbl_clear_signexponent_set_hidden(srcp1); | |
311 | Duint_from_dbl_mantissa(srcp1,srcp2,src_exponent, | |
312 | resultp1,resultp2); | |
313 | Duint_copytoptr(resultp1,resultp2,dstptr); | |
314 | ||
315 | /* check for inexact */ | |
316 | if (Dbl_isinexact_to_unsigned(srcp1,srcp2,src_exponent)) { | |
317 | if (Is_inexacttrap_enabled()) return(INEXACTEXCEPTION); | |
318 | else Set_inexactflag(); | |
319 | } | |
320 | } | |
321 | else { | |
322 | Duint_setzero(resultp1,resultp2); | |
323 | Duint_copytoptr(resultp1,resultp2,dstptr); | |
324 | ||
325 | /* check for inexact */ | |
326 | if (Dbl_isnotzero_exponentmantissa(srcp1,srcp2)) { | |
327 | if (Is_inexacttrap_enabled()) return(INEXACTEXCEPTION); | |
328 | else Set_inexactflag(); | |
329 | } | |
330 | } | |
331 | return(NOEXCEPTION); | |
332 | } |