円周率計算でスパコンの記録を20万のパソコンが破る!

| コメント(13) | トラックバック(0)
なんとなくYahooニュースを見てたらこんな記事が載っていた

http://headlines.yahoo.co.jp/hl?a=20100112-00000439-yom-sci

パリ在住のフランス人デジタルテレビ関連技術者が、円周率を2兆6999億9999万けたまで計算し、世界記録の更新を宣言した。

 筑波大学計算科学研究センターのスーパーコンピューターが昨年樹立した記録を塗り替えたことになる。

 ホームページで新記録を公表したのはファブリス・ベラールさん。筑波大のスパコンが昨年、73時間36分かけて樹立した世界記録は2兆5769億8037万けただったが、ベラールさんはこれを約1230億けた上回ったことになる。

 また、ベラールさんが計算に使ったのはスパコンではなく、価格2000ユーロ(約26万円)以下の普通のデスクトップ型パソコン。2進法による計算に103日、検算に13日をかけたという。


ベラールさん すげー、って一瞬思いかけたけどよく見ると
筑波大は73時間、つまり3日ほどで作った記録、
対してベラールさんは100日以上かけてるし・・・

スパコンに対して力業で対抗したわけですねw

もしかしたら筑波大はなにくそ~って今週7日間費やして ベラールさんなんか比較にならないような結果を出してくるかもしれませんね(笑)


ニュース書いた人の作為が・・・(笑)

ただここで終わっては普通に終わってしまうので WEB系エンジニアとしては実際に円周率の計算ってどうやるの?って調べて実際にやってみました。
円周率の定義って (円周率)=(円の外周)÷(円の直径)と小学生で習いました。

そうすると円の外周のとんでもなく正確な長さがわからないと何兆桁の計算はできない??
とりあえずWikipediaで調べてみると

http://ja.wikipedia.org/wiki/%E5%86%86%E5%91%A8%E7%8E%87#.E6.8C.87.E6.95.B0.E9.96.A2.E6.95.B0.E3.83.BB.E4.B8.89.E8.A7.92.E9.96.A2.E6.95.B0.E3.81.AB.E3.82.88.E3.82.8B.E5.AE.9A.E7.BE.A9

円周率は、元々、平面幾何学で定義された数であるが、平面幾何学に留まらず数学のいろいろな分野でその姿をあらわす。円周率 π の導出方法は多いが、どの定義によっても同じ値が得られるので、その時に応じて使いやすい定義を用いればよい。

なんと円の外周とか円の面積とか平面幾何の世界以外でも三角関数や指数関数を使った方法、積分を使った方法でも計算できるらしい。
そのうちの一つ積分を使って求める方法の展開した結果 べき乗と階数を使った式に変えられるみたいなのでせっかくその式に従ってPHPでプログラムを組んでみた。
(式はWikipedia参照のこと)

<?php
print getPai(20);

function getPai($n){
$ret = kai($n*2)
/ (pow(2,4*$n+1) * (pow(kai($n),2)) * ($n*2+1)) * 6;

if($n == 0){
return($ret);
} else {
return($ret + getPai($n-1));
}
}

function kai($n){
if($n == 0){
return(1);
} else {
return($n * kai($n-1));
}
}
?>


この出力結果が
3.1415926535898

正しい円周率が
3.1415926535897932384
ってことなので なんとこんな簡単なプログラムで円周率を12桁まで計算できてしまった!!

ここで引数の数字を大きい数字にしていけばより正確な値が出るはずだけど
PHPの取り扱える最大桁数にひっかかったり 小数の扱える最大桁数にひっかかるので
これ以上の計算をするにはより多くの桁数を扱える独自のクラスを作る必要があるみたい。
このあたりの工夫は遊びでやるにはちょっと大変なので今回はパス。

ということで5分で適当に作ったプログラムで円周率12桁まで計算!
ギリシャ時代とかならたぶん世界記録(笑)

トラックバック(0)

トラックバックURL: http://blog.romaji.net/mt4/mt_cgi/mt-tb.cgi/70

コメント(13)

Wanted to state your web site is kinda awesome. It's my job to want to hear a new challenge about this because We have the same blog inside my Country during this subject so this help´s us a lot. I have done specific searches to the matter and discovered a number of blogs but nothing can beat this.Thank you for writing a whole lot in the blog.. Greets, Fran

Really good write-up,reasonable advice. With thanks =)

I accept you ws assurance may able-bodied yield amusement there.I realize that every of them capacity ability be much wider aces a bulletin lath article, but I anticipation they could anniversary be accordant to the column and remarks.Thanks a lot

Hey! Among the best to grant a tremendous thumbs up to the good info you might have here during this post. I'll likely be coming again to the blog with regard to added soon.

Hi! You will find certainly a great deal of details that adheres to that to consider. That is a great denote bring up. I provide thoughts above as general inspiration but clearly you will discover questions including the one you talk about in which the most critical thing will likely be employed in honest good faith. I don?t know if recommendations have emerged around items like that, however i guarantee that a job is clearly known as an affordable game. Both boys and girls have the impact of simply a moment's pleasure, for the rest of their lives.

Good written story. It'll be employed to everyone who usess it, including me. Continue the best work - i most certainly will definitely continue reading posts.

Awesome website man, I wish I'm able to make such like by incorporating content like this lmao.

What a lovely day for a 369165! SCK was here

I simply needed to say thanks yet again. I do not know what I would've tried without the opinions documented by you directly on my theme. Completely was the difficult scenario in my position, nevertheless taking a look at this skilled way you treated it took me to jump over gladness. I am happy for your work and as well , expect you are aware of an amazing job you were putting in teaching many people using your web blog. Most probably you haven't encountered all of us.

Today, while I was at work, my cousin stole my iphone and tested to see if it can survive a twenty five foot drop, just so she can be a youtube sensation. My iPad is now destroyed and she has 83 views. I know this is entirely off topic but I had to share it with someone!

Hallo, always great to discover others through the hole world around my searching, I really appreciate some time it must took to create this awesome article. Cheers

Thanks so much for giving everyone an extraordinarily superb opportunity to check tips from this blog. It's usually very terrific plus jam-packed with amusement for me and my office acquaintances to search your web site at the least three times in a week to find out the newest secrets you have got. Of course, I am just usually pleased concerning the astounding hints served by you. Certain 1 areas in this post are undeniably the most beneficial we have ever had.

I usually hold my opinion to myself and don’t leave any blog comments. Couldn’t help myself after reading your posting though. Terrific document. Thanks!

コメントする

2011年9月

        1 2 3
4 5 6 7 8 9 10
11 12 13 14 15 16 17
18 19 20 21 22 23 24
25 26 27 28 29 30