単精度浮動小数点数

情報処理において、単精度浮動小数点数 は、コンピュータの数値表現の一種である。 本来の単精度浮動小数点数とは、倍精度浮動小数点数に対比して使われる用語であって倍精度浮動小数点数が無いシステムであれば単に浮動小数点数と呼ばれるべきものである。 そうして、単精度はシステムの基本となる精度であり、バイトマシンが普通になる前のワードマシンにおいて1語長分の記憶場所を占める浮動小数点数であるのが普通である。 昔のFORTRAN言語では、整数型と実数型は同じ語長を占めることを前提に言語の規格が制定されていたりもした。その記憶語1単位分のを占める実数型が単精度であり、記憶語2単位分を占める実数型は倍精度と呼ばれた。 Crayのベクトルマシンのような語長が64ビットである科学技術専用計算機では、たとえば単精度とは64ビットを意味した。

現在の浮動小数点形式の標準であるIEEE 754では、単精度は32ビット、倍精度は64ビットである。IEEE 754-2008 では単精度の形式は正式には binary32 と呼ばれている。「単」や「倍」という表現は32ビットを1ワードとする32ビットアーキテクチャを基にしている。

単精度浮動小数点数は、同じ幅の固定小数点数に比べてより広範囲な数値を表せるようになっているが、同時に精度を犠牲にしている。

ほとんどのプログラミング言語では標準の組み込み型として定義されている。C言語、C++、C#、Java、Haskellでは単精度のデータ型を float と呼ぶが、C/C++の規格ではIEEE 754に準拠することは要求されていない。MATLABでは single と呼ぶ。PascalおよびVisual Basic .NETでは Single と呼ぶ。.NET Frameworkの共通型システムではSystem.Single構造体型として定義されており、各.NET言語における組み込み型はこの型へのエイリアスである。date=2019-07。Pythonの float は倍精度であり、単精度のサポートはない。F#の単精度は float32 または single であり、倍精度は float または double である。

1. IEEE 754 での単精度浮動小数点数の形式: binary32
------------------------------------

IEEE 754 での binary32 の定義は次の通りである。

符号ビット: 1ビット
指数部の幅: 8ビット
仮数部の幅: 23ビット符号ビットは表現する数値の符号を示す。指数部は 8 ビットで、符号なし整数として見た場合は 0 から 255 の値をとり、0 の時はゼロか非正規化数であることを、1 から 254 の時は -127 のゲタ履き表現で -126 から 127 の指数を、255 の時は無限大または NaN であることを示す。

指数部が 0の場合を除き、仮数部で表現されるビットパターンのさらにひとつ上のケタに暗黙の 1 のビットがあるとみなす表現法により、通常の数の精度は、24ビット相当である。従って十進に換算したとき表現できる桁数は log10}} ≈ 7.225 桁となる。標準化されたレイアウトは以下のようになる。

符号ビットを sign、バイアスつきの指数部を e、23ビットの仮数部の各ビットを b-n}} とすると、32ビットの binary32 フォーマットで表される数値は となり、より正確に表現すると となる。

1.1. 指数部の符号化方式
--------------

単精度バイナリ浮動小数点数の指数部はオフセット表現を使って符号化されており、指数値がゼロのときのオフセット値は127である。

Emin = 01H−7FH = −126
Emax = FEH−7FH = 127
指数部バイアス = 7FH = 127指数部バイアスは、エクセスNとも言う。詳しくは符号付数値表現を参照されたい。真の指数値は、指数部の値から指数部バイアスを引いた値となる。

00H と FFH は予約された指数値である。

従って、全てのビットパターンが符号として意味がある。正の最小値は 2−149 ≈ 1.4 × 10−45 である。正の正規化数の最小値は 2−126 ≈ 1.18 × 10−38 である。表現可能な最大値は  × 2127 ≈ 3.4 × 1038 である。

1.2. 十進表現から binary32 フォーマットへの変換
-------------------------------

IEEE754 では実数値から対応する binary32 フォーマットへの変換を丸めも含めて厳密に定めている。

ここでは十進の実数値から IEEE 754 binary32 フォーマットへの変換は大まかに次のようになる。

実数値を 12.375 のように整数部と小数部で表すものとする。
整数部を二進法で表現したものに変換する。
小数部を以下に示す技法で変換する。
2つの変換結果を組み合わせ、さらに変換を加えて最終的な表現を得る。小数部の変換:

12.375 の小数部 0.375 を例とする。これを二進の小数に変換するため、小数部に2をかけて整数部をとり、さらに2をかけて整数部をとるというように繰り返していき、小数部がゼロになるか精度の上限である二進23桁になるまで続ける。

0.375 x 2 = 0.750 = 0 + 0.750 ⇒ b−1}} = 0 となるので、二進での小数部の一桁目はゼロとなる。そこでさらに 0.750 に2をかける。

0.750 x 2 = 1.500 = 1 + 0.500 ⇒ b−2}} = 1

0.500 x 2 = 1.000 = 1 + 0.000 ⇒ b−3}} = 1 となり、小数部が 0.000 となったのでここで変換は停止する。

これで 10}} は正確に二進で 2}} で表されることがわかった。十進で有限桁数の小数が必ず二進で有限桁で表せるわけではない。例えば、十進の 0.1 は二進では正確に表現できないため近似値を用いることになる。

したがって、10 = 10 + 10 = 2 + 2 = 2 となる。

