次のプログラムの出力はどうあるべきですか?
#include <stdio.h>
#include <math.h>
#include <float.h>
int main() {
int exp;
float mant = frexp(FLT_MAX, &exp);
printf("frexp(%a) = {%f, %d}\n", FLT_MAX, mant, exp);
return 0;
}
私の amd64 Linux システムでは、glibc を使って次のように表示されます。
frexp(0x1.fffffep+127) = {1.000000, 128}
IEEE 754 ウィキペディアの記事から、「最大正規化数」の指数は 127 であることがわかります。また、指数が 128 の値は ±Infinity と NaN だけであることがわかります。
frexp のマニュアル ページから、 frexp() は [0.5, 1.0) の範囲 (つまり、1.0 を除く) の値を返す必要があることを理解しています。
これらの情報に基づいて、返された仮数と指数の両方が正しくないようです。
frexp() の機能と FLT_MAX (= (2 - 2^-23) * 2^127) の値を知ると、(2 - 2^- 23) は 2 に非常に近い。
では、frexp(FLT_MAX, ...) は何を返すべきでしょうか?