Mersenne Twister in OCaml

擬似乱数生成アルゴリズム Mersenne Twister（メルセンヌ・ツイスタ）を OCaml で使えるライブラリにしました。

ダウンロード

• ダウンロード : mt_ocaml.zip

ソースコード

mt_ocaml.zip には、次の3つのファイルが含まれています。

• mersenne_twister.ml : モジュール Mersenne_twister の定義
• mt_c.c : モジュール Mersenne_twister の実装（C言語による）
• mt_test.ml : テスト用ソースコード（mt19937ar.c と同一の出力を得ます）

使い方

mersenne_twister.ml と mt_c.c をコンパイルしてリンクすることで、モジュール Mersenne_twister を使えるようになります。
モジュール Mersenne_twister では、次の型や関数が定義されています。

type t

Mersenne Twister を表す型です。t 型の値には、Mersenne Twister の内部状態が保存されています。

val init_genrand : int32 -> t

int32 型の整数値をシード値として、新しい Mersenne Twister を生成します。

val init_by_array : int32 array -> t

int32 型の配列をシード値として、新しい Mersenne Twister を生成します。

val init_by_state : int -> int32 array -> t

関数 get_state（後述）で取得した内部状態から、Mersenne Twister を復元します。
関数 get_state によって得られたもの以外を引数として与えた場合は、期待した結果が得られない可能性があります。

val genrand_int : t -> int

Mersenne Twister から新たな int 型の乱数を得ます。
乱数の範囲 : -1,073,741,824 ～ +1,073,741,823 （31ビット整数）

val genrand_int32 : t -> int32

Mersenne Twister から新たな int32 型の乱数を得ます。
乱数の範囲 : -2,147,483,648 ～ +2,147,483,647 （32ビット整数）

val genrand_int64 : t -> int64

Mersenne Twister から新たな int64 型の乱数を得ます。
乱数の範囲 : -9,223,372,036,854,775,808 ～ +9,223,372,036,854,775,807 （64ビット整数）

val genrand_real1 : t -> float

Mersenne Twister から新たな float 型の乱数を得ます。
乱数の範囲 : 0 以上 1 以下 （32ビット精度実数）

val genrand_real2 : t -> float

Mersenne Twister から新たな float 型の乱数を得ます。
乱数の範囲 : 0 以上 1 未満 （32ビット精度実数）

val genrand_real3 : t -> float

Mersenne Twister から新たな float 型の乱数を得ます。
乱数の範囲 : 0 超 1 未満 （32ビット精度実数）

val genrand_res53 : t -> float

Mersenne Twister から新たな float 型の乱数を得ます。
乱数の範囲 : 0 以上 1 未満 （53ビット精度実数）

val copy : t -> t

引数として与えたものと全く同じ内部状態を持つ Mersenne Twister を複製し、返します。

val get_state : t -> int * int32 array

Mersenne Twister の内部状態（int 値と int32 array 値）を得ます。
これを保存しておけば、前述の関数 init_by_state を用いてこの時の状態を復元できます。

使う際は、最初に1度 init_genrand や init_by_array で Mersenne Twister を生成し、以降はそれを genrand_int 等の引数に与えることで乱数を取得してください。

コンパイル・実行の仕方

OCaml 環境はインストール済とします。
以下、Windows + Visual C++ と、Linux + GCC での方法を記します。他の環境でもほぼ同様です。

インタプリタでの使い方

Windows + Visual C++ の場合

Visual Studio コマンドプロンプトで次のようにすると、インタプリタ mt_interp.exe を生成できます。ユーザが打つのは太字の部分です。
生成した mt_interp.exe は普通にインタプリタとして使用でき、かつモジュール Mersenne_twister が使えます。

Linux + GCC の場合

ターミナル画面で次のようにすると、インタプリタ mt_interp を生成できます。ユーザが打つのは太字の部分です。
生成した mt_interp.exe は普通にインタプリタとして使用でき、かつモジュール Mersenne_twister が使えます。

バイトコードコンパイラ ocamlc での使い方

Windows + Visual C++ の場合

モジュール Mersenne_twister を含む OCaml プログラム mt_test.ml を ocamlc でコンパイルするには、Visual Studio コマンドプロンプトで次のようにします。
ユーザが打つのは太字の部分です。

Linux + GCC の場合

モジュール Mersenne_twister を含む OCaml プログラム mt_test.ml を ocamlc でコンパイルするには、ターミナル画面で次のようにします。
ユーザが打つのは太字の部分です。

ネイティブコードコンパイラ ocamlopt での使い方

Windows + Visual C++ の場合

モジュール Mersenne_twister を含む OCaml プログラム mt_test.ml を ocamlopt でコンパイルするには、Visual Studio コマンドプロンプトで次のようにします。
ユーザが打つのは太字の部分です。

Linux + GCC の場合

モジュール Mersenne_twister を含む OCaml プログラム mt_test.ml を ocamlopt でコンパイルするには、ターミナル画面で次のようにします。
ユーザが打つのは太字の部分です。

注意事項

本ライブラリでは mt_c.c の中で stdint.h を include して使用していますが、コンパイラのバージョンによっては stdint.h が利用できない場合があります。
その場合は、mt_c.c 中の #include <stdint.h> を削除し、typedef unsigned long uint32_t のような記述を追加してください。
（環境によっては typedef unsigned int uint32_t としないと正常に動作しない場合もあります）

よくわからない場合は、最新版の Visual C++（Windows の場合）を使うことをお勧めします。
（Visual C++ 2010 では stdint.h が正常に使えることを確認しています）

当サイトが配布しているプログラムを利用したことにより生じる如何なる損害に対しても、その原因の如何を問わず、開発者は一切の責任を負いません。

Mersenne Twister について

Mersenne Twister（メルセンヌ・ツイスタ）は、松本眞先生（東京大学大学院数理科学研究科教授）と西村拓士先生（山形大学理学部数理科学科准教授）が開発なさった優れた擬似乱数生成アルゴリズムです。著作権はこの両名に属します。
Mersenne Twister はBSDライセンスに基づき無償かつ無保証にて頒布されるものです。

2011年7月22日　鵺