IEEE 754-2008

Der Standard IEEE 754-2008, der frühere Arbeitstitel lautete IEEE 754r, ist eine notwendig gewordene Revision des 1985 verabschiedeten Gleitkommastandards IEEE 754. Der alte Standard war sehr erfolgreich und wurde in zahlreichen Prozessoren und Programmiersprachen übernommen. Die Diskussion über die Revision begann im Jahr 2001; im Juni 2008 wurde der Standard angenommen und im August 2008 verabschiedet.[1]

Hauptziele

Die Hauptziele des verabschiedeten Standards konnten aufgeteilt werden in

  • das Zusammenführen von IEEE 754 und IEEE 854,
  • die Reduktion von Implementierungsalternativen,
  • die Entfernung von Mehrdeutigkeiten der bisherigen IEEE 754,
  • ein zusätzliches kumulierendes Produkt fused multiply-add: FMA(A,B,C) = A·B + C,
  • neben einfacher und doppelter auch Arithmetik mit halber und vierfacher Genauigkeit (zusätzlich zu 32 und 64 Bit auch 16 und 128 Bit),
  • die von der Finanzwirtschaft als notwendig erachteten Dezimalformate (IEEE 854),
  • weitere variable Formate und Austauschformate,
  • min und max mit Spezifikationen für die Spezialfälle ±0 und ±∞ sowie
  • Kosmetik: ab sofort heißt „denormalisiert“ „subnormal“

Der Standard soll Formate und Methoden für Gleitkommaarithmetik sowie eine Mindestqualität definieren.

Formate

Formate umfassen Gleitkommazahlen mit halber (16 Bit), einfacher (32 Bit), doppelter (64 Bit) sowie vierfacher (128 Bit) Genauigkeit. Das Halbformat stellt ein standardisiertes Minifloat dar. Ergänzt werden die Grundformate durch erweiterte (extended) und erweiterbare (neu!) Langzahl-Formate. Ebenfalls neu aufgenommen wurden Datenaustauschformate. Neben der 16/32/64/128-Bit-Darstellungen sind Darstellungen mit einem Vielfachen von 32 Bits definiert.

Dicht gepackte Dezimalformate (DFP, 3 Ziffern in 10 Bit) sind ebenfalls dazugekommen. Sie weichen von klassischen einzelzifferbasierten BCD-Formaten folgendermaßen ab:

  • Die Kapazität der nutzbaren Bits wird gut ausgenutzt, da 3 Dezimalziffern (000...999, 1000 genutzte Werte) in jeweils 10 Bit (0...1023, 1024 mögliche Werte) gespeichert werden. Eine solche Gruppe heißt Declet. Der Verschnitt ist gegenüber klassischen BCD-Zahlen deutlich kleiner. Die letzte Spalte der Tabelle enthält den Informationsgehalt in Bit, der nur geringfügig geringer ist als der Speicherplatz (bei d=7 Mantissenziffern und einem Exponentenwertebereich von emin - emax unter Berücksichtigung der Vorzeichenbits ).
  • Die Verarbeitung der Dezimalziffern in Dreiergruppen kommt der üblichen Gruppierungsgewohnheit (23 223 456; 24 W, 24 kW, 24 MW) entgegen.
  • Die Zahl 0 hat auch das Bitmuster „0000…0“. Allerdings hat 0 eine relativ große Kohorte.
  • Die Zahlen 0 bis 9 eines Declets haben in den 6 führenden Bits eine 0.
  • Die Zahlen 10 bis 99 eines Declets haben in den 3 führenden Bits eine 0.
  • Ungerade Zahlen in Declets können mit Hilfe eines einzelnen Bits erkannt werden.
  • Die 24 unbenutzten Bitmuster ddx11x111x mit dd = 01, 10 oder 11 können leicht identifiziert werden.
  • Die so mit Declets gepackten Zahlen (Densely Packed) sind nicht mehr binär sortierbar, im Gegensatz zu „klassischen BCD-Formaten“.
  • Statt Speicherung in Declets kann die Mantisse auch ganzzahlig binär in einem gleich großen Bitfeld gespeichert werden. Die Bitfeldaufteilung ist im Combinationfield dann anders.
  • Eine Zahl ist nicht eindeutig; mehrere Bitmuster können dieselbe Zahl bezeichnen. Die Menge der Bitmuster einer Zahl heißt Kohorte. Innerhalb einer Kohorte wurde jedoch jeweils eine kanonische Darstellung festgelegt.

