円周率の計算をしてみる

computer
  • http://www.kurims.kyoto-u.ac.jp/~ooura/pi_fft-j.html

ここにある円周率計算プログラムで計算してみたら、意外にけっこうなかなか行けてしまった。昔はSuperPiとかいうプログラムで計算していて、Celeron400を使って3355万桁を1日半ぐらいで計算した覚えがある。それに比べれば隔世の感がある。やっぱり新しいパソコンは速いのか。まあ、円周率なんていったん出してしまえばあまり使い道はないんだけど、やっぱり計算機の速度を測定したりするのにはまだまだ楽しい。

Calculation of PI using FFT and AGM, ver. LG1.1.2-MP1.5.2af.memsave
nfft= 16777216
radix= 10000
error_margin= 0.349286
mem_alloc_size= 603979896
calculating 67108864 digits of PI...
AGM iteration,	time= 71,	chksum= ffffd2d2
precision= 48,	time= 160,	chksum= fffffe29
precision= 80,	time= 220,	chksum= ffffd0c6
precision= 176,	time= 280,	chksum= ffffc78d
precision= 352,	time= 341,	chksum= ffffc7b2
precision= 688,	time= 401,	chksum= ffffe009
precision= 1392,	time= 461,	chksum= fffff660
precision= 2784,	time= 521,	chksum= ffffdd09
precision= 5584,	time= 581,	chksum= fffff686
precision= 11168,	time= 641,	chksum= fffff787
precision= 22336,	time= 701,	chksum= ffffc041
precision= 44688,	time= 761,	chksum= ffffe90c
precision= 89408,	time= 821,	chksum= ffffc429
precision= 178816,	time= 882,	chksum= ffffd40a
precision= 357648,	time= 942,	chksum= ffffc2c5
precision= 715312,	time= 1002,	chksum= ffffd770
precision= 1430640,	time= 1062,	chksum= ffffc787
precision= 2861280,	time= 1122,	chksum= ffffc1b0
precision= 5722592,	time= 1182,	chksum= ffffc7b9
precision= 11445200,	time= 1242,	chksum= ffffdec8
precision= 22890416,	time= 1302,	chksum= ffffc152
precision= 45780848,	time= 1362,	chksum= ffffd1b7
precision= 91561728,	time= 1422,	chksum= ffffdae9
Total 1509 sec. (real time),	chksum= 8eb