float.h　のテスト結果がおかしい

Question

float.h　　をテストするため、次のようなプログラムを作りました。
　float.h　に書かれている値が書き出されると思ったのですが、少し値が違います。
　FLT_MAX と DBL_MAX の値が違っています。本当にこのヘッダファイルを読んでいるかと思い、
#define FLT_MAX 3.4e+3f　　　　と変更すれば、
Max= 3.400000000000000000E+03　と出てきます。

　一致しない原因は何でしょう。
　コンパイラはルネサスのＨＥＷで、ＣＰＵはＨ８／３０５２です。

＊＊＊　float.h　＊＊＊

#define FLT_MAX 3.4028235677973364e+38f
#define FLT_MIN 1.175494351e-38F

#define DBL_MAX 1.7976931348623158e+308
#define DBL_MIN 2.2250738585072014e-308

＊＊＊　プログラム　＊＊＊

sprintf(str,"     float :%ldBytes  Max= %.18E   Min= %.18E
",sizeof(float), FLT_MAX, FLT_MIN) ;
Put_str(str);
sprintf(str,"     double:%ldBytes  Max= %.18E  Min= %.18E
",sizeof(double), DBL_MAX, DBL_MIN) ;
Put_str(str);


＊＊＊　結果　＊＊＊

float :4Bytes  Max= 3.402823466385288600E+38
　　　　　　　　Min= 1.175494350822287500E-38

double:8Bytes  Max= 1.797693134862315700E+308
　　　　　　　　　Min= 2.225073858507201400E-308

＊＊＊　まとめて整理　＊＊＊
　元の原稿は比較し易いように縦に数字が並ぶように書いているのですが、
ここに書き込むとずれてしまいます。

#define FLT_MAX 3.4028235677973364e+38f
　　　　　　　 Max= 3.402823466385288600E+38

　FLT_MIN に比べ有効桁数が多いのが気になる。
　同じ有効桁数で区切っても ...567 と ...466　とでは誤差が大きい。

　　　FLT_MAX 3.402823567 7973364e+38f
　　　　　 Max= 3.402823466 385288600E+38
　　　 FLT_MIN 1.175494351e-38F



#define FLT_MIN 1.175494351e-38F
　　　　　　　 Min= 1.175494350822287500E-38

　四捨五入なら分かる。


#define DBL_MAX 1.7976931348623158e+308
　　　　　　　 Max= 1.797693134862315700E+308
　最後が８と７で違う。しかし、DBL_MAX は８バイトで表せる値に対して四捨五入したために７が８になったと考えれば納得出来る。


#define DBL_MIN 2.2250738585072014e-308
　　　　　　　 Min= 2.225073858507201400E-308
　ピッタリ合っている


　宜しくお願いします。

jacta · Accepted Answer

> 誤差が拡大するのは分かるのですが、次の部分の誤差が大き過ぎると思います。

FLT_DIGを調べれば分かることですが、単精度浮動小数点数の精度は10進法で6桁程度しかありません。

chie65535 · Answer

＞誤差については分かるのですが、次の部分の誤差が大き過ぎると思っています。

単精度の仮数部は、省かれた１ビットを含め、２５ビットしかありません。

つまり「２の２５乗」の精度しかありません。

「２の２５乗」は「33554432」ですから「１０進数で７桁」だけが有効です（２進数の0～33554432の範囲で表わせる１０進数は0～9999999までの７桁）

「７桁だけ有効」と言うことは「整数部１桁＋小数部６桁」しか信じられないので、小数部６桁で切ると
#define FLT_MAX 3.402823 5677973364e+38f　ＨＥＷ
実験結果　 Max= 3.402823 466385288600E+38
となり、切った小数７桁以下は「無効なので無視しなければならない桁」です。

小数７桁以下を無視すると、どちらも「同じ数値」になります。

言い換えると「小数点以下７桁目が４でも５でも、２進化して単精度の浮動小数点数にコンバートすると、同一のバイナリ値が出来あがる」って事です（小数点以下７桁目が３とか２とかだと、微妙に小さい値になってしまうでしょうから駄目です。また、６とか７とかだとオーバーフローの危険があるので駄目です）

試しに
float f1;
float f2;
f1=3.4028235677973364e+38f;
f2=3.402823466385288600E+38f;
sprintf(str,"%s
",f1 == f2 ? "EQU" : "Not EQU");
Put_str(str);
と言うプログラムを書いてみましょう。

実行すると
EQU
になる筈です。

この結果により「１０進表記が異なっていても、異なるのは有効桁の外（つまり無効桁）なので、内部のバイナリ値が一緒になり、f1とf2が等しくなる」のが判るでしょう。

chie65535 · Answer

定義は合ってますが、表示方法が間違っています。

sprintfなどの「printf系関数」は「単精度実数を引数に与えると、倍精度実数に変換する」のです。

この変換の際に、４バイトから８バイトに拡張する為、大きく誤差が出ます。

しかも、１０進数表記（「0.123E+30」など、人間が理解出来る表記）にする際に、丸め誤差など、更に誤差が蓄積します。

簡単に言うと「貴方が目で見ている表示は、いろいろなキャストや変換が行われた後の結果で、元の数値とは似ても似つかないもの」なのです。

少なくとも「float⇒double」のキャストによって「全然違う値」になります。

実験するのであれば「FLT_MAX、FLT_MINを、floatの変数に代入し、変数の中身を16進ダンプする」のが正解です。

そして、16進ダンプした結果を見て、バイナリで「指数部が最大値になってて、仮数部が全部１で埋まってる」のを確認しましょう。

因みに「3.4028235677973364e+38fと書かれたソース上の文字列を、floatのバイナリ値に変換する際」に「既に誤差を含んでいる」のをお忘れ無く。

何故なら「１０進数での有限小数が、２進数では無限小数になる」からです（「0.1」は「１０進数では有限小数」ですが「２進数では循環小数」になり、0.1を２進数で正確に表現する事は不可能なのです）

つまり「3.4028235677973364e+38fと言う１０進数の数値は、２進数では正確に表現不可能」なのです。

jacta · Answer

H8のHEWは長い間触っていないので詳しいことは忘れましたが、SHのものはfloatとdoubleを常に同じ（両方単精度または倍制度）にするオプションがあったように思います。
それと同じか、同じなくてもヘッダを使い回している可能性がありますので、有効桁数以上を記述していること自体はそれほど問題ありません。

ただ、floatが単精度で、doubleが倍精度の場合、FLT_MAXはどうしても単精度の有効桁数までしか表現できません。そして、可変個数引数としてprintfに渡す際にdoubleに変換されますので、その辺りで誤差が拡大することになります。

php504 · Answer

実際のfloatで表現できる最大値以上を定義しているからでしょう
そのfloat.hのFLT_MAXは明らかにおかしいですね
誰かが書き換えたんじゃないですか
ちなみにVisualC++では
#define FLT_MAX         3.402823466e+38F
でした

float.h のテスト結果がおかしい

> 誤差が拡大するのは分かるのですが、次の部分の誤差が大き過ぎると思います。

＞誤差については分かるのですが、次の部分の誤差が大き過ぎると思っています。

定義は合ってますが、表示方法が間違っています。

H8のHEWは長い間触っていないので詳しいことは忘れましたが、SHのものはfloatとdoubleを常に同じ（両方単精度または倍制度）にするオプションがあったように思います。

実際のfloatで表現できる最大値以上を定義しているからでしょう

似たような質問が見つかりました

関連するカテゴリからQ&Aを探す

デイリーランキングこのカテゴリの人気デイリーQ&Aランキング

マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング

float.h　のテスト結果がおかしい