Signaling NaNs wurden zur Streichung vorgeschlagen (3. Februar 2003), später aber wieder in den Vorschlag aufgenommen (21. Februar 2003). Eine Signaling NaN ist eine NaN mit gesetztem Bit 7. Darstellungen von existieren und sind leicht erkennbar.

TypSpei-
cher-
bedarf
MantisseExponentInfor-
mations-
gehalt
in Bit
Bits meffektive Bits einer normali­sier­ten Zahl pBits eWertebereichWerte der Ko­horte einer nor­ma­li­sier­ten Zahl
eminemaxBias
b16 (half)016 Bit010011050000−14000015000151 ≤ E ≤ 30016
b32 (single)032 Bit02302408000−126000127001271 ≤ E ≤ 254032
b64 (double)064 Bit0520531100−1022001023010231 ≤ E ≤ 2046064
b128 (quadruple)128 Bit112113150−16382016383163831 ≤ E ≤ 32766128
b256 (octuple)256 Bit23623719−2621422621432621431 ≤ E ≤ 524286256
k = 32j
mit j ≥ 4
00k Bitk + 12 −
rnd(4·ld k)
k + 13 −
rnd(4·ld k)
rnd(4·ld k)
− 13
1 − emax2k−p−1 − 1emax1 ≤ E ≤ 2·emax00k
d32032 Bit020+5(a)07 Ziffern0600−9500960101031,83
d64064 Bit050+516 Ziffern080−38303840398063,73
d128128 Bit110+534 Ziffern12−614361446176127,53
k = 32j
mit j ≥ 1
00k Bit15 k/16 −
10
9 k/32 − 2
Ziffern
k/16 + 41 − emax3·2k/16+3emax +
p − 2
(a) 20+5 in Spalte 3 bedeutet:
  • in den 20 Bits werden 6 Dezimalstellen gespeichert (3 Stellen in jeweils 10 Bit)
  • in den 5 übrigen Bits wird gespeichert:
    • eine weitere Dezimalstelle
    • der Rest des Exponents bei Division durch 3
    • Signalisierungen für NaNs und Infs

Rundungen

Zu den vier alten IEEE-754-Rundungen kommt eine zusätzliche hinzu, so dass folgende Rundungen gefordert werden:

  • vergrößernd (in Richtung +unendlich)
  • verkleinernd (in Richtung −unendlich)
  • betragsverkleinernd (in Richtung 0)
  • bestmöglich und in der Mitte zur nächsten geraden Zahl (to next or to even)
  • bestmöglich und in der Mitte betragsvergrößernd (to next – neu in IEEE 754r, eigentlich nur die klassische Handrechnungsrundung)

Die IEEE 754-Rundung (next even) wurde schon von Carl Friedrich Gauß vorgeschlagen und vermeidet ein statistisches Ungleichgewicht bei längeren Rechnungen zu größeren Zahlen hin.

Ausnahmen

Ausnahmebedingungen und Ausnahmebehandlung werden spezifiziert.

Neue Funktionen sind Prädikatfunktionen (größer gleich) und Operatoren für Maximum und Minimum. Hier wird vor allem über die Ergebnisse bei den Sonderwerten (NaN, Inf) diskutiert.

Dezimalkodierungen

Speicherplatzbedarf
DPDGröße für äquivalente Packed BCDGewinn
032 bit07×4 + 07,58 + 1 bit = 036,58 bit0+4,48 bit
064 bit16×4 + 09,58 + 1 bit = 074,58 bit+10,48 bit
128 bit34×4 + 13,58 + 1 bit = 150,58 bit+22,48 bit

Die primäre Idee hinter der dicht gepackten Dezimaldarstellung ist, dass diese mit extrem wenig (Gatter-)Aufwand in eine klassische BCD-Darstellung für die Mantisse sowie einen binären Exponenten umkodiert werden kann, aber gleichzeitig den Speicherplatz so effizient wie möglich ausnutzt. Die eigentliche Verarbeitung findet dann im klassischen BCD-Format statt, nur beim Lesen und Schreiben von Registern ist eine Umkodierung erforderlich.

Die Kodierung von 32-bit-, 64-bit- und 128-bit-dezimalkodierten Zahlen erfolgt nach folgendem Schema. Für längere Dezimalkodierungen werden für jedes weitere 32-bit-Wort dem Exponenten 2 bit und der Mantisse 30 bit (3× 10 bit) zugeschlagen, so dass unter Beibehaltung des 5-bit-Kombinationsfeldes der Wertebereich des Exponenten sich vervierfacht und die Mantisse weitere neun Ziffern erhält.

