より良い精度

Question

a1は、第一項とすると、rは定数各用語は、次の用語及びnは、用語の数で得るために乗算されることを、等比数列であること：製品ai = a1*r**(i-1)、pnn用語のsnの合計、n用語の合計。

私はこれを計算する計算式を持っていますが、Fortran 95（Plato2）は必要な精度を認めていません。 （例：結果として-1.234E + 00567890を取得できません）。

この "巨大な"数値を使用する倍精度はどのように「十分に」使用できますか？

Answer 1

このような高い数値（例：-1.234E + 00567890）は、Fortran標準が提供する任意の組み込み数値型に対しては大きすぎます。物理的および工学的アプリケーションで使用される数値よりも大きくなります。例えば、私のgfortranは、これらのタイプをサポートしています。

huge(1.0_real32) 3.40282347E+38 
huge(1.0_real64) 1.7976931348623157E+308 
huge(1.0_real128) 1.18973149535723176508575932662800702E+4932 

限り、私ははるかに大きい固有の浮動小数点型で利用可能なFortranのコンパイラはありません知っています。

特殊な目的のために、あなた専用のライブラリが必要です。このサイトはソフトウェアに関する推奨事項ではありませんので、特にお勧めしません。それらのいくつかのリストをhttp://crd-legacy.lbl.gov/~dhbailey/mpdist/で見てください。もちろん、もっと多くのものがあります（GNU科学図書館には確かなことがあります）。