Mercurial > hg > CbC > CbC_gcc
annotate gcc/dfp.c @ 59:5b5b9ea5b220
fix
author | Shinji KONO <kono@ie.u-ryukyu.ac.jp> |
---|---|
date | Mon, 15 Feb 2010 17:22:24 +0900 |
parents | 77e2b8dfacca |
children | b7f97abdc517 |
rev | line source |
---|---|
0 | 1 /* Decimal floating point support. |
55
77e2b8dfacca
update it from 4.4.3 to 4.5.0
ryoma <e075725@ie.u-ryukyu.ac.jp>
parents:
19
diff
changeset
|
2 Copyright (C) 2005, 2006, 2007, 2008, 2009 Free Software |
77e2b8dfacca
update it from 4.4.3 to 4.5.0
ryoma <e075725@ie.u-ryukyu.ac.jp>
parents:
19
diff
changeset
|
3 Foundation, Inc. |
0 | 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 | |
9 Software Foundation; either version 3, or (at your option) any later | |
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 | |
18 along with GCC; see the file COPYING3. If not see | |
19 <http://www.gnu.org/licenses/>. */ | |
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" | |
35 #include "decimal128Local.h" | |
36 #include "decimal64.h" | |
37 #include "decimal32.h" | |
38 #include "decNumber.h" | |
39 | |
40 #ifndef WORDS_BIGENDIAN | |
41 #define WORDS_BIGENDIAN 0 | |
42 #endif | |
43 | |
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; | |
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 | |
83 decNumberFromString (&dn, s, &set); | |
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: | |
106 decNumberFromString (dn, "Infinity", &set); | |
107 break; | |
108 case rvc_nan: | |
109 if (r->signalling) | |
110 decNumberFromString (dn, "snan", &set); | |
111 else | |
112 decNumberFromString (dn, "nan", &set); | |
113 break; | |
114 case rvc_normal: | |
115 gcc_assert (r->decimal); | |
116 decimal128ToNumber ((const decimal128 *) r->sig, dn); | |
117 break; | |
118 default: | |
119 gcc_unreachable (); | |
120 } | |
121 | |
122 /* Fix up sign bit. */ | |
123 if (r->sign != decNumberIsNegative (dn)) | |
124 dn->bits ^= DECNEG; | |
125 } | |
126 | |
127 /* Encode a real into an IEEE 754 decimal32 type. */ | |
128 | |
129 void | |
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; | |
19
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
136 int32_t image; |
0 | 137 |
138 decContextDefault (&set, DEC_INIT_DECIMAL128); | |
139 set.traps = 0; | |
140 | |
55
77e2b8dfacca
update it from 4.4.3 to 4.5.0
ryoma <e075725@ie.u-ryukyu.ac.jp>
parents:
19
diff
changeset
|
141 decimal_to_decnumber (r, &dn); |
0 | 142 decimal32FromNumber (&d32, &dn, &set); |
143 | |
19
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
144 memcpy (&image, d32.bytes, sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
145 buf[0] = image; |
0 | 146 } |
147 | |
148 /* Decode an IEEE 754 decimal32 type into a real. */ | |
149 | |
150 void | |
151 decode_decimal32 (const struct real_format *fmt ATTRIBUTE_UNUSED, | |
152 REAL_VALUE_TYPE *r, const long *buf) | |
153 { | |
154 decNumber dn; | |
155 decimal32 d32; | |
156 decContext set; | |
19
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
157 int32_t image; |
0 | 158 |
159 decContextDefault (&set, DEC_INIT_DECIMAL128); | |
160 set.traps = 0; | |
161 | |
19
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
162 image = buf[0]; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
163 memcpy (&d32.bytes, &image, sizeof (int32_t)); |
0 | 164 |
165 decimal32ToNumber (&d32, &dn); | |
55
77e2b8dfacca
update it from 4.4.3 to 4.5.0
ryoma <e075725@ie.u-ryukyu.ac.jp>
parents:
19
diff
changeset
|
166 decimal_from_decnumber (r, &dn, &set); |
0 | 167 } |
168 | |
169 /* Encode a real into an IEEE 754 decimal64 type. */ | |
170 | |
171 void | |
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; | |
19
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
178 int32_t image; |
0 | 179 |
180 decContextDefault (&set, DEC_INIT_DECIMAL128); | |
181 set.traps = 0; | |
182 | |
183 decimal_to_decnumber (r, &dn); | |
184 decimal64FromNumber (&d64, &dn, &set); | |
185 | |
186 if (WORDS_BIGENDIAN == FLOAT_WORDS_BIG_ENDIAN) | |
187 { | |
19
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
188 memcpy (&image, &d64.bytes[0], sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
189 buf[0] = image; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
190 memcpy (&image, &d64.bytes[4], sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
191 buf[1] = image; |
0 | 192 } |
193 else | |
194 { | |
19
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
195 memcpy (&image, &d64.bytes[4], sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
196 buf[0] = image; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
197 memcpy (&image, &d64.bytes[0], sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
198 buf[1] = image; |
0 | 199 } |
200 } | |
201 | |
202 /* Decode an IEEE 754 decimal64 type into a real. */ | |
203 | |
204 void | |
205 decode_decimal64 (const struct real_format *fmt ATTRIBUTE_UNUSED, | |
206 REAL_VALUE_TYPE *r, const long *buf) | |
55
77e2b8dfacca
update it from 4.4.3 to 4.5.0
ryoma <e075725@ie.u-ryukyu.ac.jp>
parents:
19
diff
changeset
|
207 { |
0 | 208 decNumber dn; |
209 decimal64 d64; | |
210 decContext set; | |
19
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
211 int32_t image; |
0 | 212 |
213 decContextDefault (&set, DEC_INIT_DECIMAL128); | |
214 set.traps = 0; | |
215 | |
216 if (WORDS_BIGENDIAN == FLOAT_WORDS_BIG_ENDIAN) | |
217 { | |
19
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
218 image = buf[0]; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
219 memcpy (&d64.bytes[0], &image, sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
220 image = buf[1]; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
221 memcpy (&d64.bytes[4], &image, sizeof (int32_t)); |
0 | 222 } |
223 else | |
224 { | |
19
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
225 image = buf[1]; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
226 memcpy (&d64.bytes[0], &image, sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
227 image = buf[0]; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
228 memcpy (&d64.bytes[4], &image, sizeof (int32_t)); |
0 | 229 } |
230 | |
231 decimal64ToNumber (&d64, &dn); | |
55
77e2b8dfacca
update it from 4.4.3 to 4.5.0
ryoma <e075725@ie.u-ryukyu.ac.jp>
parents:
19
diff
changeset
|
232 decimal_from_decnumber (r, &dn, &set); |
0 | 233 } |
234 | |
235 /* Encode a real into an IEEE 754 decimal128 type. */ | |
236 | |
237 void | |
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; | |
19
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
244 int32_t image; |
0 | 245 |
246 decContextDefault (&set, DEC_INIT_DECIMAL128); | |
247 set.traps = 0; | |
248 | |
249 decimal_to_decnumber (r, &dn); | |
250 decimal128FromNumber (&d128, &dn, &set); | |
251 | |
252 if (WORDS_BIGENDIAN == FLOAT_WORDS_BIG_ENDIAN) | |
253 { | |
19
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
254 memcpy (&image, &d128.bytes[0], sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
255 buf[0] = image; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
256 memcpy (&image, &d128.bytes[4], sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
257 buf[1] = image; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
258 memcpy (&image, &d128.bytes[8], sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
259 buf[2] = image; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
260 memcpy (&image, &d128.bytes[12], sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
261 buf[3] = image; |
0 | 262 } |
263 else | |
264 { | |
19
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
265 memcpy (&image, &d128.bytes[12], sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
266 buf[0] = image; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
267 memcpy (&image, &d128.bytes[8], sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
268 buf[1] = image; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
269 memcpy (&image, &d128.bytes[4], sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
270 buf[2] = image; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
271 memcpy (&image, &d128.bytes[0], sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
272 buf[3] = image; |
0 | 273 } |
274 } | |
275 | |
276 /* Decode an IEEE 754 decimal128 type into a real. */ | |
277 | |
278 void | |
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; | |
19
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
285 int32_t image; |
0 | 286 |
287 decContextDefault (&set, DEC_INIT_DECIMAL128); | |
288 set.traps = 0; | |
289 | |
290 if (WORDS_BIGENDIAN == FLOAT_WORDS_BIG_ENDIAN) | |
291 { | |
19
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
292 image = buf[0]; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
293 memcpy (&d128.bytes[0], &image, sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
294 image = buf[1]; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
295 memcpy (&d128.bytes[4], &image, sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
296 image = buf[2]; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
297 memcpy (&d128.bytes[8], &image, sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
298 image = buf[3]; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
299 memcpy (&d128.bytes[12], &image, sizeof (int32_t)); |
0 | 300 } |
301 else | |
302 { | |
19
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
303 image = buf[3]; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
304 memcpy (&d128.bytes[0], &image, sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
305 image = buf[2]; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
306 memcpy (&d128.bytes[4], &image, sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
307 image = buf[1]; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
308 memcpy (&d128.bytes[8], &image, sizeof (int32_t)); |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
309 image = buf[0]; |
58ad6c70ea60
update gcc from 4.4.0 to 4.4.1.
kent@firefly.cr.ie.u-ryukyu.ac.jp
parents:
0
diff
changeset
|
310 memcpy (&d128.bytes[12], &image, sizeof (int32_t)); |
0 | 311 } |
312 | |
313 decimal128ToNumber (&d128, &dn); | |
55
77e2b8dfacca
update it from 4.4.3 to 4.5.0
ryoma <e075725@ie.u-ryukyu.ac.jp>
parents:
19
diff
changeset
|
314 decimal_from_decnumber (r, &dn, &set); |
0 | 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]; | |
325 const decimal128 *const d128 = (const decimal128 *) from->sig; | |
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 | |
346 representation including when one of the operands is still in the | |
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 } | |
55
77e2b8dfacca
update it from 4.4.3 to 4.5.0
ryoma <e075725@ie.u-ryukyu.ac.jp>
parents:
19
diff
changeset
|
368 |
0 | 369 /* Convert into decNumber form for comparison operation. */ |
370 decContextDefault (&set, DEC_INIT_DECIMAL128); | |
55
77e2b8dfacca
update it from 4.4.3 to 4.5.0
ryoma <e075725@ie.u-ryukyu.ac.jp>
parents:
19
diff
changeset
|
371 set.traps = 0; |
0 | 372 decimal128ToNumber ((const decimal128 *) a->sig, &dn2); |
373 decimal128ToNumber ((const decimal128 *) b->sig, &dn3); | |
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; | |
55
77e2b8dfacca
update it from 4.4.3 to 4.5.0
ryoma <e075725@ie.u-ryukyu.ac.jp>
parents:
19
diff
changeset
|
385 else |
0 | 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 | |
55
77e2b8dfacca
update it from 4.4.3 to 4.5.0
ryoma <e075725@ie.u-ryukyu.ac.jp>
parents:
19
diff
changeset
|
438 decimal_real_convert (REAL_VALUE_TYPE *r, enum machine_mode mode, |
0 | 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 | |
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) | |
462 { | |
463 const decimal128 *const d128 = (const decimal128*) r_orig->sig; | |
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); | |
55
77e2b8dfacca
update it from 4.4.3 to 4.5.0
ryoma <e075725@ie.u-ryukyu.ac.jp>
parents:
19
diff
changeset
|
487 else |
0 | 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; | |
552 decimal128ToNumber ((const decimal128 *) a->sig, &dn2); | |
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; | |
571 decimal128ToNumber ((const decimal128 *) r->sig, &dn); | |
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; | |
598 decimal128ToNumber ((const decimal128 *) r->sig, &dn); | |
599 | |
600 decNumberToIntegralValue (&dn2, &dn, &set); | |
601 decNumberZero (&dn3); | |
602 decNumberRescale (&dn, &dn2, &dn3, &set); | |
603 | |
604 /* Convert to REAL_VALUE_TYPE and call appropriate conversion | |
605 function. */ | |
606 decNumberToString (&dn, string); | |
607 real_from_string (&to, string); | |
608 real_to_integer2 (plow, phigh, &to); | |
609 } | |
610 | |
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. */ | |
614 | |
615 bool | |
616 decimal_real_arithmetic (REAL_VALUE_TYPE *r, enum tree_code code, | |
617 const REAL_VALUE_TYPE *op0, | |
618 const REAL_VALUE_TYPE *op1) | |
619 { | |
620 REAL_VALUE_TYPE a, b; | |
621 | |
622 /* If either operand is non-decimal, create temporaries. */ | |
623 if (!op0->decimal) | |
624 { | |
625 decimal_from_binary (&a, op0); | |
626 op0 = &a; | |
627 } | |
628 if (op1 && !op1->decimal) | |
629 { | |
630 decimal_from_binary (&b, op1); | |
631 op1 = &b; | |
632 } | |
633 | |
634 switch (code) | |
635 { | |
636 case PLUS_EXPR: | |
637 return decimal_do_add (r, op0, op1, 0); | |
638 | |
639 case MINUS_EXPR: | |
640 return decimal_do_add (r, op0, op1, 1); | |
641 | |
642 case MULT_EXPR: | |
643 return decimal_do_multiply (r, op0, op1); | |
644 | |
645 case RDIV_EXPR: | |
646 return decimal_do_divide (r, op0, op1); | |
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; | |
655 return false; | |
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; | |
664 return false; | |
665 | |
666 case NEGATE_EXPR: | |
667 { | |
668 *r = *op0; | |
669 /* Flip sign bit. */ | |
670 decimal128FlipSign ((decimal128 *) r->sig); | |
671 /* Keep sign field in sync. */ | |
672 r->sign ^= 1; | |
673 } | |
674 return false; | |
675 | |
676 case ABS_EXPR: | |
677 { | |
678 *r = *op0; | |
679 /* Clear sign bit. */ | |
680 decimal128ClearSign ((decimal128 *) r->sig); | |
681 /* Keep sign field in sync. */ | |
682 r->sign = 0; | |
683 } | |
684 return false; | |
685 | |
686 case FIX_TRUNC_EXPR: | |
687 decimal_do_fix_trunc (r, op0); | |
688 return false; | |
689 | |
690 default: | |
691 gcc_unreachable (); | |
692 } | |
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) | |
55
77e2b8dfacca
update it from 4.4.3 to 4.5.0
ryoma <e075725@ie.u-ryukyu.ac.jp>
parents:
19
diff
changeset
|
700 { |
0 | 701 const char *max; |
702 | |
703 switch (mode) | |
704 { | |
705 case SDmode: | |
706 max = "9.999999E96"; | |
707 break; | |
708 case DDmode: | |
709 max = "9.999999999999999E384"; | |
710 break; | |
711 case TDmode: | |
712 max = "9.999999999999999999999999999999999E6144"; | |
713 break; | |
714 default: | |
715 gcc_unreachable (); | |
716 } | |
717 | |
718 decimal_real_from_string (r, max); | |
719 if (sign) | |
720 decimal128SetSign ((decimal128 *) r->sig, 1); | |
721 } |