倍精度浮動小数点数

倍精度浮動小数点数（ばいせいどふどうしょうすうてんすう、英: double-precision floating-point number）は、64ビットの浮動小数点数表現である。 倍精度とは単精度に対する表現であり、これは32ビットを1ワードとする32ビットアーキテクチャをもとにしている。 昔[いつ?]のFORTRANでは、単精度（REAL型）よりも精度が高ければ倍精度（DOUBLE PRECISION型）を名乗ることができた（そもそもワードの長さも浮動小数点のフォーマットも機種ごとにまちまちだった）。IBMのSystem/360で採用され大型機の事実上の標準となった、指数の基数が16の浮動小数点形式は、32ビット単精度では最悪の場合の精度が十進で6桁程度となり、技術計算では倍精度以上を使わねばならないという問題があった。（注：FORTRANではREAL型が1ワード、DOUBLE PRECISION型が2ワードという前提だった） 浮動小数点算術に関する標準であるIEEE 754では、単精度は32ビット（4オクテット）、倍精度は64ビット（8オクテット）である。いずれにしろ、「倍」というのは、精度に関係する仮数部（後述）の長さが正確に2倍である、といったような意味ではなく、全体の長さが2倍である所から来ているので、実際の所「倍精度」というのはかなり大雑把な言い方に過ぎない。 倍精度浮動小数点数は単精度浮動小数点数に比べて性能や帯域幅のコストがかかるが、表現できる数値の範囲が広いためPCでもよく使われている。浮動小数点方式一般に言えることであるが（単精度の場合と同様）指数部に使われる部分があるために、同じサイズの固定長整数が表現可能な値のうち、ある範囲より絶対値が大きいものについては、その全てを正確に表現することはできない。IEEE 754:2008でbinary64として標準化されたパラメータは次の通りである。 符号ビット ( s i g n {\displaystyle sign} ): 1 ビット 指数部 ( e x p o n e n t {\displaystyle exponent} ): 11 ビット 仮数部 ( f r a c t i o n {\displaystyle fraction} ): 52 ビット 指数部が全て 0（ゼロ、および非正規化数）の場合を除き、仮数部で表現されるビットパターンのさらにひとつ上の桁に暗黙の 1 のビットがあるとみなす表現法（いわゆる「ケチ表現」）により、通常の数（正規化された数）の精度は、53ビット相当（十進に直すと約16桁、 53 log 10 ⁡ ( 2 ) ≈ 15.955 {\displaystyle 53\log _{10}(2)\approx 15.955} ）である。ビットのレイアウトは以下のようになる（具体的なレイアウトは、一般的なバイトアドレッシングのコンピュータにおいて、エンディアンの違いなどによる影響を受けることがあるが、それについては必ずしも標準化されていない）。 64ビットの「倍精度」データで表される実際の数値は、符号を s i g n {\displaystyle sign} 、バイアスのある指数部を e {\displaystyle e} 、52ビットの仮数部のビット列を b n {\displaystyle b_{n}} とすると、 = ( − 1 ) s i g n ( 1. b − 1 b − 2 ⋯ b − 52 ) 2 × 2 e − 1023 {\displaystyle =(-1)^{sign}(1.b_{-1}b_{-2}\cdots {}b_{-52})_{2}\times 2^{e-1023}} となり、より正確に表すと次のようになる。 v a l u e = ( − 1 ) s i g n ( 1 + ∑ i = 1 52   b − i 2 − i ) × 2 ( e − 1023 ) {\displaystyle value=(-1)^{sign}(1+\sum _{i=1}^{52}\ b_{-i}2^{-i})\times 2^{(e-1023)}} 252=4503599627370496 と 253=9007199254740992 の間で表現できる数値は正確に整数に対応している。253 から 254 までの範囲では、常にその2倍となるので、偶数しか表現できない。逆に 251 と 252 の間の範囲では間隔が0.5になる。 数値の小数部の間隔は 2n から 2n+1 の範囲で 2−52 から 2−53 までである。従って最も近い値への数値の丸め誤差の最大（計算機イプシロン）は 2−53 となる。 Emin (00116進) = −1022 E (50) = −973 Emax (7FE16進) = 1023 指数部バイアス (3FF16進) = 1023 指数部バイアスは、エクセスNとも言う。詳しくは符号付数値表現を参照されたい。真の指数値は、指数部の値から指数部バイアスを引いた値となる。 00016進 と 7FF16進 は予約された指数値である。 00016進 は 0（仮数部も0）と非正規化数（仮数部が0でない）を表現するのに使われる。 7FF16進 は無限大（仮数部が0）やNaN（仮数部が0でない）を表現するのに使われる。 従って、全てのビットパターンが符号として意味がある。これらの例外を除くと、倍精度浮動小数点数は次のように表される。 ( − 1 ) sign × 2 exponent − exponent bias × 1. mantissa {\displaystyle (-1)^{\text{sign}}\times 2^{{\text{exponent}}-{\text{exponent bias}}}\times 1.{\text{mantissa}}} より詳細な例として 3FD5 5555 5555 555516進 の場合は以下の通りに考える。 MSX-BASICの演算ルーチンMATHPACKの場合、同様に8バイトで表すが、IEEE 754とは異なり s（符号ビット）: 1 y（指数部の幅）: 7 x（仮数部の幅）: 56 syyy yyyy xxxx xxxx xxxx xxxx xxxx xxxx xxxx xxxx xxxx xxxx xxxx xxxx xxxx xxxx であり、指数部をバイナリ、仮数部をBCDで表現する。 そのため、有効数字は正確に10進で14桁で、指数は±63乗である。 他のパソコンのBASICはほとんどが2進での演算であった（他にBCDを採用した機種としてFP-1000がある）。 2進小数の主な利点は高速性とメモリ消費にあり、BCD小数の主な利点は現実の10進小数を誤差なく扱える事にあったと言ってよい。 2進小数は原理上0.5や0.25、0.375といったn/(2^m)以外の値、つまり0.1や0.2といった数値を正確に扱う事が苦手であり(表示上は補正されているが、10進小数を2進小数に変換すると殆どの場合循環小数となり、微小な切り捨てが発生する)高速性よりも(10進の世界での)正確性が優先されるようなケースにおいてBCD小数は強みを発揮した(各種利息計算など)。 ただし当時のスペックにおいて10進小数は相応に重い処理であり、DAAという10進補正用の専用CPU命令を持っていたZ80であっても複雑な演算になればなるほど(例えば三角関数や対数関数のような数学関数)2進小数より負荷が増していく。 MSX-BASICにおいて10進小数が採用された理由は、マイクロソフトBASICの中でも比較的後期に開発され新しい実装を試せた事、加えてホームコンピュータという需要の開拓においてより一般性を高めたいアスキー側の狙いがあったと考えられる。 CPUの浮動小数点ユニットとの相互運用を視野に入れなければならない現代の言語では、特殊な用途を除いてBCD小数を目にする機会は少ない。 IEEE 754 浮動小数点数 半精度浮動小数点数 (16bit) 単精度浮動小数点数 (32bit) 倍精度浮動小数点数 (64bit) 拡張倍精度浮動小数点数 (80bit) 四倍精度浮動小数点数 (128bit)