疑似乱数２

[1 デフォルトの名無しさん 2007/10/17(水) 22:34:59]

擬似乱数発生器について語ろうか。その２

前スレ
擬似乱数
http://pc11.2ch.net/test/read.cgi/tech/1146071975/

関連スレ
【危険】とんでもプログラム告発スレッド【悪質】
http://pc11.2ch.net/test/read.cgi/tech/1191860116/


SIMD-oriented Fast Mersenne Twister (SFMT):
http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/SFMT/index.html

[2 デフォルトの名無しさん 2007/10/17(水) 22:53:56]

で、いつになったらSFMTはboostに組み込まれるの？

[3 デフォルトの名無しさん 2007/10/17(水) 23:16:58]

俺用メモ

18* 名前：デフォルトの名無しさん [] 投稿日：2006/04/28(金) 23:53:29
軽さで言えばXorShiftとか。

unsigned long xor128(){
static unsigned long x=123456789,y=362436069,z=521288629,w=88675123;
unsigned long t;
t=(x^(x<<11));x=y;y=z;z=w; return( w=(w^(w>>19))^(t^(t>>8)) );
}

[4 デフォルトの名無しさん 2007/10/17(水) 23:35:36]

>>3
せっかくだからURLも貼っとこう
http://web.comlab.ox.ac.uk/oucl/work/richard.brent/pub/pub218.html

[5 デフォルトの名無しさん 2007/10/17(水) 23:38:04]

この手の乱数マジックナンバーでよくでてくる
123456789
なにふざけてるのかと思ってたら、
0.1234567891011121314… って超越数なんだな。
つまり123456789は、その超越数を10億倍して9桁の精度で切り落とした、
そこそこ質のよい数なのであった。
目からウロコ

[6 ヽ・´∀｀・,,）っ━━━━━━┓ 2007/10/18(木) 00:42:09]

マジックナンバーに0xDEADBEEFってよく見る

[7 デフォルトの名無しさん 2007/10/18(木) 01:51:23]

SFMTはCPUの機能に最適化させてるメルセンヌツイスタですよね？
boostはテンプレートライブラリだからそれはいつまでたっても組み込まれないと自分は思います…。

[8 デフォルトの名無しさん 2007/10/18(木) 08:32:57]

>>7意味不明

[9 7 2007/10/18(木) 10:16:21]

C++のテンプレートライブラリとして実装するものなのかなと思っただけです。

[10 デフォルトの名無しさん 2007/10/18(木) 12:38:11]

独自の乱数発生器を組み込めるのがboost::randomのよいところ。
手直しは必要だが。

っていうか、boostがテンプレートライブラリだなんて誰が言ってるの？
STLの事じゃなくて？

[11 7=9 2007/10/18(木) 14:19:16]

実際違いました…すみません。
boostの一部がSTLに移植されるとか書いてあったの見てたから勘違いかな…。
regexとかはテンプレートライブラリではなかったです。

あと、SFMTの情報斜め読みしてたから勘違いしてしまいました。申し訳ありません。

[12 デフォルトの名無しさん 2007/10/18(木) 19:08:52]

RandomはTR1に入った、つまり標準ライブラリ入り内定ではある。

Regexは、テンプレートライブラリと呼ぶかどうかはともかく、
basic_regexなどテンプレートはよく使っている。

[13 ヽ・´∀｀・,,）っ━━━━━━┓ 2007/10/19(金) 00:54:05]

regexがプリコンパイル方式なのは型が最初から決まってるからじゃん。

[14 デフォルトの名無しさん 2007/10/19(金) 09:58:49]

Xorshift RNGs

unsigned long xor128(){
static unsigned long x=123456789,y=362436069,z=521288629,w=88675123;
unsigned long t;
t=(x^(x<<11));x=y;y=z;z=w; return( w=(w^(w>>19))^(t^(t>>8)) );
}

ttp://lucille.atso-net.jp/wiki/index.php?%CD%F0%BF%F4