IEEE 754 binary32 フォーマットでは、値を という形式で表現しなければならないので、1100.011 をシフトして整数部が一桁になるようにしなければならない。この例では3桁シフトするので となる。

最終的に となる。

この結果から次が得られる。

指数部は 3 となる。
仮数部は 100011 となる。したがって、IEEE 754 binary32 フォーマットで 12.375 を表すと 0-10000010-10001100000000000000000 = 41460000H となる。

注: 68.123 を IEEE 754 binary32 フォーマットで表す場合を考えてみよう。上述の技法を適用すると 42883EF9H というビット列が得られ、最後の4ビットは 1001 となる。しかし、IEEE 754 のデフォルトの丸め方式により最終的に 42883EFAH となり、最後の4ビットは 1010 となる。

例: 十進の1を変換すると となる。ここから次のことが導かれる。

指数は 0 である。
仮数は 0 である。従って実数 1 を IEEE 754 binary32 フォーマットで表すと 0-01111111-00000000000000000000000 = 3f800000H となる。

例: 0.25 を変換すると、 となる。ここから次のことが導かれる。

指数は -2 である。
仮数は 0 である。従って実数 0.25 を IEEE 754 binary32 フォーマットで表すと 0-01111101-00000000000000000000000 = 3e800000H となる。

例: 0.375 を変換すると、 となる。ここから次のことが導かれる。

指数は -2 である。
仮数は 1 である。従って実数 0.375 を IEEE 754 binary32 フォーマットで表すと 0-01111101-10000000000000000000000 = 3ec00000H となる。

1.3. 単精度浮動小数点数の例
----------------

浮動小数点数値のビット列を十六進法で表した例を以下に示す。これには符号、バイアスを加えた指数値、仮数値が含まれている。

3f80 0000 = 1 c000 0000 = -2

7f7f ffff ≒ 3.4028234 x 1038 

0000 0000 = 0 8000 0000 = -0

7f80 0000 = 正の無限大 ff80 0000 = 負の無限大

3eaa aaab ≒ 1/3

1/3 は倍精度とは異なり、切り上げられる。これは仮数部のビット数が偶数であるため、丸める桁位置以降が 1010... となってその桁位置で 1/2 より大きくなるためである。

1.4. 単精度バイナリ形式から十進への変換
----------------------

16進の値 41c80000 を例として変換を行う。これを二進で表すと次のようになる。

41c8 000016 = 0100 0001 1100 1000 0000 0000 0000 00002

これを、符号ビット、指数部、仮数部の3つに分割する。

符号ビット: 0 指数部: 1000 00112 = 8316 = 131 仮数部: 100 1000 0000 0000 0000 00002 = 48000016

ここで仮数部に暗黙の整数ビットを加える。

仮数: 1100 1000 0000 0000 0000 00002 = C8000016

指数部の値から127を引いて実際の指数値を得る。

指数部の値: 8316 = 131 本来の指数: 131 − 127 = 4

24ビットの仮数は最上位ビットが1に対応し、次の桁が0.5、その次が0.25というように、前の桁の1/2に対応している。

bit 23 = 1 bit 22 = 0.5 bit 21 = 0.25 bit 20 = 0.125 bit 19 = 0.0625 . . bit 0 = 0.00000011920928955078125

この例では、仮数で立っているビットは bit 23、bit 22、bit 19 の3ビットだけであり、上記の対応する値を加算することで十進での仮数が得られる。

本来の仮数: 1 + 0.5 + 0.0625 = 1.5625 = C80000/223

基数 2 を指数値でべき乗したものを仮数にかけると本来の値が得られる。

1.5625 × 24 = 25

従って

41c8 0000 = 25

となる。これは次の式と等価である。 ここで は符号ビット、 は指数部、 は仮数を基数10で表したものである。

2. MSXの場合
---------

MSX-BASICの演算ルーチンMATHPACKの場合、同様に4バイトで表すが、IEEE 754とは異なり

s: 1
y: 7
x: 24syyy yyyy xxxx xxxx xxxx xxxx xxxx xxxx であり、指数部をバイナリ、仮数部をBCDで表現する。 そのため、有効数字は正確に10進で6桁で、指数は±63乗である。 他のパソコンのBASICはほとんどが2進での演算であった。

2進小数の主な利点は高速性とメモリ消費にあり、BCD小数の主な利点は現実の10進小数を誤差なく扱える事にあったと言ってよい。 2進小数は原理上0.5や0.25、0.375といったn/以外の値、つまり0.1や0.2といった数値を正確に扱う事が苦手であり高速性よりも正確性が優先されるようなケースにおいてBCD小数は強みを発揮した。

ただし当時のスペックにおいて10進小数は相応に重い処理であり、DAAという10進補正用の専用CPU命令を持っていたZ80であっても複雑な演算になればなるほど2進小数より負荷が増していく。 MSX-BASICにおいて10進小数が採用された理由は、マイクロソフトBASICの中でも比較的後期に開発され新しい実装を試せた事、加えてホームコンピュータという需要の開拓においてより一般性を高めたいアスキー側の狙いがあったと考えられる。

CPUの浮動小数点ユニットとの相互運用を視野に入れなければならない現代の言語では、特殊な用途を除いてBCD小数を目にする機会は少ない。

3. 脚注・出典
--------

4. 関連項目
-------

IEEE 754
浮動小数点数
プリミティブ型
数値的安定性

5. 外部リンク
--------

[1]
[2]
[3]
[4] : IEEEの倍精度・単精度・半精度の間での変換を行うCのソースコードがある。