[解決済み] CとC++でDBL_EPSILONを(ポータブルに)取得する方法

質問

Linux (AS 3)上でGCC 3.4を使っているのですが、以下のことを理解しようとしています。 DBL_EPSILON あるいは、少なくともまともな近似値です。 どのようにすればプログラム的にそれを得ることができますか？

どのように解決するのですか？

それは "float.h" にあるはずです。これはCとC++の標準の一部で、移植性があります（C++では非推奨ですが <cfloat> またはsbiさんの回答でquot;guaranteed"前方互換性があります)。

もし、それがなければ、あなたの倍数はIEEE64ビットなので、他の誰かのfloat.hから値を盗めばいいのです。私が見つけた最初のものはこちらです。

http://opensource.apple.com/source/gcc/gcc-937.2/float.h

#define DBL_EPSILON 2.2204460492503131e-16

この値は私にはほぼ正しいように見えますが、もしコンパイラで確認したい場合は (1.0 + DBL_EPSILON) != 1.0 && (1.0 + DBL_EPSILON/2) == 1.0

Edit: "programmatically"の意味がよくわからないのですが。これは標準的な定数で、計算することは想定されていません。しかし、このようなことは可能だと思います。DBL_EPSILON は、1.0 の表現の最後の精度のビットで 1 を表す 0.5 の何乗かになるように、再び IEEE 表現か何かを仮定しています。

double getDblEpsilon(void) {
    double d = 1;
    while (1.0 + d/2 != 1.0) {
        d = d/2;
    }
    return d;
}

コンパイラの設定によっては、中間結果の精度が double の結果が小さくなってしまいます。 d よりも DBL_EPSILON . コンパイラのマニュアルを確認するか、強制的に 1.0 + d/2 を格納し、リロードして実際の double オブジェクトと比較する前に 1.0 . 大雑把に言うと、PCではコンパイラがx86 FPU命令（高精度）を使うか、新しいx64浮動小数点演算（倍精度）を使うかによって変わります。