OLD | NEW |
1 /**************************************************************** | 1 /**************************************************************** |
2 * | 2 * |
3 * The author of this software is David M. Gay. | 3 * The author of this software is David M. Gay. |
4 * | 4 * |
5 * Copyright (c) 1991, 2000, 2001 by Lucent Technologies. | 5 * Copyright (c) 1991, 2000, 2001 by Lucent Technologies. |
6 * | 6 * |
7 * Permission to use, copy, modify, and distribute this software for any | 7 * Permission to use, copy, modify, and distribute this software for any |
8 * purpose without fee is hereby granted, provided that this entire notice | 8 * purpose without fee is hereby granted, provided that this entire notice |
9 * is included in all copies of any software which is or includes a copy | 9 * is included in all copies of any software which is or includes a copy |
10 * or modification of this software and in all copies of the supporting | 10 * or modification of this software and in all copies of the supporting |
(...skipping 3551 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
3562 */ | 3562 */ |
3563 | 3563 |
3564 int bbits, b2, b5, be, dig, i, ieps, ilim, ilim0, ilim1, | 3564 int bbits, b2, b5, be, dig, i, ieps, ilim, ilim0, ilim1, |
3565 j, j1, k, k0, k_check, leftright, m2, m5, s2, s5, | 3565 j, j1, k, k0, k_check, leftright, m2, m5, s2, s5, |
3566 spec_case, try_quick; | 3566 spec_case, try_quick; |
3567 Long L; | 3567 Long L; |
3568 #ifndef Sudden_Underflow | 3568 #ifndef Sudden_Underflow |
3569 int denorm; | 3569 int denorm; |
3570 ULong x; | 3570 ULong x; |
3571 #endif | 3571 #endif |
3572 » Bigint *b, *b1, *delta, *mlo, *mhi, *S; | 3572 » Bigint *b, *b1, *delta, *mlo = 0, *mhi, *S; // clang (no idea if correc
t) |
3573 U d2, eps, u; | 3573 U d2, eps, u; |
3574 double ds; | 3574 double ds; |
3575 char *s, *s0; | 3575 char *s, *s0; |
3576 #ifdef SET_INEXACT | 3576 #ifdef SET_INEXACT |
3577 int inexact, oldinexact; | 3577 int inexact, oldinexact; |
3578 #endif | 3578 #endif |
3579 #ifdef Honor_FLT_ROUNDS /*{*/ | 3579 #ifdef Honor_FLT_ROUNDS /*{*/ |
3580 int Rounding; | 3580 int Rounding; |
3581 #ifdef Trust_FLT_ROUNDS /*{{ only define this if FLT_ROUNDS really works! */ | 3581 #ifdef Trust_FLT_ROUNDS /*{{ only define this if FLT_ROUNDS really works! */ |
3582 Rounding = Flt_Rounds; | 3582 Rounding = Flt_Rounds; |
(...skipping 617 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
4200 #endif | 4200 #endif |
4201 Bfree(b); | 4201 Bfree(b); |
4202 *s = 0; | 4202 *s = 0; |
4203 *decpt = k + 1; | 4203 *decpt = k + 1; |
4204 if (rve) | 4204 if (rve) |
4205 *rve = s; | 4205 *rve = s; |
4206 return s0; | 4206 return s0; |
4207 } | 4207 } |
4208 | 4208 |
4209 } // namespace dmg_fp | 4209 } // namespace dmg_fp |
OLD | NEW |