]>
Commit | Line | Data |
---|---|---|
909e2256 | 1 | /* Decimal floating point support. |
95aac11b BE |
2 | Copyright (C) 2005, 2006, 2007, 2008, 2009 Free Software |
3 | Foundation, Inc. | |
909e2256 JG |
4 | |
5 | This file is part of GCC. | |
6 | ||
7 | GCC is free software; you can redistribute it and/or modify it under | |
8 | the terms of the GNU General Public License as published by the Free | |
9dcd6f09 | 9 | Software Foundation; either version 3, or (at your option) any later |
909e2256 JG |
10 | version. |
11 | ||
12 | GCC is distributed in the hope that it will be useful, but WITHOUT ANY | |
13 | WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
14 | FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |
15 | for more details. | |
16 | ||
17 | You should have received a copy of the GNU General Public License | |
9dcd6f09 NC |
18 | along with GCC; see the file COPYING3. If not see |
19 | <http://www.gnu.org/licenses/>. */ | |
909e2256 JG |
20 | |
21 | #include "config.h" | |
22 | #include "system.h" | |
23 | #include "coretypes.h" | |
24 | #include "tm.h" | |
25 | #include "tree.h" | |
26 | #include "toplev.h" | |
27 | #include "real.h" | |
28 | #include "tm_p.h" | |
29 | #include "dfp.h" | |
30 | ||
31 | /* The order of the following headers is important for making sure | |
32 | decNumber structure is large enough to hold decimal128 digits. */ | |
33 | ||
34 | #include "decimal128.h" | |
2533577f | 35 | #include "decimal128Local.h" |
909e2256 JG |
36 | #include "decimal64.h" |
37 | #include "decimal32.h" | |
38 | #include "decNumber.h" | |
39 | ||
bc6d4c3f JM |
40 | #ifndef WORDS_BIGENDIAN |
41 | #define WORDS_BIGENDIAN 0 | |
42 | #endif | |
43 | ||
909e2256 JG |
44 | /* Initialize R (a real with the decimal flag set) from DN. Can |
45 | utilize status passed in via CONTEXT, if a previous operation had | |
46 | interesting status. */ | |
47 | ||
48 | static void | |
49 | decimal_from_decnumber (REAL_VALUE_TYPE *r, decNumber *dn, decContext *context) | |
50 | { | |
51 | memset (r, 0, sizeof (REAL_VALUE_TYPE)); | |
52 | ||
53 | r->cl = rvc_normal; | |
909e2256 JG |
54 | if (decNumberIsNaN (dn)) |
55 | r->cl = rvc_nan; | |
56 | if (decNumberIsInfinite (dn)) | |
57 | r->cl = rvc_inf; | |
58 | if (context->status & DEC_Overflow) | |
59 | r->cl = rvc_inf; | |
60 | if (decNumberIsNegative (dn)) | |
61 | r->sign = 1; | |
62 | r->decimal = 1; | |
63 | ||
64 | if (r->cl != rvc_normal) | |
65 | return; | |
66 | ||
67 | decContextDefault (context, DEC_INIT_DECIMAL128); | |
68 | context->traps = 0; | |
69 | ||
70 | decimal128FromNumber ((decimal128 *) r->sig, dn, context); | |
71 | } | |
72 | ||
73 | /* Create decimal encoded R from string S. */ | |
74 | ||
75 | void | |
76 | decimal_real_from_string (REAL_VALUE_TYPE *r, const char *s) | |
77 | { | |
78 | decNumber dn; | |
79 | decContext set; | |
80 | decContextDefault (&set, DEC_INIT_DECIMAL128); | |
81 | set.traps = 0; | |
82 | ||
5f754896 | 83 | decNumberFromString (&dn, s, &set); |
909e2256 JG |
84 | |
85 | /* It would be more efficient to store directly in decNumber format, | |
86 | but that is impractical from current data structure size. | |
87 | Encoding as a decimal128 is much more compact. */ | |
88 | decimal_from_decnumber (r, &dn, &set); | |
89 | } | |
90 | ||
91 | /* Initialize a decNumber from a REAL_VALUE_TYPE. */ | |
92 | ||
93 | static void | |
94 | decimal_to_decnumber (const REAL_VALUE_TYPE *r, decNumber *dn) | |
95 | { | |
96 | decContext set; | |
97 | decContextDefault (&set, DEC_INIT_DECIMAL128); | |
98 | set.traps = 0; | |
99 | ||
100 | switch (r->cl) | |
101 | { | |
102 | case rvc_zero: | |
103 | decNumberZero (dn); | |
104 | break; | |
105 | case rvc_inf: | |
5f754896 | 106 | decNumberFromString (dn, "Infinity", &set); |
909e2256 JG |
107 | break; |
108 | case rvc_nan: | |
109 | if (r->signalling) | |
5f754896 | 110 | decNumberFromString (dn, "snan", &set); |
909e2256 | 111 | else |
5f754896 | 112 | decNumberFromString (dn, "nan", &set); |
909e2256 JG |
113 | break; |
114 | case rvc_normal: | |
115 | gcc_assert (r->decimal); | |
5f754896 | 116 | decimal128ToNumber ((const decimal128 *) r->sig, dn); |
909e2256 JG |
117 | break; |
118 | default: | |
119 | gcc_unreachable (); | |
120 | } | |
121 | ||
122 | /* Fix up sign bit. */ | |
123 | if (r->sign != decNumberIsNegative (dn)) | |
f64ad1d3 | 124 | dn->bits ^= DECNEG; |
909e2256 JG |
125 | } |
126 | ||
7292b8e4 | 127 | /* Encode a real into an IEEE 754 decimal32 type. */ |
909e2256 | 128 | |
83f676b3 | 129 | void |
909e2256 JG |
130 | encode_decimal32 (const struct real_format *fmt ATTRIBUTE_UNUSED, |
131 | long *buf, const REAL_VALUE_TYPE *r) | |
132 | { | |
133 | decNumber dn; | |
134 | decimal32 d32; | |
135 | decContext set; | |
5a5c6435 | 136 | int32_t image; |
909e2256 JG |
137 | |
138 | decContextDefault (&set, DEC_INIT_DECIMAL128); | |
139 | set.traps = 0; | |
140 | ||
b8698a0f | 141 | decimal_to_decnumber (r, &dn); |
909e2256 JG |
142 | decimal32FromNumber (&d32, &dn, &set); |
143 | ||
5a5c6435 JJ |
144 | memcpy (&image, d32.bytes, sizeof (int32_t)); |
145 | buf[0] = image; | |
909e2256 JG |
146 | } |
147 | ||
7292b8e4 | 148 | /* Decode an IEEE 754 decimal32 type into a real. */ |
909e2256 | 149 | |
83f676b3 RS |
150 | void |
151 | decode_decimal32 (const struct real_format *fmt ATTRIBUTE_UNUSED, | |
152 | REAL_VALUE_TYPE *r, const long *buf) | |
909e2256 JG |
153 | { |
154 | decNumber dn; | |
155 | decimal32 d32; | |
156 | decContext set; | |
5a5c6435 | 157 | int32_t image; |
909e2256 JG |
158 | |
159 | decContextDefault (&set, DEC_INIT_DECIMAL128); | |
160 | set.traps = 0; | |
161 | ||
5a5c6435 JJ |
162 | image = buf[0]; |
163 | memcpy (&d32.bytes, &image, sizeof (int32_t)); | |
909e2256 JG |
164 | |
165 | decimal32ToNumber (&d32, &dn); | |
b8698a0f | 166 | decimal_from_decnumber (r, &dn, &set); |
909e2256 JG |
167 | } |
168 | ||
7292b8e4 | 169 | /* Encode a real into an IEEE 754 decimal64 type. */ |
909e2256 | 170 | |
83f676b3 | 171 | void |
909e2256 JG |
172 | encode_decimal64 (const struct real_format *fmt ATTRIBUTE_UNUSED, |
173 | long *buf, const REAL_VALUE_TYPE *r) | |
174 | { | |
175 | decNumber dn; | |
176 | decimal64 d64; | |
177 | decContext set; | |
5a5c6435 | 178 | int32_t image; |
909e2256 JG |
179 | |
180 | decContextDefault (&set, DEC_INIT_DECIMAL128); | |
181 | set.traps = 0; | |
182 | ||
183 | decimal_to_decnumber (r, &dn); | |
184 | decimal64FromNumber (&d64, &dn, &set); | |
185 | ||
bc6d4c3f JM |
186 | if (WORDS_BIGENDIAN == FLOAT_WORDS_BIG_ENDIAN) |
187 | { | |
5a5c6435 JJ |
188 | memcpy (&image, &d64.bytes[0], sizeof (int32_t)); |
189 | buf[0] = image; | |
190 | memcpy (&image, &d64.bytes[4], sizeof (int32_t)); | |
191 | buf[1] = image; | |
bc6d4c3f JM |
192 | } |
193 | else | |
194 | { | |
5a5c6435 JJ |
195 | memcpy (&image, &d64.bytes[4], sizeof (int32_t)); |
196 | buf[0] = image; | |
197 | memcpy (&image, &d64.bytes[0], sizeof (int32_t)); | |
198 | buf[1] = image; | |
bc6d4c3f | 199 | } |
909e2256 JG |
200 | } |
201 | ||
7292b8e4 | 202 | /* Decode an IEEE 754 decimal64 type into a real. */ |
909e2256 | 203 | |
83f676b3 | 204 | void |
909e2256 JG |
205 | decode_decimal64 (const struct real_format *fmt ATTRIBUTE_UNUSED, |
206 | REAL_VALUE_TYPE *r, const long *buf) | |
b8698a0f | 207 | { |
909e2256 JG |
208 | decNumber dn; |
209 | decimal64 d64; | |
210 | decContext set; | |
5a5c6435 | 211 | int32_t image; |
909e2256 JG |
212 | |
213 | decContextDefault (&set, DEC_INIT_DECIMAL128); | |
214 | set.traps = 0; | |
215 | ||
bc6d4c3f JM |
216 | if (WORDS_BIGENDIAN == FLOAT_WORDS_BIG_ENDIAN) |
217 | { | |
5a5c6435 JJ |
218 | image = buf[0]; |
219 | memcpy (&d64.bytes[0], &image, sizeof (int32_t)); | |
220 | image = buf[1]; | |
221 | memcpy (&d64.bytes[4], &image, sizeof (int32_t)); | |
bc6d4c3f JM |
222 | } |
223 | else | |
224 | { | |
5a5c6435 JJ |
225 | image = buf[1]; |
226 | memcpy (&d64.bytes[0], &image, sizeof (int32_t)); | |
227 | image = buf[0]; | |
228 | memcpy (&d64.bytes[4], &image, sizeof (int32_t)); | |
bc6d4c3f | 229 | } |
909e2256 JG |
230 | |
231 | decimal64ToNumber (&d64, &dn); | |
b8698a0f | 232 | decimal_from_decnumber (r, &dn, &set); |
909e2256 JG |
233 | } |
234 | ||
7292b8e4 | 235 | /* Encode a real into an IEEE 754 decimal128 type. */ |
909e2256 | 236 | |
83f676b3 | 237 | void |
909e2256 JG |
238 | encode_decimal128 (const struct real_format *fmt ATTRIBUTE_UNUSED, |
239 | long *buf, const REAL_VALUE_TYPE *r) | |
240 | { | |
241 | decNumber dn; | |
242 | decContext set; | |
243 | decimal128 d128; | |
5a5c6435 | 244 | int32_t image; |
909e2256 JG |
245 | |
246 | decContextDefault (&set, DEC_INIT_DECIMAL128); | |
247 | set.traps = 0; | |
248 | ||
249 | decimal_to_decnumber (r, &dn); | |
250 | decimal128FromNumber (&d128, &dn, &set); | |
251 | ||
bc6d4c3f JM |
252 | if (WORDS_BIGENDIAN == FLOAT_WORDS_BIG_ENDIAN) |
253 | { | |
5a5c6435 JJ |
254 | memcpy (&image, &d128.bytes[0], sizeof (int32_t)); |
255 | buf[0] = image; | |
256 | memcpy (&image, &d128.bytes[4], sizeof (int32_t)); | |
257 | buf[1] = image; | |
258 | memcpy (&image, &d128.bytes[8], sizeof (int32_t)); | |
259 | buf[2] = image; | |
260 | memcpy (&image, &d128.bytes[12], sizeof (int32_t)); | |
261 | buf[3] = image; | |
bc6d4c3f JM |
262 | } |
263 | else | |
264 | { | |
5a5c6435 JJ |
265 | memcpy (&image, &d128.bytes[12], sizeof (int32_t)); |
266 | buf[0] = image; | |
267 | memcpy (&image, &d128.bytes[8], sizeof (int32_t)); | |
268 | buf[1] = image; | |
269 | memcpy (&image, &d128.bytes[4], sizeof (int32_t)); | |
270 | buf[2] = image; | |
271 | memcpy (&image, &d128.bytes[0], sizeof (int32_t)); | |
272 | buf[3] = image; | |
bc6d4c3f | 273 | } |
909e2256 JG |
274 | } |
275 | ||
7292b8e4 | 276 | /* Decode an IEEE 754 decimal128 type into a real. */ |
909e2256 | 277 | |
83f676b3 | 278 | void |
909e2256 JG |
279 | decode_decimal128 (const struct real_format *fmt ATTRIBUTE_UNUSED, |
280 | REAL_VALUE_TYPE *r, const long *buf) | |
281 | { | |
282 | decNumber dn; | |
283 | decimal128 d128; | |
284 | decContext set; | |
5a5c6435 | 285 | int32_t image; |
909e2256 JG |
286 | |
287 | decContextDefault (&set, DEC_INIT_DECIMAL128); | |
288 | set.traps = 0; | |
289 | ||
bc6d4c3f JM |
290 | if (WORDS_BIGENDIAN == FLOAT_WORDS_BIG_ENDIAN) |
291 | { | |
5a5c6435 JJ |
292 | image = buf[0]; |
293 | memcpy (&d128.bytes[0], &image, sizeof (int32_t)); | |
294 | image = buf[1]; | |
295 | memcpy (&d128.bytes[4], &image, sizeof (int32_t)); | |
296 | image = buf[2]; | |
297 | memcpy (&d128.bytes[8], &image, sizeof (int32_t)); | |
298 | image = buf[3]; | |
299 | memcpy (&d128.bytes[12], &image, sizeof (int32_t)); | |
bc6d4c3f JM |
300 | } |
301 | else | |
302 | { | |
5a5c6435 JJ |
303 | image = buf[3]; |
304 | memcpy (&d128.bytes[0], &image, sizeof (int32_t)); | |
305 | image = buf[2]; | |
306 | memcpy (&d128.bytes[4], &image, sizeof (int32_t)); | |
307 | image = buf[1]; | |
308 | memcpy (&d128.bytes[8], &image, sizeof (int32_t)); | |
309 | image = buf[0]; | |
310 | memcpy (&d128.bytes[12], &image, sizeof (int32_t)); | |
bc6d4c3f | 311 | } |
909e2256 JG |
312 | |
313 | decimal128ToNumber (&d128, &dn); | |
b8698a0f | 314 | decimal_from_decnumber (r, &dn, &set); |
909e2256 JG |
315 | } |
316 | ||
317 | /* Helper function to convert from a binary real internal | |
318 | representation. */ | |
319 | ||
320 | static void | |
321 | decimal_to_binary (REAL_VALUE_TYPE *to, const REAL_VALUE_TYPE *from, | |
322 | enum machine_mode mode) | |
323 | { | |
324 | char string[256]; | |
5f754896 | 325 | const decimal128 *const d128 = (const decimal128 *) from->sig; |
909e2256 JG |
326 | |
327 | decimal128ToString (d128, string); | |
328 | real_from_string3 (to, string, mode); | |
329 | } | |
330 | ||
331 | ||
332 | /* Helper function to convert from a binary real internal | |
333 | representation. */ | |
334 | ||
335 | static void | |
336 | decimal_from_binary (REAL_VALUE_TYPE *to, const REAL_VALUE_TYPE *from) | |
337 | { | |
338 | char string[256]; | |
339 | ||
340 | /* We convert to string, then to decNumber then to decimal128. */ | |
341 | real_to_decimal (string, from, sizeof (string), 0, 1); | |
342 | decimal_real_from_string (to, string); | |
343 | } | |
344 | ||
345 | /* Helper function to real.c:do_compare() to handle decimal internal | |
6416ae7f | 346 | representation including when one of the operands is still in the |
909e2256 JG |
347 | binary internal representation. */ |
348 | ||
349 | int | |
350 | decimal_do_compare (const REAL_VALUE_TYPE *a, const REAL_VALUE_TYPE *b, | |
351 | int nan_result) | |
352 | { | |
353 | decContext set; | |
354 | decNumber dn, dn2, dn3; | |
355 | REAL_VALUE_TYPE a1, b1; | |
356 | ||
357 | /* If either operand is non-decimal, create temporary versions. */ | |
358 | if (!a->decimal) | |
359 | { | |
360 | decimal_from_binary (&a1, a); | |
361 | a = &a1; | |
362 | } | |
363 | if (!b->decimal) | |
364 | { | |
365 | decimal_from_binary (&b1, b); | |
366 | b = &b1; | |
367 | } | |
b8698a0f | 368 | |
909e2256 JG |
369 | /* Convert into decNumber form for comparison operation. */ |
370 | decContextDefault (&set, DEC_INIT_DECIMAL128); | |
b8698a0f | 371 | set.traps = 0; |
5f754896 KG |
372 | decimal128ToNumber ((const decimal128 *) a->sig, &dn2); |
373 | decimal128ToNumber ((const decimal128 *) b->sig, &dn3); | |
909e2256 JG |
374 | |
375 | /* Finally, do the comparison. */ | |
376 | decNumberCompare (&dn, &dn2, &dn3, &set); | |
377 | ||
378 | /* Return the comparison result. */ | |
379 | if (decNumberIsNaN (&dn)) | |
380 | return nan_result; | |
381 | else if (decNumberIsZero (&dn)) | |
382 | return 0; | |
383 | else if (decNumberIsNegative (&dn)) | |
384 | return -1; | |
b8698a0f | 385 | else |
909e2256 JG |
386 | return 1; |
387 | } | |
388 | ||
389 | /* Helper to round_for_format, handling decimal float types. */ | |
390 | ||
391 | void | |
392 | decimal_round_for_format (const struct real_format *fmt, REAL_VALUE_TYPE *r) | |
393 | { | |
394 | decNumber dn; | |
395 | decContext set; | |
396 | ||
397 | /* Real encoding occurs later. */ | |
398 | if (r->cl != rvc_normal) | |
399 | return; | |
400 | ||
401 | decContextDefault (&set, DEC_INIT_DECIMAL128); | |
402 | set.traps = 0; | |
403 | decimal128ToNumber ((decimal128 *) r->sig, &dn); | |
404 | ||
405 | if (fmt == &decimal_quad_format) | |
406 | { | |
407 | /* The internal format is already in this format. */ | |
408 | return; | |
409 | } | |
410 | else if (fmt == &decimal_single_format) | |
411 | { | |
412 | decimal32 d32; | |
413 | decContextDefault (&set, DEC_INIT_DECIMAL32); | |
414 | set.traps = 0; | |
415 | ||
416 | decimal32FromNumber (&d32, &dn, &set); | |
417 | decimal32ToNumber (&d32, &dn); | |
418 | } | |
419 | else if (fmt == &decimal_double_format) | |
420 | { | |
421 | decimal64 d64; | |
422 | decContextDefault (&set, DEC_INIT_DECIMAL64); | |
423 | set.traps = 0; | |
424 | ||
425 | decimal64FromNumber (&d64, &dn, &set); | |
426 | decimal64ToNumber (&d64, &dn); | |
427 | } | |
428 | else | |
429 | gcc_unreachable (); | |
430 | ||
431 | decimal_from_decnumber (r, &dn, &set); | |
432 | } | |
433 | ||
434 | /* Extend or truncate to a new mode. Handles conversions between | |
435 | binary and decimal types. */ | |
436 | ||
437 | void | |
b8698a0f | 438 | decimal_real_convert (REAL_VALUE_TYPE *r, enum machine_mode mode, |
909e2256 JG |
439 | const REAL_VALUE_TYPE *a) |
440 | { | |
441 | const struct real_format *fmt = REAL_MODE_FORMAT (mode); | |
442 | ||
443 | if (a->decimal && fmt->b == 10) | |
444 | return; | |
445 | if (a->decimal) | |
446 | decimal_to_binary (r, a, mode); | |
447 | else | |
448 | decimal_from_binary (r, a); | |
449 | } | |
450 | ||
451 | /* Render R_ORIG as a decimal floating point constant. Emit DIGITS | |
452 | significant digits in the result, bounded by BUF_SIZE. If DIGITS | |
453 | is 0, choose the maximum for the representation. If | |
454 | CROP_TRAILING_ZEROS, strip trailing zeros. Currently, not honoring | |
455 | DIGITS or CROP_TRAILING_ZEROS. */ | |
456 | ||
83f676b3 RS |
457 | void |
458 | decimal_real_to_decimal (char *str, const REAL_VALUE_TYPE *r_orig, | |
459 | size_t buf_size, | |
460 | size_t digits ATTRIBUTE_UNUSED, | |
461 | int crop_trailing_zeros ATTRIBUTE_UNUSED) | |
909e2256 | 462 | { |
5f754896 | 463 | const decimal128 *const d128 = (const decimal128*) r_orig->sig; |
909e2256 JG |
464 | |
465 | /* decimal128ToString requires space for at least 24 characters; | |
466 | Require two more for suffix. */ | |
467 | gcc_assert (buf_size >= 24); | |
468 | decimal128ToString (d128, str); | |
469 | } | |
470 | ||
471 | static bool | |
472 | decimal_do_add (REAL_VALUE_TYPE *r, const REAL_VALUE_TYPE *op0, | |
473 | const REAL_VALUE_TYPE *op1, int subtract_p) | |
474 | { | |
475 | decNumber dn; | |
476 | decContext set; | |
477 | decNumber dn2, dn3; | |
478 | ||
479 | decimal_to_decnumber (op0, &dn2); | |
480 | decimal_to_decnumber (op1, &dn3); | |
481 | ||
482 | decContextDefault (&set, DEC_INIT_DECIMAL128); | |
483 | set.traps = 0; | |
484 | ||
485 | if (subtract_p) | |
486 | decNumberSubtract (&dn, &dn2, &dn3, &set); | |
b8698a0f | 487 | else |
909e2256 JG |
488 | decNumberAdd (&dn, &dn2, &dn3, &set); |
489 | ||
490 | decimal_from_decnumber (r, &dn, &set); | |
491 | ||
492 | /* Return true, if inexact. */ | |
493 | return (set.status & DEC_Inexact); | |
494 | } | |
495 | ||
496 | /* Compute R = OP0 * OP1. */ | |
497 | ||
498 | static bool | |
499 | decimal_do_multiply (REAL_VALUE_TYPE *r, const REAL_VALUE_TYPE *op0, | |
500 | const REAL_VALUE_TYPE *op1) | |
501 | { | |
502 | decContext set; | |
503 | decNumber dn, dn2, dn3; | |
504 | ||
505 | decimal_to_decnumber (op0, &dn2); | |
506 | decimal_to_decnumber (op1, &dn3); | |
507 | ||
508 | decContextDefault (&set, DEC_INIT_DECIMAL128); | |
509 | set.traps = 0; | |
510 | ||
511 | decNumberMultiply (&dn, &dn2, &dn3, &set); | |
512 | decimal_from_decnumber (r, &dn, &set); | |
513 | ||
514 | /* Return true, if inexact. */ | |
515 | return (set.status & DEC_Inexact); | |
516 | } | |
517 | ||
518 | /* Compute R = OP0 / OP1. */ | |
519 | ||
520 | static bool | |
521 | decimal_do_divide (REAL_VALUE_TYPE *r, const REAL_VALUE_TYPE *op0, | |
522 | const REAL_VALUE_TYPE *op1) | |
523 | { | |
524 | decContext set; | |
525 | decNumber dn, dn2, dn3; | |
526 | ||
527 | decimal_to_decnumber (op0, &dn2); | |
528 | decimal_to_decnumber (op1, &dn3); | |
529 | ||
530 | decContextDefault (&set, DEC_INIT_DECIMAL128); | |
531 | set.traps = 0; | |
532 | ||
533 | decNumberDivide (&dn, &dn2, &dn3, &set); | |
534 | decimal_from_decnumber (r, &dn, &set); | |
535 | ||
536 | /* Return true, if inexact. */ | |
537 | return (set.status & DEC_Inexact); | |
538 | } | |
539 | ||
540 | /* Set R to A truncated to an integral value toward zero (decimal | |
541 | floating point). */ | |
542 | ||
543 | void | |
544 | decimal_do_fix_trunc (REAL_VALUE_TYPE *r, const REAL_VALUE_TYPE *a) | |
545 | { | |
546 | decNumber dn, dn2; | |
547 | decContext set; | |
548 | ||
549 | decContextDefault (&set, DEC_INIT_DECIMAL128); | |
550 | set.traps = 0; | |
551 | set.round = DEC_ROUND_DOWN; | |
5f754896 | 552 | decimal128ToNumber ((const decimal128 *) a->sig, &dn2); |
909e2256 JG |
553 | |
554 | decNumberToIntegralValue (&dn, &dn2, &set); | |
555 | decimal_from_decnumber (r, &dn, &set); | |
556 | } | |
557 | ||
558 | /* Render decimal float value R as an integer. */ | |
559 | ||
560 | HOST_WIDE_INT | |
561 | decimal_real_to_integer (const REAL_VALUE_TYPE *r) | |
562 | { | |
563 | decContext set; | |
564 | decNumber dn, dn2, dn3; | |
565 | REAL_VALUE_TYPE to; | |
566 | char string[256]; | |
567 | ||
568 | decContextDefault (&set, DEC_INIT_DECIMAL128); | |
569 | set.traps = 0; | |
570 | set.round = DEC_ROUND_DOWN; | |
5f754896 | 571 | decimal128ToNumber ((const decimal128 *) r->sig, &dn); |
909e2256 JG |
572 | |
573 | decNumberToIntegralValue (&dn2, &dn, &set); | |
574 | decNumberZero (&dn3); | |
575 | decNumberRescale (&dn, &dn2, &dn3, &set); | |
576 | ||
577 | /* Convert to REAL_VALUE_TYPE and call appropriate conversion | |
578 | function. */ | |
579 | decNumberToString (&dn, string); | |
580 | real_from_string (&to, string); | |
581 | return real_to_integer (&to); | |
582 | } | |
583 | ||
584 | /* Likewise, but to an integer pair, HI+LOW. */ | |
585 | ||
586 | void | |
587 | decimal_real_to_integer2 (HOST_WIDE_INT *plow, HOST_WIDE_INT *phigh, | |
588 | const REAL_VALUE_TYPE *r) | |
589 | { | |
590 | decContext set; | |
591 | decNumber dn, dn2, dn3; | |
592 | REAL_VALUE_TYPE to; | |
593 | char string[256]; | |
594 | ||
595 | decContextDefault (&set, DEC_INIT_DECIMAL128); | |
596 | set.traps = 0; | |
597 | set.round = DEC_ROUND_DOWN; | |
5f754896 | 598 | decimal128ToNumber ((const decimal128 *) r->sig, &dn); |
909e2256 JG |
599 | |
600 | decNumberToIntegralValue (&dn2, &dn, &set); | |
601 | decNumberZero (&dn3); | |
602 | decNumberRescale (&dn, &dn2, &dn3, &set); | |
603 | ||
fa10beec | 604 | /* Convert to REAL_VALUE_TYPE and call appropriate conversion |
909e2256 JG |
605 | function. */ |
606 | decNumberToString (&dn, string); | |
607 | real_from_string (&to, string); | |
608 | real_to_integer2 (plow, phigh, &to); | |
609 | } | |
610 | ||
0b59f49d BE |
611 | /* Perform the decimal floating point operation described by CODE. |
612 | For a unary operation, OP1 will be NULL. This function returns | |
613 | true if the result may be inexact due to loss of precision. */ | |
909e2256 JG |
614 | |
615 | bool | |
0b59f49d | 616 | decimal_real_arithmetic (REAL_VALUE_TYPE *r, enum tree_code code, |
909e2256 JG |
617 | const REAL_VALUE_TYPE *op0, |
618 | const REAL_VALUE_TYPE *op1) | |
619 | { | |
0b59f49d | 620 | REAL_VALUE_TYPE a, b; |
909e2256 | 621 | |
0b59f49d | 622 | /* If either operand is non-decimal, create temporaries. */ |
909e2256 JG |
623 | if (!op0->decimal) |
624 | { | |
0b59f49d BE |
625 | decimal_from_binary (&a, op0); |
626 | op0 = &a; | |
909e2256 JG |
627 | } |
628 | if (op1 && !op1->decimal) | |
629 | { | |
0b59f49d BE |
630 | decimal_from_binary (&b, op1); |
631 | op1 = &b; | |
909e2256 JG |
632 | } |
633 | ||
634 | switch (code) | |
635 | { | |
636 | case PLUS_EXPR: | |
0b59f49d | 637 | return decimal_do_add (r, op0, op1, 0); |
909e2256 JG |
638 | |
639 | case MINUS_EXPR: | |
0b59f49d | 640 | return decimal_do_add (r, op0, op1, 1); |
909e2256 JG |
641 | |
642 | case MULT_EXPR: | |
0b59f49d | 643 | return decimal_do_multiply (r, op0, op1); |
909e2256 JG |
644 | |
645 | case RDIV_EXPR: | |
0b59f49d | 646 | return decimal_do_divide (r, op0, op1); |
909e2256 JG |
647 | |
648 | case MIN_EXPR: | |
649 | if (op1->cl == rvc_nan) | |
650 | *r = *op1; | |
651 | else if (real_compare (UNLT_EXPR, op0, op1)) | |
652 | *r = *op0; | |
653 | else | |
654 | *r = *op1; | |
0b59f49d | 655 | return false; |
909e2256 JG |
656 | |
657 | case MAX_EXPR: | |
658 | if (op1->cl == rvc_nan) | |
659 | *r = *op1; | |
660 | else if (real_compare (LT_EXPR, op0, op1)) | |
661 | *r = *op1; | |
662 | else | |
663 | *r = *op0; | |
0b59f49d | 664 | return false; |
909e2256 JG |
665 | |
666 | case NEGATE_EXPR: | |
667 | { | |
909e2256 | 668 | *r = *op0; |
79b87c74 MM |
669 | /* Flip sign bit. */ |
670 | decimal128FlipSign ((decimal128 *) r->sig); | |
909e2256 JG |
671 | /* Keep sign field in sync. */ |
672 | r->sign ^= 1; | |
673 | } | |
0b59f49d | 674 | return false; |
909e2256 JG |
675 | |
676 | case ABS_EXPR: | |
677 | { | |
909e2256 | 678 | *r = *op0; |
79b87c74 MM |
679 | /* Clear sign bit. */ |
680 | decimal128ClearSign ((decimal128 *) r->sig); | |
909e2256 JG |
681 | /* Keep sign field in sync. */ |
682 | r->sign = 0; | |
683 | } | |
0b59f49d | 684 | return false; |
909e2256 JG |
685 | |
686 | case FIX_TRUNC_EXPR: | |
687 | decimal_do_fix_trunc (r, op0); | |
0b59f49d | 688 | return false; |
909e2256 JG |
689 | |
690 | default: | |
691 | gcc_unreachable (); | |
692 | } | |
909e2256 JG |
693 | } |
694 | ||
695 | /* Fills R with the largest finite value representable in mode MODE. | |
696 | If SIGN is nonzero, R is set to the most negative finite value. */ | |
697 | ||
698 | void | |
699 | decimal_real_maxval (REAL_VALUE_TYPE *r, int sign, enum machine_mode mode) | |
b8698a0f | 700 | { |
5f754896 | 701 | const char *max; |
909e2256 JG |
702 | |
703 | switch (mode) | |
704 | { | |
705 | case SDmode: | |
5f754896 | 706 | max = "9.999999E96"; |
909e2256 JG |
707 | break; |
708 | case DDmode: | |
5f754896 | 709 | max = "9.999999999999999E384"; |
909e2256 JG |
710 | break; |
711 | case TDmode: | |
5f754896 | 712 | max = "9.999999999999999999999999999999999E6144"; |
909e2256 JG |
713 | break; |
714 | default: | |
715 | gcc_unreachable (); | |
716 | } | |
717 | ||
718 | decimal_real_from_string (r, max); | |
719 | if (sign) | |
79b87c74 | 720 | decimal128SetSign ((decimal128 *) r->sig, 1); |
909e2256 | 721 | } |