RepoMirrors
/
musl
mirror of git://git.musl-libc.org/musl
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

math: support invalid ld80 representations in fpclassify 

apparently gnulib requires invalid long double representations
to be handled correctly in printf so we classify them according
to how the fpu treats them: bad inf is nan, bad nan is nan,
bad normal is nan and bad subnormal/zero is minimal normal
This commit is contained in:
Szabolcs Nagy 2013-09-05 18:05:07 +00:00
parent f4d9bfb3a9
commit f657fe4b9f
1 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
src/math/__fpclassifyl.c
View File

@ -3,13 +3,15 @@
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
#elif LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384 #elif LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384
/* invalid representations (top bit of u.i.m is wrong) are not handled */
int __fpclassifyl(long double x) int __fpclassifyl(long double x)
{ {
union ldshape u = {x}; union ldshape u = {x};
int e = u.i.se & 0x7fff; int e = u.i.se & 0x7fff;
if (!e) int msb = u.i.m>>63;
if (!e && !msb)
return u.i.m ? FP_SUBNORMAL : FP_ZERO; return u.i.m ? FP_SUBNORMAL : FP_ZERO;
if (!msb)
return FP_NAN;
if (e == 0x7fff) if (e == 0x7fff)
return u.i.m << 1 ? FP_NAN : FP_INFINITE; return u.i.m << 1 ? FP_NAN : FP_INFINITE;
return FP_NORMAL; return FP_NORMAL;