1*48fb7bfaSmrg /* Public domain. */ 2*48fb7bfaSmrg typedef int SItype __attribute__ ((mode (SI))); 3*48fb7bfaSmrg typedef unsigned int USItype __attribute__ ((mode (SI))); 4*48fb7bfaSmrg typedef float SFtype __attribute__ ((mode (SF))); 5*48fb7bfaSmrg 6*48fb7bfaSmrg SFtype __floatunsisf(USItype u)7*48fb7bfaSmrg__floatunsisf (USItype u) 8*48fb7bfaSmrg { 9*48fb7bfaSmrg SItype s = (SItype) u; 10*48fb7bfaSmrg if (s < 0) 11*48fb7bfaSmrg { 12*48fb7bfaSmrg /* As in expand_float, compute (u & 1) | (u >> 1) to ensure 13*48fb7bfaSmrg correct rounding if a nonzero bit is shifted out. */ 14*48fb7bfaSmrg return (SFtype) 2.0 * (SFtype) (SItype) ((u & 1) | (u >> 1)); 15*48fb7bfaSmrg } 16*48fb7bfaSmrg else 17*48fb7bfaSmrg return (SFtype) s; 18*48fb7bfaSmrg } 19