1 /* Copyright (C) 2007-2023 Free Software Foundation, Inc.
2
3 This file is part of GCC.
4
5 GCC is free software; you can redistribute it and/or modify it under
6 the terms of the GNU General Public License as published by the Free
7 Software Foundation; either version 3, or (at your option) any later
8 version.
9
10 GCC is distributed in the hope that it will be useful, but WITHOUT ANY
11 WARRANTY; without even the implied warranty of MERCHANTABILITY or
12 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
13 for more details.
14
15 Under Section 7 of GPL version 3, you are granted additional
16 permissions described in the GCC Runtime Library Exception, version
17 3.1, as published by the Free Software Foundation.
18
19 You should have received a copy of the GNU General Public License and
20 a copy of the GCC Runtime Library Exception along with this program;
21 see the files COPYING3 and COPYING.RUNTIME respectively. If not, see
22 <http://www.gnu.org/licenses/>. */
23
24 /*****************************************************************************
25 * BID128_to_string
26 ****************************************************************************/
27
28 #define BID_128RES
29 #include <stdio.h>
30 #include "bid_internal.h"
31 #include "bid128_2_str.h"
32 #include "bid128_2_str_macros.h"
33
34 extern int bid128_coeff_2_string (UINT64 X_hi, UINT64 X_lo,
35 char *char_ptr);
36
37 #if DECIMAL_CALL_BY_REFERENCE
38
39 void
40 bid128_to_string (char *str,
41 UINT128 *
42 px _EXC_FLAGS_PARAM _EXC_MASKS_PARAM
43 _EXC_INFO_PARAM) {
44 UINT128 x;
45 #else
46
47 void
48 bid128_to_string (char *str, UINT128 x
49 _EXC_FLAGS_PARAM _EXC_MASKS_PARAM _EXC_INFO_PARAM) {
50 #endif
51 UINT64 x_sign;
52 UINT64 x_exp;
53 int exp; // unbiased exponent
54 // Note: C1.w[1], C1.w[0] represent x_signif_hi, x_signif_lo (all are UINT64)
55 int ind;
56 UINT128 C1;
57 unsigned int k = 0; // pointer in the string
58 unsigned int d0, d123;
59 UINT64 HI_18Dig, LO_18Dig, Tmp;
60 UINT32 MiDi[12], *ptr;
61 char *c_ptr_start, *c_ptr;
62 int midi_ind, k_lcv, len;
63
64 #if DECIMAL_CALL_BY_REFERENCE
65 x = *px;
66 #endif
67
68 BID_SWAP128(x);
69 // check for NaN or Infinity
70 if ((x.w[1] & MASK_SPECIAL) == MASK_SPECIAL) {
71 // x is special
72 if ((x.w[1] & MASK_NAN) == MASK_NAN) { // x is NAN
73 if ((x.w[1] & MASK_SNAN) == MASK_SNAN) { // x is SNAN
74 // set invalid flag
75 str[0] = ((SINT64)x.w[1]<0)? '-':'+';
76 str[1] = 'S';
77 str[2] = 'N';
78 str[3] = 'a';
79 str[4] = 'N';
80 str[5] = '\0';
81 } else { // x is QNaN
82 str[0] = ((SINT64)x.w[1]<0)? '-':'+';
83 str[1] = 'Q';
84 str[2] = 'N';
85 str[3] = 'a';
86 str[4] = 'N';
87 str[5] = '\0';
88 }
89 } else { // x is not a NaN, so it must be infinity
90 if ((x.w[1] & MASK_SIGN) == 0x0ull) { // x is +inf
91 str[0] = '+';
92 str[1] = 'I';
93 str[2] = 'n';
94 str[3] = 'f';
95 str[4] = '\0';
96 } else { // x is -inf
97 str[0] = '-';
98 str[1] = 'I';
99 str[2] = 'n';
100 str[3] = 'f';
101 str[4] = '\0';
102 }
103 }
104 return;
105 } else if (((x.w[1] & MASK_COEFF) == 0x0ull) && (x.w[0] == 0x0ull)) {
106 // x is 0
107 len = 0;
108
109 //determine if +/-
110 if (x.w[1] & MASK_SIGN)
111 str[len++] = '-';
112 else
113 str[len++] = '+';
114 str[len++] = '0';
115 str[len++] = 'E';
116
117 // extract the exponent and print
118 exp = (int) (((x.w[1] & MASK_EXP) >> 49) - 6176);
119 if(exp > (((0x5ffe)>>1) - (6176))) {
120 exp = (int) ((((x.w[1]<<2) & MASK_EXP) >> 49) - 6176);
121 }
122 if (exp >= 0) {
123 str[len++] = '+';
124 len += sprintf (str + len, "%u", exp);// should not use sprintf (should
125 // use sophisticated algorithm, since we know range of exp is limited)
126 str[len++] = '\0';
127 } else {
128 len += sprintf (str + len, "%d", exp);// should not use sprintf (should
129 // use sophisticated algorithm, since we know range of exp is limited)
130 str[len++] = '\0';
131 }
132 return;
133 } else { // x is not special and is not zero
134 // unpack x
135 x_sign = x.w[1] & MASK_SIGN;// 0 for positive, MASK_SIGN for negative
136 x_exp = x.w[1] & MASK_EXP;// biased and shifted left 49 bit positions
137 if ((x.w[1] & 0x6000000000000000ull) == 0x6000000000000000ull)
138 x_exp = (x.w[1]<<2) & MASK_EXP;// biased and shifted left 49 bit positions
139 C1.w[1] = x.w[1] & MASK_COEFF;
140 C1.w[0] = x.w[0];
141 exp = (x_exp >> 49) - 6176;
142
143 // determine sign's representation as a char
144 if (x_sign)
145 str[k++] = '-';// negative number
146 else
147 str[k++] = '+';// positive number
148
149 // determine coefficient's representation as a decimal string
150
151 // if zero or non-canonical, set coefficient to '0'
152 if ((C1.w[1] > 0x0001ed09bead87c0ull) ||
153 (C1.w[1] == 0x0001ed09bead87c0ull &&
154 (C1.w[0] > 0x378d8e63ffffffffull)) ||
155 ((x.w[1] & 0x6000000000000000ull) == 0x6000000000000000ull) ||
156 ((C1.w[1] == 0) && (C1.w[0] == 0))) {
157 str[k++] = '0';
158 } else {
159 /* ****************************************************
160 This takes a bid coefficient in C1.w[1],C1.w[0]
161 and put the converted character sequence at location
162 starting at &(str[k]). The function returns the number
163 of MiDi returned. Note that the character sequence
164 does not have leading zeros EXCEPT when the input is of
165 zero value. It will then output 1 character '0'
166 The algorithm essentailly tries first to get a sequence of
167 Millenial Digits "MiDi" and then uses table lookup to get the
168 character strings of these MiDis.
169 **************************************************** */
170 /* Algorithm first decompose possibly 34 digits in hi and lo
171 18 digits. (The high can have at most 16 digits). It then
172 uses macro that handle 18 digit portions.
173 The first step is to get hi and lo such that
174 2^(64) C1.w[1] + C1.w[0] = hi * 10^18 + lo, 0 <= lo < 10^18.
175 We use a table lookup method to obtain the hi and lo 18 digits.
176 [C1.w[1],C1.w[0]] = c_8 2^(107) + c_7 2^(101) + ... + c_0 2^(59) + d
177 where 0 <= d < 2^59 and each c_j has 6 bits. Because d fits in
178 18 digits, we set hi = 0, and lo = d to begin with.
179 We then retrieve from a table, for j = 0, 1, ..., 8
180 that gives us A and B where c_j 2^(59+6j) = A * 10^18 + B.
181 hi += A ; lo += B; After each accumulation into lo, we normalize
182 immediately. So at the end, we have the decomposition as we need. */
183
184 Tmp = C1.w[0] >> 59;
185 LO_18Dig = (C1.w[0] << 5) >> 5;
186 Tmp += (C1.w[1] << 5);
187 HI_18Dig = 0;
188 k_lcv = 0;
189 // Tmp = {C1.w[1]{49:0}, C1.w[0]{63:59}}
190 // Lo_18Dig = {C1.w[0]{58:0}}
191
192 while (Tmp) {
193 midi_ind = (int) (Tmp & 0x000000000000003FLL);
194 midi_ind <<= 1;
195 Tmp >>= 6;
196 HI_18Dig += mod10_18_tbl[k_lcv][midi_ind++];
197 LO_18Dig += mod10_18_tbl[k_lcv++][midi_ind];
198 __L0_Normalize_10to18 (HI_18Dig, LO_18Dig);
199 }
200 ptr = MiDi;
201 if (HI_18Dig == 0LL) {
202 __L1_Split_MiDi_6_Lead (LO_18Dig, ptr);
203 } else {
204 __L1_Split_MiDi_6_Lead (HI_18Dig, ptr);
205 __L1_Split_MiDi_6 (LO_18Dig, ptr);
206 }
207 len = ptr - MiDi;
208 c_ptr_start = &(str[k]);
209 c_ptr = c_ptr_start;
210
211 /* now convert the MiDi into character strings */
212 __L0_MiDi2Str_Lead (MiDi[0], c_ptr);
213 for (k_lcv = 1; k_lcv < len; k_lcv++) {
214 __L0_MiDi2Str (MiDi[k_lcv], c_ptr);
215 }
216 k = k + (c_ptr - c_ptr_start);
217 }
218
219 // print E and sign of exponent
220 str[k++] = 'E';
221 if (exp < 0) {
222 exp = -exp;
223 str[k++] = '-';
224 } else {
225 str[k++] = '+';
226 }
227
228 // determine exponent's representation as a decimal string
229 // d0 = exp / 1000;
230 // Use Property 1
231 d0 = (exp * 0x418a) >> 24;// 0x418a * 2^-24 = (10^(-3))RP,15
232 d123 = exp - 1000 * d0;
233
234 if (d0) { // 1000 <= exp <= 6144 => 4 digits to return
235 str[k++] = d0 + 0x30;// ASCII for decimal digit d0
236 ind = 3 * d123;
237 str[k++] = char_table3[ind];
238 str[k++] = char_table3[ind + 1];
239 str[k++] = char_table3[ind + 2];
240 } else { // 0 <= exp <= 999 => d0 = 0
241 if (d123 < 10) { // 0 <= exp <= 9 => 1 digit to return
242 str[k++] = d123 + 0x30;// ASCII
243 } else if (d123 < 100) { // 10 <= exp <= 99 => 2 digits to return
244 ind = 2 * (d123 - 10);
245 str[k++] = char_table2[ind];
246 str[k++] = char_table2[ind + 1];
247 } else { // 100 <= exp <= 999 => 3 digits to return
248 ind = 3 * d123;
249 str[k++] = char_table3[ind];
250 str[k++] = char_table3[ind + 1];
251 str[k++] = char_table3[ind + 2];
252 }
253 }
254 str[k] = '\0';
255
256 }
257 return;
258
259 }
260
261
262 #define MAX_FORMAT_DIGITS_128 34
263 #define MAX_STRING_DIGITS_128 100
264 #define MAX_SEARCH MAX_STRING_DIGITS_128-MAX_FORMAT_DIGITS_128-1
265
266
267 #if DECIMAL_CALL_BY_REFERENCE
268
269 void
270 bid128_from_string (UINT128 * pres,
271 char *ps _RND_MODE_PARAM _EXC_FLAGS_PARAM
272 _EXC_MASKS_PARAM _EXC_INFO_PARAM) {
273 #else
274
275 UINT128
276 bid128_from_string (char *ps _RND_MODE_PARAM _EXC_FLAGS_PARAM
277 _EXC_MASKS_PARAM _EXC_INFO_PARAM) {
278 #endif
279 UINT128 CX, res;
280 UINT64 sign_x, coeff_high, coeff_low, coeff2, coeff_l2, carry = 0x0ull,
281 scale_high, right_radix_leading_zeros;
282 int ndigits_before, ndigits_after, ndigits_total, dec_expon, sgn_exp,
283 i, d2, rdx_pt_enc;
284 char c, buffer[MAX_STRING_DIGITS_128];
285 int save_rnd_mode;
286 int save_fpsf;
287
288 #if DECIMAL_CALL_BY_REFERENCE
289 #if !DECIMAL_GLOBAL_ROUNDING
290 _IDEC_round rnd_mode = *prnd_mode;
291 #endif
292 #endif
293
294 save_rnd_mode = rnd_mode; // dummy
295 save_fpsf = *pfpsf; // dummy
296
297 right_radix_leading_zeros = rdx_pt_enc = 0;
298
299 // if null string, return NaN
300 if (!ps) {
301 res.w[1] = 0x7c00000000000000ull;
302 res.w[0] = 0;
303 BID_RETURN (res);
304 }
305 // eliminate leading white space
306 while ((*ps == ' ') || (*ps == '\t'))
307 ps++;
308
309 // c gets first character
310 c = *ps;
311
312
313 // if c is null or not equal to a (radix point, negative sign,
314 // positive sign, or number) it might be SNaN, sNaN, Infinity
315 if (!c
316 || (c != '.' && c != '-' && c != '+'
317 && ((unsigned) (c - '0') > 9))) {
318 res.w[0] = 0;
319 // Infinity?
320 if ((tolower_macro (ps[0]) == 'i' && tolower_macro (ps[1]) == 'n'
321 && tolower_macro (ps[2]) == 'f')
322 && (!ps[3]
323 || (tolower_macro (ps[3]) == 'i'
324 && tolower_macro (ps[4]) == 'n'
325 && tolower_macro (ps[5]) == 'i'
326 && tolower_macro (ps[6]) == 't'
327 && tolower_macro (ps[7]) == 'y' && !ps[8])
328 )) {
329 res.w[1] = 0x7800000000000000ull;
330 BID_RETURN (res);
331 }
332 // return sNaN
333 if (tolower_macro (ps[0]) == 's' && tolower_macro (ps[1]) == 'n' &&
334 tolower_macro (ps[2]) == 'a' && tolower_macro (ps[3]) == 'n') {
335 // case insensitive check for snan
336 res.w[1] = 0x7e00000000000000ull;
337 BID_RETURN (res);
338 } else {
339 // return qNaN
340 res.w[1] = 0x7c00000000000000ull;
341 BID_RETURN (res);
342 }
343 }
344 // if +Inf, -Inf, +Infinity, or -Infinity (case insensitive check for inf)
345 if ((tolower_macro (ps[1]) == 'i' && tolower_macro (ps[2]) == 'n' &&
346 tolower_macro (ps[3]) == 'f') && (!ps[4] ||
347 (tolower_macro (ps[4]) == 'i' && tolower_macro (ps[5]) == 'n' &&
348 tolower_macro (ps[6]) == 'i' && tolower_macro (ps[7]) == 't' &&
349 tolower_macro (ps[8]) == 'y' && !ps[9]))) { // ci check for infinity
350 res.w[0] = 0;
351
352 if (c == '+')
353 res.w[1] = 0x7800000000000000ull;
354 else if (c == '-')
355 res.w[1] = 0xf800000000000000ull;
356 else
357 res.w[1] = 0x7c00000000000000ull;
358
359 BID_RETURN (res);
360 }
361 // if +sNaN, +SNaN, -sNaN, or -SNaN
362 if (tolower_macro (ps[1]) == 's' && tolower_macro (ps[2]) == 'n'
363 && tolower_macro (ps[3]) == 'a' && tolower_macro (ps[4]) == 'n') {
364 res.w[0] = 0;
365 if (c == '-')
366 res.w[1] = 0xfe00000000000000ull;
367 else
368 res.w[1] = 0x7e00000000000000ull;
369 BID_RETURN (res);
370 }
371 // set up sign_x to be OR'ed with the upper word later
372 if (c == '-')
373 sign_x = 0x8000000000000000ull;
374 else
375 sign_x = 0;
376
377 // go to next character if leading sign
378 if (c == '-' || c == '+')
379 ps++;
380
381 c = *ps;
382
383 // if c isn't a decimal point or a decimal digit, return NaN
384 if (c != '.' && ((unsigned) (c - '0') > 9)) {
385 res.w[1] = 0x7c00000000000000ull | sign_x;
386 res.w[0] = 0;
387 BID_RETURN (res);
388 }
389 // detect zero (and eliminate/ignore leading zeros)
390 if (*(ps) == '0') {
391
392 // if all numbers are zeros (with possibly 1 radix point, the number is zero
393 // should catch cases such as: 000.0
394 while (*ps == '0') {
395
396 ps++;
397
398 // for numbers such as 0.0000000000000000000000000000000000001001,
399 // we want to count the leading zeros
400 if (rdx_pt_enc) {
401 right_radix_leading_zeros++;
402 }
403 // if this character is a radix point, make sure we haven't already
404 // encountered one
405 if (*(ps) == '.') {
406 if (rdx_pt_enc == 0) {
407 rdx_pt_enc = 1;
408 // if this is the first radix point, and the next character is NULL,
409 // we have a zero
410 if (!*(ps + 1)) {
411 res.w[1] =
412 (0x3040000000000000ull -
413 (right_radix_leading_zeros << 49)) | sign_x;
414 res.w[0] = 0;
415 BID_RETURN (res);
416 }
417 ps = ps + 1;
418 } else {
419 // if 2 radix points, return NaN
420 res.w[1] = 0x7c00000000000000ull | sign_x;
421 res.w[0] = 0;
422 BID_RETURN (res);
423 }
424 } else if (!*(ps)) {
425 //res.w[1] = 0x3040000000000000ull | sign_x;
426 res.w[1] =
427 (0x3040000000000000ull -
428 (right_radix_leading_zeros << 49)) | sign_x;
429 res.w[0] = 0;
430 BID_RETURN (res);
431 }
432 }
433 }
434
435 c = *ps;
436
437 // initialize local variables
438 ndigits_before = ndigits_after = ndigits_total = 0;
439 sgn_exp = 0;
440 // pstart_coefficient = ps;
441
442 if (!rdx_pt_enc) {
443 // investigate string (before radix point)
444 while ((unsigned) (c - '0') <= 9
445 && ndigits_before < MAX_STRING_DIGITS_128) {
446 buffer[ndigits_before] = c;
447 ps++;
448 c = *ps;
449 ndigits_before++;
450 }
451
452 ndigits_total = ndigits_before;
453 if (c == '.') {
454 ps++;
455 if ((c = *ps)) {
456
457 // investigate string (after radix point)
458 while ((unsigned) (c - '0') <= 9
459 && ndigits_total < MAX_STRING_DIGITS_128) {
460 buffer[ndigits_total] = c;
461 ps++;
462 c = *ps;
463 ndigits_total++;
464 }
465 ndigits_after = ndigits_total - ndigits_before;
466 }
467 }
468 } else {
469 // we encountered a radix point while detecting zeros
470 //if (c = *ps){
471
472 c = *ps;
473 ndigits_total = 0;
474 // investigate string (after radix point)
475 while ((unsigned) (c - '0') <= 9
476 && ndigits_total < MAX_STRING_DIGITS_128) {
477 buffer[ndigits_total] = c;
478 ps++;
479 c = *ps;
480 ndigits_total++;
481 }
482 ndigits_after = ndigits_total - ndigits_before;
483 }
484
485 // get exponent
486 dec_expon = 0;
487 if (ndigits_total < MAX_STRING_DIGITS_128) {
488 if (c) {
489 if (c != 'e' && c != 'E') {
490 // return NaN
491 res.w[1] = 0x7c00000000000000ull;
492 res.w[0] = 0;
493 BID_RETURN (res);
494 }
495 ps++;
496 c = *ps;
497
498 if (((unsigned) (c - '0') > 9)
499 && ((c != '+' && c != '-') || (unsigned) (ps[1] - '0') > 9)) {
500 // return NaN
501 res.w[1] = 0x7c00000000000000ull;
502 res.w[0] = 0;
503 BID_RETURN (res);
504 }
505
506 if (c == '-') {
507 sgn_exp = -1;
508 ps++;
509 c = *ps;
510 } else if (c == '+') {
511 ps++;
512 c = *ps;
513 }
514
515 dec_expon = c - '0';
516 i = 1;
517 ps++;
518 c = *ps - '0';
519 while (((unsigned) c) <= 9 && i < 7) {
520 d2 = dec_expon + dec_expon;
521 dec_expon = (d2 << 2) + d2 + c;
522 ps++;
523 c = *ps - '0';
524 i++;
525 }
526 }
527
528 dec_expon = (dec_expon + sgn_exp) ^ sgn_exp;
529 }
530
531
532 if (ndigits_total <= MAX_FORMAT_DIGITS_128) {
533 dec_expon +=
534 DECIMAL_EXPONENT_BIAS_128 - ndigits_after -
535 right_radix_leading_zeros;
536 if (dec_expon < 0) {
537 res.w[1] = 0 | sign_x;
538 res.w[0] = 0;
539 }
540 if (ndigits_total == 0) {
541 CX.w[0] = 0;
542 CX.w[1] = 0;
543 } else if (ndigits_total <= 19) {
544 coeff_high = buffer[0] - '0';
545 for (i = 1; i < ndigits_total; i++) {
546 coeff2 = coeff_high + coeff_high;
547 coeff_high = (coeff2 << 2) + coeff2 + buffer[i] - '0';
548 }
549 CX.w[0] = coeff_high;
550 CX.w[1] = 0;
551 } else {
552 coeff_high = buffer[0] - '0';
553 for (i = 1; i < ndigits_total - 17; i++) {
554 coeff2 = coeff_high + coeff_high;
555 coeff_high = (coeff2 << 2) + coeff2 + buffer[i] - '0';
556 }
557 coeff_low = buffer[i] - '0';
558 i++;
559 for (; i < ndigits_total; i++) {
560 coeff_l2 = coeff_low + coeff_low;
561 coeff_low = (coeff_l2 << 2) + coeff_l2 + buffer[i] - '0';
562 }
563 // now form the coefficient as coeff_high*10^19+coeff_low+carry
564 scale_high = 100000000000000000ull;
565 __mul_64x64_to_128_fast (CX, coeff_high, scale_high);
566
567 CX.w[0] += coeff_low;
568 if (CX.w[0] < coeff_low)
569 CX.w[1]++;
570 }
571 get_BID128 (&res, sign_x, dec_expon, CX,&rnd_mode,pfpsf);
572 BID_RETURN (res);
573 } else {
574 // simply round using the digits that were read
575
576 dec_expon +=
577 ndigits_before + DECIMAL_EXPONENT_BIAS_128 -
578 MAX_FORMAT_DIGITS_128 - right_radix_leading_zeros;
579
580 if (dec_expon < 0) {
581 res.w[1] = 0 | sign_x;
582 res.w[0] = 0;
583 }
584
585 coeff_high = buffer[0] - '0';
586 for (i = 1; i < MAX_FORMAT_DIGITS_128 - 17; i++) {
587 coeff2 = coeff_high + coeff_high;
588 coeff_high = (coeff2 << 2) + coeff2 + buffer[i] - '0';
589 }
590 coeff_low = buffer[i] - '0';
591 i++;
592 for (; i < MAX_FORMAT_DIGITS_128; i++) {
593 coeff_l2 = coeff_low + coeff_low;
594 coeff_low = (coeff_l2 << 2) + coeff_l2 + buffer[i] - '0';
595 }
596 switch(rnd_mode) {
597 case ROUNDING_TO_NEAREST:
598 carry = ((unsigned) ('4' - buffer[i])) >> 31;
599 if ((buffer[i] == '5' && !(coeff_low & 1)) || dec_expon < 0) {
600 if (dec_expon >= 0) {
601 carry = 0;
602 i++;
603 }
604 for (; i < ndigits_total; i++) {
605 if (buffer[i] > '0') {
606 carry = 1;
607 break;
608 }
609 }
610 }
611 break;
612
613 case ROUNDING_DOWN:
614 if(sign_x)
615 for (; i < ndigits_total; i++) {
616 if (buffer[i] > '0') {
617 carry = 1;
618 break;
619 }
620 }
621 break;
622 case ROUNDING_UP:
623 if(!sign_x)
624 for (; i < ndigits_total; i++) {
625 if (buffer[i] > '0') {
626 carry = 1;
627 break;
628 }
629 }
630 break;
631 case ROUNDING_TO_ZERO:
632 carry=0;
633 break;
634 case ROUNDING_TIES_AWAY:
635 carry = ((unsigned) ('4' - buffer[i])) >> 31;
636 if (dec_expon < 0) {
637 for (; i < ndigits_total; i++) {
638 if (buffer[i] > '0') {
639 carry = 1;
640 break;
641 }
642 }
643 }
644 break;
645
646
647 }
648 // now form the coefficient as coeff_high*10^17+coeff_low+carry
649 scale_high = 100000000000000000ull;
650 if (dec_expon < 0) {
651 if (dec_expon > -MAX_FORMAT_DIGITS_128) {
652 scale_high = 1000000000000000000ull;
653 coeff_low = (coeff_low << 3) + (coeff_low << 1);
654 dec_expon--;
655 }
656 if (dec_expon == -MAX_FORMAT_DIGITS_128
657 && coeff_high > 50000000000000000ull)
658 carry = 0;
659 }
660
661 __mul_64x64_to_128_fast (CX, coeff_high, scale_high);
662
663 coeff_low += carry;
664 CX.w[0] += coeff_low;
665 if (CX.w[0] < coeff_low)
666 CX.w[1]++;
667
668
669 get_BID128(&res, sign_x, dec_expon, CX, &rnd_mode, pfpsf);
670 BID_RETURN (res);
671 }
672 }