drew
/
acadia-newlib
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

2006-08-03 J"orn Rennecke <joern.rennecke@st.com> 

* dtoa.c (dtoa_r): Handle the denormalized _DOUBLE_IS_32BITS case
        properly.
This commit is contained in:
Jeff Johnston 2006-08-03 16:56:05 +00:00
parent fb43e6237c
commit 079b45dc00
2 changed files with 13 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
newlib/ChangeLog
View File

@ -1,3 +1,8 @@
2006-08-03 J"orn Rennecke <joern.rennecke@st.com>
* dtoa.c (dtoa_r): Handle the denormalized _DOUBLE_IS_32BITS case
properly.
2006-08-01 Jeff Johnston <jjohnstn@redhat.com> 2006-08-01 Jeff Johnston <jjohnstn@redhat.com>
* libc/Makefile.am: Add libc_TEXINFOS dependency on * libc/Makefile.am: Add libc_TEXINFOS dependency on

8
newlib/libc/stdlib/dtoa.c
View File

@ -338,15 +338,23 @@ _DEFUN (_dtoa_r,
/* d is denormalized */ /* d is denormalized */
i = bbits + be + (Bias + (P - 1) - 1); i = bbits + be + (Bias + (P - 1) - 1);
#if defined (_DOUBLE_IS_32BITS)
x = word0 (d) << (32 - i);
#else
x = (i > 32) ? (word0 (d) << (64 - i)) | (word1 (d) >> (i - 32)) x = (i > 32) ? (word0 (d) << (64 - i)) | (word1 (d) >> (i - 32))
: (word1 (d) << (32 - i)); : (word1 (d) << (32 - i));
#endif
d2.d = x; d2.d = x;
word0 (d2) -= 31 * Exp_msk1; /* adjust exponent */ word0 (d2) -= 31 * Exp_msk1; /* adjust exponent */
i -= (Bias + (P - 1) - 1) + 1; i -= (Bias + (P - 1) - 1) + 1;
denorm = 1; denorm = 1;
} }
#endif #endif
#if defined (_DOUBLE_IS_32BITS)
ds = (d2.d - 1.5) * 0.289529651 + 0.176091269 + i * 0.30103001;
#else
ds = (d2.d - 1.5) * 0.289529654602168 + 0.1760912590558 + i * 0.301029995663981; ds = (d2.d - 1.5) * 0.289529654602168 + 0.1760912590558 + i * 0.301029995663981;
#endif
k = (int) ds; k = (int) ds;
if (ds < 0. && ds != k) if (ds < 0. && ds != k)
k--; /* want k = floor(ds) */ k--; /* want k = floor(ds) */