|
|
◆ GenrandExp_sub()
| Sub GenrandExp_sub |
( |
Ret As |
Double, |
|
|
Optional Theta As |
Double = 1 |
|
) |
| |
53 ビット実数指数乱数 (Ziggurat 法) (メルセンヌ・ツイスター) (サブルーチン版)
- 目的
- 指数分布 Exp(θ) に従う指数乱数 (53 ビット実数) を生成する.
メルセンヌ・ツイスターの乱数生成ルーチン GenrandRes53() により一様乱数を生成し, Ziggurat 法により指数乱数に変換する.
- 引数
-
| [out] | Ret | 指数乱数 (53 ビット実数). |
| [in] | Theta | (省略可)
指数分布 Exp(θ) のパラメータθ. (省略時 = 1) |
- 参考文献
- Mersenne Twister Home Page
- 四辻哲章「確率分布乱数生成法」第2版, プレアデス出版, 2019.
|