昔の記事を自分から見に行くということはほとんどないが,ときどき誰かが読んでくれたという記録に目がとまり,一体どんなことを書いたのかなと確認することがある。
「ハーディ・ラマヌジャン数」もそのようなものだった。
それは,R で書いたプログラムであったが,今回 Julia で書いて,速度を比較してみた。
Julia のプログラムは
function ramanujan4()
m = floor(Int, 9999^(1/3))
x = (1:m) .^ 3
z = freqtable([x[i] + x[j] for i=1:m, j=1:m if i < j])
z[z .>= 2]
end
Benchmark によれば,実行時間の中央値は 56 μs であった。R は 29 ms だったので,Julia は R より,500 倍速いということだ。
※コメント投稿者のブログIDはブログ作成者のみに通知されます