rand.rs

//type F = f32;
//type U = u32;
type F = f64;
type U = u64;

pub fn sample_a(i: U) -> F {
    let d = 2.0 / F::EPSILON;
    (i % d as U) as F / d
}

pub fn sample_b(i: U) -> F {
    let b = F::MANTISSA_DIGITS - 1;
    let i = i % (2 << b);
    F::from_bits(match i.checked_ilog2() {
        None => 0,
        Some(e) => U::from(F::MAX_EXP as u32 - b - 2 + e) << b | (i ^ 1 << e) << (b - e),
    })
}