spec_math/utils/
mod.rs

1pub fn frexp(x: f64) -> (f64, i32) {
2    let mut y = x.to_bits();
3    let ee = ((y >> 52) & 0x7ff) as i32;
4
5    if ee == 0 {
6        if x != 0.0 {
7            let x1p64 = f64::from_bits(0x43f0000000000000);
8            let (x, e) = frexp(x * x1p64);
9            (x, e - 64)
10        } else {
11            (x, 0)
12        }
13    } else if ee == 0x7ff {
14        (x, 0)
15    } else {
16        let e = ee - 0x3fe;
17        y &= 0x800fffffffffffff;
18        y |= 0x3fe0000000000000;
19        
20        (f64::from_bits(y), e)
21    }
22}