1

我不想使用 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 的正确值

4

1 回答 1

3

我意识到我犯了一个非常愚蠢的错误。System.Random.Mersenne.Pure64 中的函数 randomDouble 返回一个已经在 [0,1] 范围内的双精度数(文档根本没有提到),因此我所有的数字都非常小,因为我将它们除以最大双精度数.

于 2015-07-20T23:38:19.307 回答