Current Path : /usr/src/contrib/gcc/config/ |
FreeBSD hs32.drive.ne.jp 9.1-RELEASE FreeBSD 9.1-RELEASE #1: Wed Jan 14 12:18:08 JST 2015 root@hs32.drive.ne.jp:/sys/amd64/compile/hs32 amd64 |
Current File : //usr/src/contrib/gcc/config/floatdisf.c |
/* Public domain. */ typedef int DItype __attribute__ ((mode (DI))); typedef unsigned int UDItype __attribute__ ((mode (DI))); typedef int SItype __attribute__ ((mode (SI))); typedef unsigned int USItype __attribute__ ((mode (SI))); typedef float SFtype __attribute__ ((mode (SF))); typedef float DFtype __attribute__ ((mode (DF))); SFtype __floatdisf (DItype); SFtype __floatdisf (DItype u) { /* Protect against double-rounding error. Represent any low-order bits, that might be truncated by a bit that won't be lost. The bit can go in anywhere below the rounding position of SFtype. A fixed mask and bit position handles all usual configurations. */ if (53 < (sizeof (DItype) * 8) && 53 > ((sizeof (DItype) * 8) - 53 + 24)) { if (!(- ((DItype) 1 << 53) < u && u < ((DItype) 1 << 53))) { if ((UDItype) u & (((UDItype) 1 << (sizeof (DItype) * 8 - 53)) - 1)) { u &= ~ (((UDItype) 1 << (sizeof (DItype) * 8 - 53)) - 1); u |= (UDItype) 1 << (sizeof (DItype) * 8 - 53); } } } /* Do the calculation in a wider type so that we don't lose any of the precision of the high word while multiplying it. */ DFtype f = (SItype) (u >> (sizeof (SItype) * 8)); f *= 0x1p32f; f += (USItype) u; return (SFtype) f; }