| pub fn modf(x: f64) -> (f64, f64) { |
| let rv2: f64; |
| let mut u = x.to_bits(); |
| let mask: u64; |
| let e = ((u >> 52 & 0x7ff) as i32) - 0x3ff; |
| |
| /* no fractional part */ |
| if e >= 52 { |
| rv2 = x; |
| if e == 0x400 && (u << 12) != 0 { |
| /* nan */ |
| return (x, rv2); |
| } |
| u &= 1 << 63; |
| return (f64::from_bits(u), rv2); |
| } |
| |
| /* no integral part*/ |
| if e < 0 { |
| u &= 1 << 63; |
| rv2 = f64::from_bits(u); |
| return (x, rv2); |
| } |
| |
| mask = ((!0) >> 12) >> e; |
| if (u & mask) == 0 { |
| rv2 = x; |
| u &= 1 << 63; |
| return (f64::from_bits(u), rv2); |
| } |
| u &= !mask; |
| rv2 = f64::from_bits(u); |
| return (x - rv2, rv2); |
| } |
