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