Mastodawn

tetuさん、お世話になっています。
ついに、ついに、ついに、………

Julia と gcc の死闘、ついに終止符が打たれるのか？！？！？！

40億回のモンテカルロループにかかった時間

Julia帝国軍 → 4.2秒

gcc dSFMT OpenMP FMA 連合軍 → 1.9秒

ついに、gcc連合軍の圧倒的勝利！

しかし、たくさんの血が流された……………

OMP_NUM_THREADS の値によってどれだけ速さが変わるかに興味がある人は添付画像を見て下さい。

CPU: Intel(R) Core(TM) i7-4770HQ CPU @ 2.20GHz

ちなみにマクロ定義を除いたJulia側のコードは次の一行だけです。

findpi_macrosum(n::Int) = (4.0/n) * (＠sum i 1:n ifelse(rand()^2 + rand()^2 <= 1.0, 1, 0))

この1行だけで円周率をモンテカルロで計算できる。