自分にわかりやすく書いた備忘録です。出力はじめにこの記事は前回のプログラムを改良したものについて扱います。ソースコードの意味を知りたい方は、前回の記事を先にご覧ください。doubleの精度が限界でした。そこで、精度を理論上、メモリの限界まで上げることができるライブラリである、"GMP"というものを使用しました。これにより、8GBメモリ環境では余裕で前回の精度を超えることができました。library"の略で、多倍長の計算を簡単に実装できるライブラリとなっています。ちなみに、これはUnix系のOSでのみ動作します。Windowsで動かすには、Cygwinを入れるか、GMPの代わりにMPIRというものを入れる必要があります。幸運にも、私のPCは、LinuxとWindowsをデュアルブートしていたので、Linux環境で計算しました。Windows環境をご利用の方は、Cygwinを入れるのをお勧めします。3.40GHzRAM:8.00GBちなみに、これはサブPCです。GMPに移植この作業が本当に大変でした...。初期化する関数を使わなければならない。四則演算はすべて専用の関数で行う。メモリを最後に開放しなければならない。計算ごとに中間状態を保存するtemp変数を使わなければならない。最後、どうしてもセグフォになると思って、ググってみると、変数の初期化、開放関数を複数の変数に対して同時に使う時は、最後の引数にNULLを入れないといけないとか。自分の調べた限りでは、英語のサイトしか出てこなかったので、余計厄介でした。リファレンスを読めってことですね。コードを移植した結果こうなりました。0;}やっていることは、前回の記事とほぼ同じです。p2を経由しないように変更したり、出力先をテキストファイルにしたり。GMPに興味のある方は、こんなサイトなどを参考にしてみてください。スポンサーリンク出力0.3141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117067982148086513282306647093844609550582231725359408128481117450284102701938521105559644622948954930381964428810975665933446128475648233786783165271201909145648566923460348610454326648213393607260249141273724587006606315588174881520920962829254091715364367892590360011330530548820466521384146951941511609433057270365759591953092186117381932611793105118548074462379962749567351885752724891227938183011949129833673362440656643086021394946395224737190702179860943702770539217176293176752384674818467669405132000568127145263560827785771342757789609173637178721468440901224953430146549585371050792279689258923542019956112129021960864034418159813629774771309960518707211349999998372978049951059731732816096318595024459455346908302642522308253344685035261931188171010003137838752886587533208381420617177669147303598253490428755468731159562863882353787593751957781857780532171226806613001927876611195909216420199e1これが生のデータです。最後についているe1は、10の1乗という意味なので、正しく3.14...になっていることがわかります。天下のWikiから、正しいπのデータを引用してきます。...比較してみると、正しく1000桁(正確には998桁)計算できていることがわかります。プログラムのパラメータを書き換えれば、メモリの限り精度をいくらでも上げられます。ちなみに、1000桁程度では、enterを叩いた瞬間に計算が終わるほどの速さです。機会があれば、高みを目指していきたいものです。ご覧いただき、ありがとうございました。前回はライプニッツの公式編を書きました。それでは収束が遅く、精度がなかなか出にくいということがわかりました。円周率といえば、理系の人間ならだれでも(?)憧れを抱く無理数ですね。公開鍵暗号の1つである「RSA暗号」について理解を深めるため、実装してみました。円周率は理系が崇拝する数値で、さまざまな求め方が数学者たちによって研究されてきました。コメントをキャンセルメールアドレスが公開されることはありません。Tanidaizのプログラミング備忘録・日記.。