FormatVorzeichenKombinationsfeldrestl. Exponentrestliche Mantisse
032 bit1 bit5 bits06 bits020 bits
sm m m m mxxxxxxbbbbbbbbbb bbbbbbbbbb
064 bit1 bit5 bits08 bits050 bits
sm m m m mxxxxxxxxbbbbbbbbbb bbbbbbbbbb bbbbbbbbbb bbbbbbbbbb bbbbbbbbbb
128 bit1 bit5 bits12 bits110 bits
sm m m m mxxxxxxxxxxxxbbbbbbbbbb bbbbbbbbbb bbbbbbbbbb bbbbbbbbbb bbbbbbbbbb
bbbbbbbbbb bbbbbbbbbb bbbbbbbbbb bbbbbbbbbb bbbbbbbbbb
bbbbbbbbbb
0: positiv
1: negativ
Kodierung der MSBs
nach Tabelle 1
binäre
Kodierung
Jedes Declet ist nach Tabelle 2 kodiert und liefert drei weitere Ziffern.
VorzeichenZiffer 1Expon. MSBExpon. LSBZiffer 2Ziffer 3Ziffer 4Ziffer 5Ziffer 6Ziffer 7Ziffer 8Ziffer 9...

Die Zahl besteht aus

  • einem Vorzeichen: dieses wird im Vorzeichenbit s gespeichert.
  • einem Exponenten, der seinen Wertebereich von emin ... emax unter Zuhilfenahme eines Bias auf die Werte 0 ... 3 · 2e − 1 = (0 ... 2) · 2e + (0 ... 2e − 1) abbildet. Die oberen drei Zustände werden im Kombinationsfeld, die restlichen e bit binär im restlichen Exponenten gespeichert.
  • einer Mantisse, die aus p = 3 · n + 1 Ziffern besteht. Die höchstwertige Ziffer wird im Kombinationsfeld, die restlichen 3 · n Ziffern werden in Dreiergruppen in der restlichen Mantisse gespeichert.

Zur Dekodierung und Kodierung werden folgende Kodiertabellen benötigt:

Tabelle 1: Kodierregeln für das Kombinationsfeld der MSBs des Exponenten und der Mantisse
KombinationsfeldMSBs desKod.
Wert
Beschreibung
m4m3m2m1m0Exp.Mant.
00abc000abc(0-7)Ziffer bis 7
01abc010abc
10abc100abc
1100c00100c(8-9)Ziffer größer 7
1101c01100c
1110c10100c
11110±Infinity
11111NaN
Bemerkung
Das Vorzeichenbit von NaNs wird ignoriert. Das MSB des restlichen Exponenten bestimmt, ob das NAN quiet oder signaling ist.
Tabelle 2: Kodierregeln für die Declets der dichtgepackten dezimalen Ziffern der restlichen Mantisse[2]
DPD kodierter WertDezimalziffern
b9b8b7b6b5b4b3b2b1b0d2d1d0Kodierter WertBeschreibung
abcdef0ghi0abc0def0ghi(0–7) (0–7) (0–7)drei Ziffern bis 7
abcdef100i0abc0def100i(0–7) (0–7) (8–9)zwei Ziffern bis 7,
eine größer 7
abcghf101i0abc100f0ghi(0–7) (8–9) (0–7)
ghcdef110i100c0def0ghi(8–9) (0–7) (0–7)
ghc00f111i100c100f0ghi(8–9) (8–9) (0–7)eine Ziffer bis 7,
zwei Ziffern größer 7
dec01f111i100c0def100i(8–9) (0–7) (8–9)
abc10f111i0abc100f100i(0–7) (8–9) (8–9)
??c11f111i100c100f100i(8–9) (8–9) (8–9)drei Ziffern größer 7
Hinweis
Da im Gegensatz zur Binärdarstellung, in der durch Normalisierung und Weglassen des MSBs eine Normalisierung erzwungen wird, keine Normalisierung erzwungen wird und die Ziffer 0 als höchstwertige Ziffer verfügbar ist, sind Zahlen nicht eindeutig kodierbar.

Dezimale Gleitkommazahlen in der Praxis

