アカウント名:
パスワード:
高橋先生のWebに詳細データが有りました。
http://www.hpcs.is.tsukuba.ac.jp/~daisuke/pi-j.html [tsukuba.ac.jp]
今までの記録である2002年の1兆2411億桁の記録は、マチンの公式 [wikipedia.org]で代表される逆正接関数を使った円周率の公式を、後保範先生らによって考案されたDRM法(Binary Splitting法の一種)で求めたものでしたが、今回はそれ以前の記録で使われていたガウス・ルジャンドルアルゴリズム [wikipedia.org]に戻ったようですね。
FFTによる乗算の畳み込みは、この桁数になると倍精度浮動小数点で精度保証をするのはかなりきつい(一個の変数で表せる桁数が極端に減る)と思うのですが、その辺は力技なのか、既知の手法があるのか、気になりますね。
もちろん乗算用の整数演算FFTがあります。libgmp4 辺りを見るといいんじゃないでしょうか。確か各種論文へのリファレンスもあったと思います。
より多くのコメントがこの議論にあるかもしれませんが、JavaScriptが有効ではない環境を使用している場合、クラシックなコメントシステム(D1)に設定を変更する必要があります。
Stableって古いって意味だっけ? -- Debian初級
アルゴリズム (スコア:4, 参考になる)
高橋先生のWebに詳細データが有りました。
http://www.hpcs.is.tsukuba.ac.jp/~daisuke/pi-j.html [tsukuba.ac.jp]
今までの記録である2002年の1兆2411億桁の記録は、マチンの公式 [wikipedia.org]で代表される逆正接関数を使った円周率の公式を、後保範先生らによって考案されたDRM法(Binary Splitting法の一種)で求めたものでしたが、今回はそれ以前の記録で使われていたガウス・ルジャンドルアルゴリズム [wikipedia.org]に戻ったようですね。
Re:アルゴリズム (スコア:1)
この桁数になると倍精度浮動小数点で精度保証をするのはかなりきつい
(一個の変数で表せる桁数が極端に減る)と思うのですが、
その辺は力技なのか、既知の手法があるのか、気になりますね。
Re:アルゴリズム (スコア:3, 興味深い)
もちろん乗算用の整数演算FFTがあります。
libgmp4 辺りを見るといいんじゃないでしょうか。
確か各種論文へのリファレンスもあったと思います。