o5r/plan9port
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

fix nan64

Browse source
This commit is contained in:
rsc 2007-05-10 04:18:22 +00:00
parent e54f9a4ad2
commit 35920e96a1
1 changed files with 11 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

15
src/lib9/fmt/nan64.c
View file

@ -26,11 +26,18 @@ __NaN(void)
int int
__isNaN(double d) __isNaN(double d)
{ {
/*
* Used to just say x = *(uvlong*)&d,
* but gcc miscompiles that!
*/
union {
uvlong i;
double f;
} u;
uvlong x; uvlong x;
double *p;
u.f = d;
p = &d; x = u.i;
x = *(uvlong*)p;
/* IEEE 754: exponent bits 0x7FF and non-zero mantissa */ /* IEEE 754: exponent bits 0x7FF and non-zero mantissa */
return (x&uvinf) == uvinf && (x&~uvneginf) != 0; return (x&uvinf) == uvinf && (x&~uvneginf) != 0;
} }