sfadd.c revision 1.1.4.2 1 1.1.4.2 gehenna /* $NetBSD: sfadd.c,v 1.1.4.2 2002/07/14 17:47:48 gehenna Exp $ */
2 1.1.4.2 gehenna
3 1.1.4.2 gehenna /* $OpenBSD: sfadd.c,v 1.4 2001/03/29 03:58:19 mickey Exp $ */
4 1.1.4.2 gehenna
5 1.1.4.2 gehenna /*
6 1.1.4.2 gehenna * Copyright 1996 1995 by Open Software Foundation, Inc.
7 1.1.4.2 gehenna * All Rights Reserved
8 1.1.4.2 gehenna *
9 1.1.4.2 gehenna * Permission to use, copy, modify, and distribute this software and
10 1.1.4.2 gehenna * its documentation for any purpose and without fee is hereby granted,
11 1.1.4.2 gehenna * provided that the above copyright notice appears in all copies and
12 1.1.4.2 gehenna * that both the copyright notice and this permission notice appear in
13 1.1.4.2 gehenna * supporting documentation.
14 1.1.4.2 gehenna *
15 1.1.4.2 gehenna * OSF DISCLAIMS ALL WARRANTIES WITH REGARD TO THIS SOFTWARE
16 1.1.4.2 gehenna * INCLUDING ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
17 1.1.4.2 gehenna * FOR A PARTICULAR PURPOSE.
18 1.1.4.2 gehenna *
19 1.1.4.2 gehenna * IN NO EVENT SHALL OSF BE LIABLE FOR ANY SPECIAL, INDIRECT, OR
20 1.1.4.2 gehenna * CONSEQUENTIAL DAMAGES OR ANY DAMAGES WHATSOEVER RESULTING FROM
21 1.1.4.2 gehenna * LOSS OF USE, DATA OR PROFITS, WHETHER IN ACTION OF CONTRACT,
22 1.1.4.2 gehenna * NEGLIGENCE, OR OTHER TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION
23 1.1.4.2 gehenna * WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
24 1.1.4.2 gehenna *
25 1.1.4.2 gehenna */
26 1.1.4.2 gehenna /*
27 1.1.4.2 gehenna * pmk1.1
28 1.1.4.2 gehenna */
29 1.1.4.2 gehenna /*
30 1.1.4.2 gehenna * (c) Copyright 1986 HEWLETT-PACKARD COMPANY
31 1.1.4.2 gehenna *
32 1.1.4.2 gehenna * To anyone who acknowledges that this file is provided "AS IS"
33 1.1.4.2 gehenna * without any express or implied warranty:
34 1.1.4.2 gehenna * permission to use, copy, modify, and distribute this file
35 1.1.4.2 gehenna * for any purpose is hereby granted without fee, provided that
36 1.1.4.2 gehenna * the above copyright notice and this notice appears in all
37 1.1.4.2 gehenna * copies, and that the name of Hewlett-Packard Company not be
38 1.1.4.2 gehenna * used in advertising or publicity pertaining to distribution
39 1.1.4.2 gehenna * of the software without specific, written prior permission.
40 1.1.4.2 gehenna * Hewlett-Packard Company makes no representations about the
41 1.1.4.2 gehenna * suitability of this software for any purpose.
42 1.1.4.2 gehenna */
43 1.1.4.2 gehenna
44 1.1.4.2 gehenna #include "../spmath/float.h"
45 1.1.4.2 gehenna #include "../spmath/sgl_float.h"
46 1.1.4.2 gehenna
47 1.1.4.2 gehenna /*
48 1.1.4.2 gehenna * Single_add: add two single precision values.
49 1.1.4.2 gehenna */
50 1.1.4.2 gehenna int
51 1.1.4.2 gehenna sgl_fadd(leftptr, rightptr, dstptr, status)
52 1.1.4.2 gehenna sgl_floating_point *leftptr, *rightptr, *dstptr;
53 1.1.4.2 gehenna unsigned int *status;
54 1.1.4.2 gehenna {
55 1.1.4.2 gehenna register unsigned int left, right, result, extent;
56 1.1.4.2 gehenna register unsigned int signless_upper_left, signless_upper_right, save;
57 1.1.4.2 gehenna
58 1.1.4.2 gehenna
59 1.1.4.2 gehenna register int result_exponent, right_exponent, diff_exponent;
60 1.1.4.2 gehenna register int sign_save, jumpsize;
61 1.1.4.2 gehenna register int inexact = FALSE;
62 1.1.4.2 gehenna register int underflowtrap;
63 1.1.4.2 gehenna
64 1.1.4.2 gehenna /* Create local copies of the numbers */
65 1.1.4.2 gehenna left = *leftptr;
66 1.1.4.2 gehenna right = *rightptr;
67 1.1.4.2 gehenna
68 1.1.4.2 gehenna /* A zero "save" helps discover equal operands (for later), *
69 1.1.4.2 gehenna * and is used in swapping operands (if needed). */
70 1.1.4.2 gehenna Sgl_xortointp1(left,right,/*to*/save);
71 1.1.4.2 gehenna
72 1.1.4.2 gehenna /*
73 1.1.4.2 gehenna * check first operand for NaN's or infinity
74 1.1.4.2 gehenna */
75 1.1.4.2 gehenna if ((result_exponent = Sgl_exponent(left)) == SGL_INFINITY_EXPONENT)
76 1.1.4.2 gehenna {
77 1.1.4.2 gehenna if (Sgl_iszero_mantissa(left))
78 1.1.4.2 gehenna {
79 1.1.4.2 gehenna if (Sgl_isnotnan(right))
80 1.1.4.2 gehenna {
81 1.1.4.2 gehenna if (Sgl_isinfinity(right) && save!=0)
82 1.1.4.2 gehenna {
83 1.1.4.2 gehenna /*
84 1.1.4.2 gehenna * invalid since operands are opposite signed infinity's
85 1.1.4.2 gehenna */
86 1.1.4.2 gehenna if (Is_invalidtrap_enabled()) return(INVALIDEXCEPTION);
87 1.1.4.2 gehenna Set_invalidflag();
88 1.1.4.2 gehenna Sgl_makequietnan(result);
89 1.1.4.2 gehenna *dstptr = result;
90 1.1.4.2 gehenna return(NOEXCEPTION);
91 1.1.4.2 gehenna }
92 1.1.4.2 gehenna /*
93 1.1.4.2 gehenna * return infinity
94 1.1.4.2 gehenna */
95 1.1.4.2 gehenna *dstptr = left;
96 1.1.4.2 gehenna return(NOEXCEPTION);
97 1.1.4.2 gehenna }
98 1.1.4.2 gehenna }
99 1.1.4.2 gehenna else
100 1.1.4.2 gehenna {
101 1.1.4.2 gehenna /*
102 1.1.4.2 gehenna * is NaN; signaling or quiet?
103 1.1.4.2 gehenna */
104 1.1.4.2 gehenna if (Sgl_isone_signaling(left))
105 1.1.4.2 gehenna {
106 1.1.4.2 gehenna /* trap if INVALIDTRAP enabled */
107 1.1.4.2 gehenna if (Is_invalidtrap_enabled()) return(INVALIDEXCEPTION);
108 1.1.4.2 gehenna /* make NaN quiet */
109 1.1.4.2 gehenna Set_invalidflag();
110 1.1.4.2 gehenna Sgl_set_quiet(left);
111 1.1.4.2 gehenna }
112 1.1.4.2 gehenna /*
113 1.1.4.2 gehenna * is second operand a signaling NaN?
114 1.1.4.2 gehenna */
115 1.1.4.2 gehenna else if (Sgl_is_signalingnan(right))
116 1.1.4.2 gehenna {
117 1.1.4.2 gehenna /* trap if INVALIDTRAP enabled */
118 1.1.4.2 gehenna if (Is_invalidtrap_enabled()) return(INVALIDEXCEPTION);
119 1.1.4.2 gehenna /* make NaN quiet */
120 1.1.4.2 gehenna Set_invalidflag();
121 1.1.4.2 gehenna Sgl_set_quiet(right);
122 1.1.4.2 gehenna *dstptr = right;
123 1.1.4.2 gehenna return(NOEXCEPTION);
124 1.1.4.2 gehenna }
125 1.1.4.2 gehenna /*
126 1.1.4.2 gehenna * return quiet NaN
127 1.1.4.2 gehenna */
128 1.1.4.2 gehenna *dstptr = left;
129 1.1.4.2 gehenna return(NOEXCEPTION);
130 1.1.4.2 gehenna }
131 1.1.4.2 gehenna } /* End left NaN or Infinity processing */
132 1.1.4.2 gehenna /*
133 1.1.4.2 gehenna * check second operand for NaN's or infinity
134 1.1.4.2 gehenna */
135 1.1.4.2 gehenna if (Sgl_isinfinity_exponent(right))
136 1.1.4.2 gehenna {
137 1.1.4.2 gehenna if (Sgl_iszero_mantissa(right))
138 1.1.4.2 gehenna {
139 1.1.4.2 gehenna /* return infinity */
140 1.1.4.2 gehenna *dstptr = right;
141 1.1.4.2 gehenna return(NOEXCEPTION);
142 1.1.4.2 gehenna }
143 1.1.4.2 gehenna /*
144 1.1.4.2 gehenna * is NaN; signaling or quiet?
145 1.1.4.2 gehenna */
146 1.1.4.2 gehenna if (Sgl_isone_signaling(right))
147 1.1.4.2 gehenna {
148 1.1.4.2 gehenna /* trap if INVALIDTRAP enabled */
149 1.1.4.2 gehenna if (Is_invalidtrap_enabled()) return(INVALIDEXCEPTION);
150 1.1.4.2 gehenna /* make NaN quiet */
151 1.1.4.2 gehenna Set_invalidflag();
152 1.1.4.2 gehenna Sgl_set_quiet(right);
153 1.1.4.2 gehenna }
154 1.1.4.2 gehenna /*
155 1.1.4.2 gehenna * return quiet NaN
156 1.1.4.2 gehenna */
157 1.1.4.2 gehenna *dstptr = right;
158 1.1.4.2 gehenna return(NOEXCEPTION);
159 1.1.4.2 gehenna } /* End right NaN or Infinity processing */
160 1.1.4.2 gehenna
161 1.1.4.2 gehenna /* Invariant: Must be dealing with finite numbers */
162 1.1.4.2 gehenna
163 1.1.4.2 gehenna /* Compare operands by removing the sign */
164 1.1.4.2 gehenna Sgl_copytoint_exponentmantissa(left,signless_upper_left);
165 1.1.4.2 gehenna Sgl_copytoint_exponentmantissa(right,signless_upper_right);
166 1.1.4.2 gehenna
167 1.1.4.2 gehenna /* sign difference selects add or sub operation. */
168 1.1.4.2 gehenna if(Sgl_ismagnitudeless(signless_upper_left,signless_upper_right))
169 1.1.4.2 gehenna {
170 1.1.4.2 gehenna /* Set the left operand to the larger one by XOR swap *
171 1.1.4.2 gehenna * First finish the first word using "save" */
172 1.1.4.2 gehenna Sgl_xorfromintp1(save,right,/*to*/right);
173 1.1.4.2 gehenna Sgl_xorfromintp1(save,left,/*to*/left);
174 1.1.4.2 gehenna result_exponent = Sgl_exponent(left);
175 1.1.4.2 gehenna }
176 1.1.4.2 gehenna /* Invariant: left is not smaller than right. */
177 1.1.4.2 gehenna
178 1.1.4.2 gehenna if((right_exponent = Sgl_exponent(right)) == 0)
179 1.1.4.2 gehenna {
180 1.1.4.2 gehenna /* Denormalized operands. First look for zeroes */
181 1.1.4.2 gehenna if(Sgl_iszero_mantissa(right))
182 1.1.4.2 gehenna {
183 1.1.4.2 gehenna /* right is zero */
184 1.1.4.2 gehenna if(Sgl_iszero_exponentmantissa(left))
185 1.1.4.2 gehenna {
186 1.1.4.2 gehenna /* Both operands are zeros */
187 1.1.4.2 gehenna if(Is_rounding_mode(ROUNDMINUS))
188 1.1.4.2 gehenna {
189 1.1.4.2 gehenna Sgl_or_signs(left,/*with*/right);
190 1.1.4.2 gehenna }
191 1.1.4.2 gehenna else
192 1.1.4.2 gehenna {
193 1.1.4.2 gehenna Sgl_and_signs(left,/*with*/right);
194 1.1.4.2 gehenna }
195 1.1.4.2 gehenna }
196 1.1.4.2 gehenna else
197 1.1.4.2 gehenna {
198 1.1.4.2 gehenna /* Left is not a zero and must be the result. Trapped
199 1.1.4.2 gehenna * underflows are signaled if left is denormalized. Result
200 1.1.4.2 gehenna * is always exact. */
201 1.1.4.2 gehenna if( (result_exponent == 0) && Is_underflowtrap_enabled() )
202 1.1.4.2 gehenna {
203 1.1.4.2 gehenna /* need to normalize results mantissa */
204 1.1.4.2 gehenna sign_save = Sgl_signextendedsign(left);
205 1.1.4.2 gehenna Sgl_leftshiftby1(left);
206 1.1.4.2 gehenna Sgl_normalize(left,result_exponent);
207 1.1.4.2 gehenna Sgl_set_sign(left,/*using*/sign_save);
208 1.1.4.2 gehenna Sgl_setwrapped_exponent(left,result_exponent,unfl);
209 1.1.4.2 gehenna *dstptr = left;
210 1.1.4.2 gehenna return(UNDERFLOWEXCEPTION);
211 1.1.4.2 gehenna }
212 1.1.4.2 gehenna }
213 1.1.4.2 gehenna *dstptr = left;
214 1.1.4.2 gehenna return(NOEXCEPTION);
215 1.1.4.2 gehenna }
216 1.1.4.2 gehenna
217 1.1.4.2 gehenna /* Neither are zeroes */
218 1.1.4.2 gehenna Sgl_clear_sign(right); /* Exponent is already cleared */
219 1.1.4.2 gehenna if(result_exponent == 0 )
220 1.1.4.2 gehenna {
221 1.1.4.2 gehenna /* Both operands are denormalized. The result must be exact
222 1.1.4.2 gehenna * and is simply calculated. A sum could become normalized and a
223 1.1.4.2 gehenna * difference could cancel to a true zero. */
224 1.1.4.2 gehenna if( (/*signed*/int) save < 0 )
225 1.1.4.2 gehenna {
226 1.1.4.2 gehenna Sgl_subtract(left,/*minus*/right,/*into*/result);
227 1.1.4.2 gehenna if(Sgl_iszero_mantissa(result))
228 1.1.4.2 gehenna {
229 1.1.4.2 gehenna if(Is_rounding_mode(ROUNDMINUS))
230 1.1.4.2 gehenna {
231 1.1.4.2 gehenna Sgl_setone_sign(result);
232 1.1.4.2 gehenna }
233 1.1.4.2 gehenna else
234 1.1.4.2 gehenna {
235 1.1.4.2 gehenna Sgl_setzero_sign(result);
236 1.1.4.2 gehenna }
237 1.1.4.2 gehenna *dstptr = result;
238 1.1.4.2 gehenna return(NOEXCEPTION);
239 1.1.4.2 gehenna }
240 1.1.4.2 gehenna }
241 1.1.4.2 gehenna else
242 1.1.4.2 gehenna {
243 1.1.4.2 gehenna Sgl_addition(left,right,/*into*/result);
244 1.1.4.2 gehenna if(Sgl_isone_hidden(result))
245 1.1.4.2 gehenna {
246 1.1.4.2 gehenna *dstptr = result;
247 1.1.4.2 gehenna return(NOEXCEPTION);
248 1.1.4.2 gehenna }
249 1.1.4.2 gehenna }
250 1.1.4.2 gehenna if(Is_underflowtrap_enabled())
251 1.1.4.2 gehenna {
252 1.1.4.2 gehenna /* need to normalize result */
253 1.1.4.2 gehenna sign_save = Sgl_signextendedsign(result);
254 1.1.4.2 gehenna Sgl_leftshiftby1(result);
255 1.1.4.2 gehenna Sgl_normalize(result,result_exponent);
256 1.1.4.2 gehenna Sgl_set_sign(result,/*using*/sign_save);
257 1.1.4.2 gehenna Sgl_setwrapped_exponent(result,result_exponent,unfl);
258 1.1.4.2 gehenna *dstptr = result;
259 1.1.4.2 gehenna return(UNDERFLOWEXCEPTION);
260 1.1.4.2 gehenna }
261 1.1.4.2 gehenna *dstptr = result;
262 1.1.4.2 gehenna return(NOEXCEPTION);
263 1.1.4.2 gehenna }
264 1.1.4.2 gehenna right_exponent = 1; /* Set exponent to reflect different bias
265 1.1.4.2 gehenna * with denomalized numbers. */
266 1.1.4.2 gehenna }
267 1.1.4.2 gehenna else
268 1.1.4.2 gehenna {
269 1.1.4.2 gehenna Sgl_clear_signexponent_set_hidden(right);
270 1.1.4.2 gehenna }
271 1.1.4.2 gehenna Sgl_clear_exponent_set_hidden(left);
272 1.1.4.2 gehenna diff_exponent = result_exponent - right_exponent;
273 1.1.4.2 gehenna
274 1.1.4.2 gehenna /*
275 1.1.4.2 gehenna * Special case alignment of operands that would force alignment
276 1.1.4.2 gehenna * beyond the extent of the extension. A further optimization
277 1.1.4.2 gehenna * could special case this but only reduces the path length for this
278 1.1.4.2 gehenna * infrequent case.
279 1.1.4.2 gehenna */
280 1.1.4.2 gehenna if(diff_exponent > SGL_THRESHOLD)
281 1.1.4.2 gehenna {
282 1.1.4.2 gehenna diff_exponent = SGL_THRESHOLD;
283 1.1.4.2 gehenna }
284 1.1.4.2 gehenna
285 1.1.4.2 gehenna /* Align right operand by shifting to right */
286 1.1.4.2 gehenna Sgl_right_align(/*operand*/right,/*shifted by*/diff_exponent,
287 1.1.4.2 gehenna /*and lower to*/extent);
288 1.1.4.2 gehenna
289 1.1.4.2 gehenna /* Treat sum and difference of the operands separately. */
290 1.1.4.2 gehenna if( (/*signed*/int) save < 0 )
291 1.1.4.2 gehenna {
292 1.1.4.2 gehenna /*
293 1.1.4.2 gehenna * Difference of the two operands. Their can be no overflow. A
294 1.1.4.2 gehenna * borrow can occur out of the hidden bit and force a post
295 1.1.4.2 gehenna * normalization phase.
296 1.1.4.2 gehenna */
297 1.1.4.2 gehenna Sgl_subtract_withextension(left,/*minus*/right,/*with*/extent,/*into*/result);
298 1.1.4.2 gehenna if(Sgl_iszero_hidden(result))
299 1.1.4.2 gehenna {
300 1.1.4.2 gehenna /* Handle normalization */
301 1.1.4.2 gehenna /* A straight foward algorithm would now shift the result
302 1.1.4.2 gehenna * and extension left until the hidden bit becomes one. Not
303 1.1.4.2 gehenna * all of the extension bits need participate in the shift.
304 1.1.4.2 gehenna * Only the two most significant bits (round and guard) are
305 1.1.4.2 gehenna * needed. If only a single shift is needed then the guard
306 1.1.4.2 gehenna * bit becomes a significant low order bit and the extension
307 1.1.4.2 gehenna * must participate in the rounding. If more than a single
308 1.1.4.2 gehenna * shift is needed, then all bits to the right of the guard
309 1.1.4.2 gehenna * bit are zeros, and the guard bit may or may not be zero. */
310 1.1.4.2 gehenna sign_save = Sgl_signextendedsign(result);
311 1.1.4.2 gehenna Sgl_leftshiftby1_withextent(result,extent,result);
312 1.1.4.2 gehenna
313 1.1.4.2 gehenna /* Need to check for a zero result. The sign and exponent
314 1.1.4.2 gehenna * fields have already been zeroed. The more efficient test
315 1.1.4.2 gehenna * of the full object can be used.
316 1.1.4.2 gehenna */
317 1.1.4.2 gehenna if(Sgl_iszero(result))
318 1.1.4.2 gehenna /* Must have been "x-x" or "x+(-x)". */
319 1.1.4.2 gehenna {
320 1.1.4.2 gehenna if(Is_rounding_mode(ROUNDMINUS)) Sgl_setone_sign(result);
321 1.1.4.2 gehenna *dstptr = result;
322 1.1.4.2 gehenna return(NOEXCEPTION);
323 1.1.4.2 gehenna }
324 1.1.4.2 gehenna result_exponent--;
325 1.1.4.2 gehenna /* Look to see if normalization is finished. */
326 1.1.4.2 gehenna if(Sgl_isone_hidden(result))
327 1.1.4.2 gehenna {
328 1.1.4.2 gehenna if(result_exponent==0)
329 1.1.4.2 gehenna {
330 1.1.4.2 gehenna /* Denormalized, exponent should be zero. Left operand *
331 1.1.4.2 gehenna * was normalized, so extent (guard, round) was zero */
332 1.1.4.2 gehenna goto underflow;
333 1.1.4.2 gehenna }
334 1.1.4.2 gehenna else
335 1.1.4.2 gehenna {
336 1.1.4.2 gehenna /* No further normalization is needed. */
337 1.1.4.2 gehenna Sgl_set_sign(result,/*using*/sign_save);
338 1.1.4.2 gehenna Ext_leftshiftby1(extent);
339 1.1.4.2 gehenna goto round;
340 1.1.4.2 gehenna }
341 1.1.4.2 gehenna }
342 1.1.4.2 gehenna
343 1.1.4.2 gehenna /* Check for denormalized, exponent should be zero. Left *
344 1.1.4.2 gehenna * operand was normalized, so extent (guard, round) was zero */
345 1.1.4.2 gehenna if(!(underflowtrap = Is_underflowtrap_enabled()) &&
346 1.1.4.2 gehenna result_exponent==0) goto underflow;
347 1.1.4.2 gehenna
348 1.1.4.2 gehenna /* Shift extension to complete one bit of normalization and
349 1.1.4.2 gehenna * update exponent. */
350 1.1.4.2 gehenna Ext_leftshiftby1(extent);
351 1.1.4.2 gehenna
352 1.1.4.2 gehenna /* Discover first one bit to determine shift amount. Use a
353 1.1.4.2 gehenna * modified binary search. We have already shifted the result
354 1.1.4.2 gehenna * one position right and still not found a one so the remainder
355 1.1.4.2 gehenna * of the extension must be zero and simplifies rounding. */
356 1.1.4.2 gehenna /* Scan bytes */
357 1.1.4.2 gehenna while(Sgl_iszero_hiddenhigh7mantissa(result))
358 1.1.4.2 gehenna {
359 1.1.4.2 gehenna Sgl_leftshiftby8(result);
360 1.1.4.2 gehenna if((result_exponent -= 8) <= 0 && !underflowtrap)
361 1.1.4.2 gehenna goto underflow;
362 1.1.4.2 gehenna }
363 1.1.4.2 gehenna /* Now narrow it down to the nibble */
364 1.1.4.2 gehenna if(Sgl_iszero_hiddenhigh3mantissa(result))
365 1.1.4.2 gehenna {
366 1.1.4.2 gehenna /* The lower nibble contains the normalizing one */
367 1.1.4.2 gehenna Sgl_leftshiftby4(result);
368 1.1.4.2 gehenna if((result_exponent -= 4) <= 0 && !underflowtrap)
369 1.1.4.2 gehenna goto underflow;
370 1.1.4.2 gehenna }
371 1.1.4.2 gehenna /* Select case were first bit is set (already normalized)
372 1.1.4.2 gehenna * otherwise select the proper shift. */
373 1.1.4.2 gehenna if((jumpsize = Sgl_hiddenhigh3mantissa(result)) > 7)
374 1.1.4.2 gehenna {
375 1.1.4.2 gehenna /* Already normalized */
376 1.1.4.2 gehenna if(result_exponent <= 0) goto underflow;
377 1.1.4.2 gehenna Sgl_set_sign(result,/*using*/sign_save);
378 1.1.4.2 gehenna Sgl_set_exponent(result,/*using*/result_exponent);
379 1.1.4.2 gehenna *dstptr = result;
380 1.1.4.2 gehenna return(NOEXCEPTION);
381 1.1.4.2 gehenna }
382 1.1.4.2 gehenna Sgl_sethigh4bits(result,/*using*/sign_save);
383 1.1.4.2 gehenna switch(jumpsize)
384 1.1.4.2 gehenna {
385 1.1.4.2 gehenna case 1:
386 1.1.4.2 gehenna {
387 1.1.4.2 gehenna Sgl_leftshiftby3(result);
388 1.1.4.2 gehenna result_exponent -= 3;
389 1.1.4.2 gehenna break;
390 1.1.4.2 gehenna }
391 1.1.4.2 gehenna case 2:
392 1.1.4.2 gehenna case 3:
393 1.1.4.2 gehenna {
394 1.1.4.2 gehenna Sgl_leftshiftby2(result);
395 1.1.4.2 gehenna result_exponent -= 2;
396 1.1.4.2 gehenna break;
397 1.1.4.2 gehenna }
398 1.1.4.2 gehenna case 4:
399 1.1.4.2 gehenna case 5:
400 1.1.4.2 gehenna case 6:
401 1.1.4.2 gehenna case 7:
402 1.1.4.2 gehenna {
403 1.1.4.2 gehenna Sgl_leftshiftby1(result);
404 1.1.4.2 gehenna result_exponent -= 1;
405 1.1.4.2 gehenna break;
406 1.1.4.2 gehenna }
407 1.1.4.2 gehenna }
408 1.1.4.2 gehenna if(result_exponent > 0)
409 1.1.4.2 gehenna {
410 1.1.4.2 gehenna Sgl_set_exponent(result,/*using*/result_exponent);
411 1.1.4.2 gehenna *dstptr = result;
412 1.1.4.2 gehenna return(NOEXCEPTION); /* Sign bit is already set */
413 1.1.4.2 gehenna }
414 1.1.4.2 gehenna /* Fixup potential underflows */
415 1.1.4.2 gehenna underflow:
416 1.1.4.2 gehenna if(Is_underflowtrap_enabled())
417 1.1.4.2 gehenna {
418 1.1.4.2 gehenna Sgl_set_sign(result,sign_save);
419 1.1.4.2 gehenna Sgl_setwrapped_exponent(result,result_exponent,unfl);
420 1.1.4.2 gehenna *dstptr = result;
421 1.1.4.2 gehenna /* inexact = FALSE; */
422 1.1.4.2 gehenna return(UNDERFLOWEXCEPTION);
423 1.1.4.2 gehenna }
424 1.1.4.2 gehenna /*
425 1.1.4.2 gehenna * Since we cannot get an inexact denormalized result,
426 1.1.4.2 gehenna * we can now return.
427 1.1.4.2 gehenna */
428 1.1.4.2 gehenna Sgl_right_align(result,/*by*/(1-result_exponent),extent);
429 1.1.4.2 gehenna Sgl_clear_signexponent(result);
430 1.1.4.2 gehenna Sgl_set_sign(result,sign_save);
431 1.1.4.2 gehenna *dstptr = result;
432 1.1.4.2 gehenna return(NOEXCEPTION);
433 1.1.4.2 gehenna } /* end if(hidden...)... */
434 1.1.4.2 gehenna /* Fall through and round */
435 1.1.4.2 gehenna } /* end if(save < 0)... */
436 1.1.4.2 gehenna else
437 1.1.4.2 gehenna {
438 1.1.4.2 gehenna /* Add magnitudes */
439 1.1.4.2 gehenna Sgl_addition(left,right,/*to*/result);
440 1.1.4.2 gehenna if(Sgl_isone_hiddenoverflow(result))
441 1.1.4.2 gehenna {
442 1.1.4.2 gehenna /* Prenormalization required. */
443 1.1.4.2 gehenna Sgl_rightshiftby1_withextent(result,extent,extent);
444 1.1.4.2 gehenna Sgl_arithrightshiftby1(result);
445 1.1.4.2 gehenna result_exponent++;
446 1.1.4.2 gehenna } /* end if hiddenoverflow... */
447 1.1.4.2 gehenna } /* end else ...add magnitudes... */
448 1.1.4.2 gehenna
449 1.1.4.2 gehenna /* Round the result. If the extension is all zeros,then the result is
450 1.1.4.2 gehenna * exact. Otherwise round in the correct direction. No underflow is
451 1.1.4.2 gehenna * possible. If a postnormalization is necessary, then the mantissa is
452 1.1.4.2 gehenna * all zeros so no shift is needed. */
453 1.1.4.2 gehenna round:
454 1.1.4.2 gehenna if(Ext_isnotzero(extent))
455 1.1.4.2 gehenna {
456 1.1.4.2 gehenna inexact = TRUE;
457 1.1.4.2 gehenna switch(Rounding_mode())
458 1.1.4.2 gehenna {
459 1.1.4.2 gehenna case ROUNDNEAREST: /* The default. */
460 1.1.4.2 gehenna if(Ext_isone_sign(extent))
461 1.1.4.2 gehenna {
462 1.1.4.2 gehenna /* at least 1/2 ulp */
463 1.1.4.2 gehenna if(Ext_isnotzero_lower(extent) ||
464 1.1.4.2 gehenna Sgl_isone_lowmantissa(result))
465 1.1.4.2 gehenna {
466 1.1.4.2 gehenna /* either exactly half way and odd or more than 1/2ulp */
467 1.1.4.2 gehenna Sgl_increment(result);
468 1.1.4.2 gehenna }
469 1.1.4.2 gehenna }
470 1.1.4.2 gehenna break;
471 1.1.4.2 gehenna
472 1.1.4.2 gehenna case ROUNDPLUS:
473 1.1.4.2 gehenna if(Sgl_iszero_sign(result))
474 1.1.4.2 gehenna {
475 1.1.4.2 gehenna /* Round up positive results */
476 1.1.4.2 gehenna Sgl_increment(result);
477 1.1.4.2 gehenna }
478 1.1.4.2 gehenna break;
479 1.1.4.2 gehenna
480 1.1.4.2 gehenna case ROUNDMINUS:
481 1.1.4.2 gehenna if(Sgl_isone_sign(result))
482 1.1.4.2 gehenna {
483 1.1.4.2 gehenna /* Round down negative results */
484 1.1.4.2 gehenna Sgl_increment(result);
485 1.1.4.2 gehenna }
486 1.1.4.2 gehenna
487 1.1.4.2 gehenna case ROUNDZERO:;
488 1.1.4.2 gehenna /* truncate is simple */
489 1.1.4.2 gehenna } /* end switch... */
490 1.1.4.2 gehenna if(Sgl_isone_hiddenoverflow(result)) result_exponent++;
491 1.1.4.2 gehenna }
492 1.1.4.2 gehenna if(result_exponent == SGL_INFINITY_EXPONENT)
493 1.1.4.2 gehenna {
494 1.1.4.2 gehenna /* Overflow */
495 1.1.4.2 gehenna if(Is_overflowtrap_enabled())
496 1.1.4.2 gehenna {
497 1.1.4.2 gehenna Sgl_setwrapped_exponent(result,result_exponent,ovfl);
498 1.1.4.2 gehenna *dstptr = result;
499 1.1.4.2 gehenna if (inexact) {
500 1.1.4.2 gehenna if (Is_inexacttrap_enabled())
501 1.1.4.2 gehenna return(OVERFLOWEXCEPTION | INEXACTEXCEPTION);
502 1.1.4.2 gehenna else Set_inexactflag();
503 1.1.4.2 gehenna }
504 1.1.4.2 gehenna return(OVERFLOWEXCEPTION);
505 1.1.4.2 gehenna }
506 1.1.4.2 gehenna else
507 1.1.4.2 gehenna {
508 1.1.4.2 gehenna Set_overflowflag();
509 1.1.4.2 gehenna inexact = TRUE;
510 1.1.4.2 gehenna Sgl_setoverflow(result);
511 1.1.4.2 gehenna }
512 1.1.4.2 gehenna }
513 1.1.4.2 gehenna else Sgl_set_exponent(result,result_exponent);
514 1.1.4.2 gehenna *dstptr = result;
515 1.1.4.2 gehenna if(inexact) {
516 1.1.4.2 gehenna if(Is_inexacttrap_enabled()) return(INEXACTEXCEPTION);
517 1.1.4.2 gehenna else Set_inexactflag();
518 1.1.4.2 gehenna }
519 1.1.4.2 gehenna return(NOEXCEPTION);
520 1.1.4.2 gehenna }
521