逆数を求めるプログラム

巨大な桁同士の高精度な除算をするには、普通に除算をするよりも、逆数を求め、それを掛けた方が遥かに速く計算できる。つまり、Ａ÷ＢはＡ×(１／Ｂ)として計算する。ここでＢの逆数、つまり、１／Ｂはニュートン法を用いると乗算のみで計算が可能。

このプログラムは、巨大桁ＡとＢを乱数で求めた上で、Ｂの逆数を求め、Ａに掛けることによりほぼ１になることを確かめる除算検証プログラムです。
また、ＦＦＴ部分は大浦さん公開の汎用FFTパッケージを使用します。
尚、intは32ビット、long doubleは80ビットタイプを想定してます。

ちなみに、このプログラムを円周率計算ソフトpt278wに使用しました。

参考にしたサイト
逆数と平方根を求める高次収束アルゴリズム

関連
円周率計算ソフト Windows版
FFTを用いて多倍長乗算をしてみた
円周率１００万桁を計算する

人気ベスト２０(ページ内)
トップページに戻る
作成日2007/2/8,最終更新日2008/04/11