並列化の試行やJITの効果テストなどのために，Forをたくさん回す計算例が必要になる。 そのような例としてモンテカルロ法を置いておく。 ToDo 並列計算のテスト用としては，モンテカルロ法のような単純に並列化される(データ共有が無い)ものだけでなく，適時…

Rによるモンテカルロ法入門のメモ，その１。 イントロ 統計的推論では最適化と積分の計算，特に実用上は数値計算が必要になる。 数値積分の手法としてモンテカルロ積分がある。シンプソン法などの区分法は多次元の場合に区分を細かくしようとすると計算量が…

timeモジュール(time.perf_counter) 使いどころは特定のコードブロックに対する計測。精度はマイクロ秒程度。精度を求める場合はこれ。(逆にcProfile/line_profileなどのプロファイリングは基本的にクロック程度の精度,1/100秒程度しかなく，ボトルネックの…