メルセンヌツイスタの高速版（SFMT)がExcel VBAで簡単に使える！

いつの間にか、メルセンヌツイスタの高速版（SIMD-oriented Fast Mersenne Twister (SFMT)）が開発されていた！2倍速いとのこと。

http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/SFMT/index-jp.html

しかも、倍精度浮動小数点に特化したDouble precision SIMD-oriented Fast Mersenne    Twister (dSFMT)も！私がしたいシミュレーション、乱数は正規乱数が多いのでこれはいい。

このｄSFMTを和田さんがdll化、VBAその他で使えるようにしてくれています。すばらしい！

ありがとうございます。

http://www001.upp.so-net.ne.jp/isaku/

ここからたどっていって

http://www001.upp.so-net.ne.jp/isaku/rand2.html

のlibZD.zipです。

では使わせてもらいます。

zipを解凍して、libZD.dllをパスが通った場所に置く。例えばc:\windows。

で、まずは正規乱数を30000個生成してみる。こんな感じで。

Option Explicit

Private Sub 実行_Click()
    Dim s As Long
    Dim i As Integer
   
    s = 1
   
    Call InitMt(s)
    Application.ScreenUpdating = False
   
    For i = 1 To 30000
        Worksheets("Sheet1").Cells(i + 2, 3) = NextNormal
    Next i
   
    Application.ScreenUpdating = True
   
End Sub

簡単。で、これをExcelのデータ分析ツールでヒストグラムにすると、、、

きれいな釣鐘型。ガウス乱数がものすごい速さで生成された。これは十分実用的ではないだろうか。これを使っていろいろやってみよう。