Die Probleme von dezimalen Gleitkommazahlen sind unter anderem:

  • Sowohl im Binär- wie im Dezimalformat sind die meisten Zahlen nicht präzise darstellbar. Nach wenigen Rechenschritten sind die meisten Berechnungen unpräzise. Eine Währungsumrechnung oder das Abziehen der Umsatzsteuer reicht aus.
  • Für die meisten angegebenen Probleme gibt es einfachere und gleichzeitig leistungsfähigere Lösungen. Für Finanzaufgaben steht unter .NET z. B. der Datentyp System.Decimal zur Verfügung, der Ganzzahlen mit Beträgen bis 79.228.162.514.264.337.593.543.950.335 exakt darstellen kann.
  • Sie stellt eine weitere Fehlerquelle für Hardware (zusätzliche Logik) und Software (Konvertierfehler) dar.

Die Ergebnisse sind:

  • Dezimale Gleitkommazahlen sind standardisiert, aber auch nach 15 Jahren nicht in fester Hardware verfügbar. Man kann sie in Software, in FPGAs und in ASICs implementieren, aber selbst darüber halten sich die Publikationen in Grenzen und sind meist auf Addition und Subtraktion beschränkt.
  • Die Dezimalformate werden hauptsächlich von der Finanzwirtschaft gefordert, aber sobald man genauer hinschaut, nicht benötigt. Festkommadarstellungen auf Basis der kleinsten Verrechnungseinheit und 64-bit-Ganzzahlen decken gegenüber Decimal64 einen 922× so großen Wertebereich exakt ab (−92.233.720.368.547.758,08...+92.233.720.368.547.758,07 gegenüber −99.999.999.999.999,99...+99.999.999.999.999,99). Sie können allerdings keine noch größeren Werte mit dann verminderter Genauigkeit darstellen noch können sie kleinere Beträge genauer darstellen.

Sinnvoll sind sie:

  • uneingeschränkt als Austauschformate, wenn die genaue Repräsentation von Dezimalwerten erforderlich ist.

Hier prallen zwei gegensätzliche Standpunkte aufeinander.

  • Auf der einen Seite werden die Speicher-, Rechenzeit- und Kosten-Vorteile, sowie die gleichmäßigere Zahlenverteilung eines dualen Formates herausgestellt.
  • Auf der anderen Seite wird argumentiert, dass exakte Ergebnisse (meist sind Ergebnisse wie bei Handrechnungen gemeint) nur mit Dezimalarithmetik möglich sind und in Zeiten schneller Prozessoren und billiger Speicher die Nachteile nicht mehr ins Gewicht fallen.

William Kahan hat behauptet, dass duale Arithmetik in Zukunft kaum noch eine Rolle spielen wird.

“Why is decimal floating-point hardware a good idea anyway? Because it can help our industry avoid errors designed not to be found.”

„Warum ist dezimale Gleitkommahardware auf jeden Fall eine gute Idee? Weil sie unserer Industrie hilft die Fehler zu vermeiden, die verfahrensbedingt nicht gefunden werden können.“

William Kahan: Floating-Point Arithmetic Besieged by “Business Decisions”[3]

Er übersieht dabei aber, dass

  • Gepackte Dezimalformate zusätzliche Chipfläche benötigen, eine geringere Effizienz aufweisen und langsamer sind.
  • Rechenleistung jeder Größenordnung neue Aufgabenbereiche eröffnet und es trotzdem immer wieder Aufgaben geben wird, für die sie nicht ausreicht.
  • Es niemals so viel Rechenleistung geben wird, dass man freiwillig auf diese verzichten würde.
  • Je komplexer die Rechnung, desto weniger interessiert es jemanden, ob diese dezimal exakt darstellbar ist. Nur wenigen auserwählten Zahlen wird die Ehre zuteil, von einem Menschen im Dezimalsystem eingetippt zu werden oder von einem Menschen im Dezimalsystem gelesen zu werden.

Weblinks

Einzelnachweise

  1. IEEE 754-2008: Standard for Floating-Point Arithmetic. IEEE Standards Association, 2008, doi:10.1109/IEEESTD.2008.4610935
  2. Michael F. Cowlishaw: A Summary of Densely Packed Decimal encoding. IBM, 13. Februar 2007, archiviert vom Original am 24. September 2015; abgerufen am 7. Februar 2016.Vorlage:Cite web/temporär
  3. William Kahan: Floating-Point Arithmetic Besieged by “Business Decisions”. (PDF; 174 kB) IEEE-Sponsored ARITH 17 Symposium on Computer Arithmetic, 5. Juli 2005, S. 6 von 28, abgerufen am 19. Februar 2020 (englisch).