diff options
Diffstat (limited to 'vendor/libm-0.1.4/src/math/modf.rs')
-rw-r--r-- | vendor/libm-0.1.4/src/math/modf.rs | 34 |
1 files changed, 34 insertions, 0 deletions
diff --git a/vendor/libm-0.1.4/src/math/modf.rs b/vendor/libm-0.1.4/src/math/modf.rs new file mode 100644 index 000000000..bcab33a81 --- /dev/null +++ b/vendor/libm-0.1.4/src/math/modf.rs @@ -0,0 +1,34 @@ +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); +} |