将u64伪随机数转换为0-1区间f64时结果异常的问题排查

我正在开发一个项目，需要一个极快的伪随机数生成器。目前已实现xorshift算法，可生成伪随机u64，但需要将其转换为0到1区间内的浮点数。

我参考了相关资料，甚至用了和rand库完全相同的实现方式，但结果却不一致，始终无法得到预期的0-1区间结果。

我的代码实现：

let seeds: [u64; 64] = core::array::from_fn(|i| i as u64);

let bitshift12 = u64x64::splat(12);
let bitshift25 = u64x64::splat(25);
let bitshift27 = u64x64::splat(27);

let bitshift52 = u64x64::splat(52);

let mut random_states = Simd::from(seeds);

random_states ^= random_states >> bitshift12;
random_states ^= random_states << bitshift25;
random_states ^= random_states >> bitshift27;

random_states = random_states | ((u64x64::splat(1023) + u64x64::splat(0)) << bitshift52);

let mut generated = Simd::<f64, 64>::from_bits(random_states);

println!("{:?}", generated);

输出结果：

[1.0, 1.0000000074505808, 1.0000000149011616, 1.0000000223517425, 1.0000000298023235, 1.0000000372529039, ...]

错误分析与修正 #

你的核心问题在于构造f64的指数部分设置错误：

• f64的指数偏移量是1023，你设置指数为1023，等价于2^(1023-1023) = 2^0 = 1，再加上随机尾数，最终得到的是1.xxxx格式的浮点数，自然会大于等于1.0。

快速修正（得到[0.5, 1.0)区间）

把指数部分改为1023 - 1 = 1022，对应2^-1 = 0.5，这样生成的数值范围就是[0.5, 1.0)：

// 将原来的1023替换为1022
random_states = random_states | (u64x64::splat(1022) << bitshift52);
let mut generated = Simd::<f64, 64>::from_bits(random_states);

扩展到[0.0, 1.0)区间

如果需要覆盖完整的0到1区间，可通过额外的随机位调整：

// 先生成[0.5, 1.0)的浮点数
random_states = random_states | (u64x64::splat(1022) << bitshift52);
let mut generated = Simd::<f64, 64>::from_bits(random_states);
// 取随机数的最低位，决定是否减去0.5，得到[0.0, 0.5)或[0.5, 1.0)
let offset = Simd::from_bits((random_states & u64x64::splat(1)) << bitshift52);
generated = generated - offset;

对应rand库的标准实现（直接得到[0.0, 1.0)）

取随机数的高53位作为尾数，将指数设为1023 - 53 = 970，这是rand库中Open01分布的实现方式：

// 截取高53位作为尾数，设置对应指数
random_states = (random_states >> (64 - 53)) | (u64x64::splat(970) << 52);
let generated = Simd::<f64, 64>::from_bits(random_states);

另外补充：你当前使用连续整数作为种子，初始生成的随机数会有明显线性趋势，建议使用更分散的随机种子。

内容的提问来源于stack exchange，提问作者BlueOyster