我不想使用 System.Random,因为它的速度要慢得多,而且我需要生成数百万个随机浮点数。我也不能使用 System.Random.MWC 因为它不纯。
我自己尝试编写此代码,但我的解决方案没有生成统一的范围,所有值都非常接近 0。
randomRMSP :: (Float,Float) -> PureMT -> (Float,PureMT)
randomRMSP (lo,hi) rng =
let (f,rng') = first double2Float (randomDouble rng)
in (lo + f / (maxFloat + 1) * (hi - lo + 1),rng')
maxRealFloat :: RealFloat a => a -> a
maxRealFloat a = encodeFloat m n where
b = floatRadix a
e = floatDigits a
(_, e') = floatRange a
m = b ^ e - 1
n = e' - e
我很确定 maxRealFloat 函数是正确的,因为它根据 Wikipedia 返回了 Floats 和 Doubles 的正确值