[0.0, 1.0) の乱数を得るための“本当の”方法

レイトレ合宿９(*)のセミナー発表スライドです。 * https://sites.google.com/view/rtcamp9/home - 2023/09/08 “除算法2”追記。（@Reputelessさんありがとうございました）

[t-tanaka]

結論は，乱数を扱う用途でdoubleやfloatを用いてはならない，じゃないか？ ほとんどの問題は整数を使った解法に変換できるはず。

[hase0510]

クソリプなんすけど浮動小数じゃなくて浮動小数点数なんすよ。小数点の位置が浮動。floating-point number。

[yorkfield]

exponentを決定するときに0付近の選択確率は下げる話が44ページ前後にある。

[rryu]

浮動小数点数は表せる値の密度が均等ではないということが効いて思っていたよりめんどくさい話だった。

[altar]

表現可能な全ての値を出現させつつ一様にできるけど、厳密にやると最悪(最小値)のケースで2の指数部乗の乱数bitが必要という話。期待値は1/2の等比級数で高々2bitなので、生成時間の秘匿とか考えない限り普通に高速では

[Windymelt]

思ってたよりずっと面白かった、奥深い

[kyopeco]

物理的な話だと思ったらそうじゃなかった。

[takilog]

ほえー面白かった

[khwarizmi]

すげえ。力作。ただ、思ったのはこの関数の小さい方の品質がクリティカルに効くのはexp分布なんだろうけど、ナイーブな逆関数法ではない別実装を用意しとくのが賢いって話なのでは？みたいな

[aya_momo]

途中からどうすればいいのか分かった。

[knok]

True random numberを得る話ではなかった

[mohno]

うわっ、めんどくさ。floatなんか使うな……とはいかないのかな、SIMDとかあるし。（現実に困ったんだろうけど）1.0が出てきて困るなら捨てる、で終わらせそう。

[n_231]

コンピュータに自動コイントスマシーンを組み込もう

[softstone]

０近辺のすべての表現が出現するようにすると、１近辺では精度不均衡になって特定の値に偏るような気がするけどその辺のフォローがよく分からなかった。半分寝てるので気のせいかもしれない。