プログラムが止まる。 *** glibc detected *** ./simulation1-4.exe: double free or corruption (!prev): 0x0000000000518350 ***こんなメッセージを残して落ちてしまう。 どうやら、ポインタとかのメモリ解放をしないと行けないらしい。 調べんと。参考 追…

メルセンヌツイスターを使ったときのコンパイル SIMD-oriented Fast Mersenne Twister (SFMT)の場合。 Dオプションで、臨時にマクロを定義している。 gcc -DMEXP=607 -o sample1 SFMT/SFMT.c sample1.c メルセンヌツイスターのコンパイル法 gccのオプション…

Cのプログラム、できるかな実験 Testという実行プログラムが一つ。 インプットファイルや(rの読み込み)アウトプットファイル(aの上書き)が重複していても、大丈夫。具体的には、以下の２つを同時に実行させても大丈夫、ということ。 ./Test [Input1(r)] [Inp…

一様乱数法による生成で、genrand_res53()が[0,1)の一様乱数に対応する。 int Poisson(double lambda){ int k=0; lambda=exp(lambda)*(genrand_res53()); while (lambda>1) { lambda*=(genrand_res53()); k++; } return k; } しかし、lambdaの値が大きくなる…