筑波大の高橋さんらが円周率を2兆5769億8037万桁まで計算→アルゴリズムをkeisan.casio.jpにUP!
円周率の計算といえば東大の金田さんグループだと思っていましたが、今回の記録更新は筑波大学の高橋さんグループでした。
http://www.yomiuri.co.jp/science/news/20090817-OYT1T00638.htm
でもやはり金田さんと一緒に仕事をされてて以前の計算にも関わっていた方(金田研出身)で、共書も多いみたいですよ。
ご自身のwebサイトによると、アルゴリズムは主計算にガウス・ルジャンドルアルゴリズム、検証計算にボールウェインの 4 次の収束アルゴリズムを用いられたということです。
http://www.hpcs.is.tsukuba.ac.jp/~daisuke/pi-j.html
ガウス・ルジャンドルの説明はこちら。
http://en.wikipedia.org/wiki/Gauss-Legendre_algorithm
いわゆる算術幾何平均法(AGM)。
a0=1
b0=1/√2
t0=1/4
p0=1
を初期値として、
an+1 = (an + bn) / 2
bn+1 = √(anbn)
tn+1 = tn - pn * (an+1 - an)^2
pn+1 = 2*pn
として、
π≒(an + bn)^2 / (4*tn)
で計算したもの。これだけで求まるのは不思議ですね。
早速高精度計算サイトkeisan.casio.jpにUPしてみた。(これで101個目)
リンクはこちら。
最初のほうの計算はこんな感じでかなり収束がはやい。
0:2.9142135623730950488016887242096980785696718753769
1:3.1405792505221682483113312689758233117734402375129
2:3.1415926462135422821493444319826957743144372233456
3:3.1415926535897932382795127748018639743812255048354
4:3.1415926535897932384626433832795028841971146782836
5:3.1415926535897932384626433832795028841971693993751
ボールウェインの方はこちら。
http://en.wikipedia.org/wiki/Borwein%27s_algorithm
両方とも反復計算アルゴリズムなんですね。私の知っているコンピュータでπの計算法といえば、マチンの式とかその変形、ラマヌジャンの式とかなんですが、だいぶ違ってますね。並列コンピュータで解きやすいのかな。
※追記 この文献などの内容が使われたのかな?
http://ci.nii.ac.jp/naid/110002932029/
"分散メモリ型並列計算機による円周率の高精度計算"
« 大久保のイズミヤの看板は普通 | トップページ | 豊中稲荷神社に行ってきた。 »
「パソコン・インターネット」カテゴリの記事
- 高周波回路シミュレータQucsStudioを使ってみる(その3)Mixed Mode S parameterを計算(2018.10.12)
- 高周波回路シミュレータQucsStudioを使ってみる(その2)SパラメータのTouchStoneフォーマットで出力するには?(2018.10.11)
- 高周波回路シミュレータQucsStudioを使ってみる(その1)まずは何をさておきμの文字化けだけには注意。(2018.10.10)
- 円の弧長,弦長,矢高,半径のどれか2つを与えて残りを計算(カシオの高精度計算サイト自作式)で180°以上、複数解に対応。(2018.10.09)
- macbook proをmacOS Mojaveにアップデート。せっかくなんでダークモードにしてみる。(2018.09.26)
「学問・資格」カテゴリの記事
- 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2021年2月16日 IEEE Microwave Magazineの特集はオールデジタルのRFID、Microwave JournalはEバンド ミリ波通信に衛星や気球を使う話、アメリカの半導体企業がバイデンに投資を迫る、(2021.02.17)
- カオスを生じる電気回路、Chua’s circuitをLTspiceで回路シミュレーションしてみる。(2021.02.19)
- Labyrinth Chaos(迷宮カオス)を生むThomas-Rössler方程式のパラメータbを色々変えて、Python+Scipyでルンゲクッタ8次のDOP853(Dormand&Prince)を使って計算してGIFアニメ(2021.02.16)
- フィッツヒュー・南雲 (FitzHugh-Nagumo) 方程式をPython+Scipyでルンゲクッタ8次のDOP853(Dormand Prince)で計算。(2021.02.23)
- 「水晶振動子の等価回路計算」をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP! インピーダンスの大きさと位相がグラフ化できる。(2021.02.12)
「日記・コラム・つぶやき」カテゴリの記事
- ExcelのLET関数+SEQUENCE関数で数値計算シリーズ(その4)ワンライナーでシンプソンの積分公式でクロソイド曲線を計算する。(2021.02.25)
- ExcelのLET関数+SEQUENCE関数で数値計算シリーズ(その3?)ワンライナー(1行というか 1セル)で数値積分のシンプソンの公式を計算する。(2021.02.24)
- 新型コロナウイルス、中国、日本、韓国、アメリカ、ドイツ、フランス、イギリスでの感染者数を指数関数&ロジスティック関数&Log-Logプロットでべき関数フィッティングした(2/21更新)さすがに日本も増加率は鈍化してきた。しかし中国が2/7から更新データがない、、、(2021.02.22)
- 点Pのあおり運転をGeogebraで3次元GIFアニメにしてみた。(虚構新聞の記事:辺上で急接近や急減速 あおり運転の疑いで点Pを逮捕を見て)(2021.02.18)
- 新型コロナウイルス、中国、日本、韓国、アメリカ、ドイツ、フランス、イギリスでの感染者数を指数関数&ロジスティック関数&Log-Logプロットでべき関数フィッティングした(2/14更新)増加率はさすがに減少したものの、まだ日本が一番増加率が高い。(2021.02.15)
コメント