A model of the relative error in floating point multiplication is developed and is analyzed stochastically for various choices of computer design parameters. These parameters include the base, the type of rounding rule, the number of guard digits, and whether the post-arithmetic normalization shift (if needed) is done before or after rounding. Under the assumption of logarithmic distribution for the fraction (mantissa), the major stochastic conclusions are:1.The average relative error in multiplication increases as the base increases.2.This error is minimized by selecting the machine base to be binary (better yet, binary with a hidden bit) and is rather large for machines with base 16.3.The classical relative error bounds are pessimistic. The average overestimation by those bounds increases as the base increases.ZusammenfassungEs wird ein Modell für den relativen Fehler bei der Gleitkomma-Multiplikation entwickelt unf für verschiedene Kombinationen der Arithmetikparameter stochastisch analysiert. Die Parameter sind die Basis, die Rundungsart, die Anzahl der Schutzstellen, und ob die Normalisierung vor oder nach dem Runden erfolgt. Bei einer angenommenen logarithmischen Verteilung für die Mantisse kommt man zu folgenden Schlüssen:1.Der durchschnittliche relative Fehler bei der Multiplikation wächst mit der Basis.2.Dieser Fehler wird minimal für die Basis 2 (am besten mit verborgenem ersten Bit) und recht groß für die Basis 16.3.Die klassischen Schranken für den relativen Fehler sind pessimistisch. Ihre durchschnittliche Überschätzung wächst mit der Basis.
[1]
Donald Ervin Knuth,et al.
The Art of Computer Programming, Volume II: Seminumerical Algorithms
,
1970
.
[2]
Peter R. Turner.
The Distribution of Leading Significant Digits
,
1982
.
[3]
Richard Goodman,et al.
Effect of guard digits and normalization options on floating point multiplication
,
1977,
Computing.
[4]
Seppo Linnainmaa,et al.
Improved Trailing Digits Estimates Applied to Optimal Computer Arithmetic
,
1979,
JACM.
[5]
Donald Ervin Knuth,et al.
The Art of Computer Programming
,
1968
.
[6]
Richard Goodman,et al.
Convergence Estimates for the Distribution of Trailing Digits
,
1976,
JACM.
[7]
P. Henrici.
Discrete Variable Methods in Ordinary Differential Equations
,
1962
.
[8]
Nai-Kuan Tsao,et al.
On the distributions of significant digits and roundoff errors
,
1974,
CACM.
[9]
Bede Liu,et al.
On Local Roundoff Errors in Floating-Point Arithmetic
,
1973,
JACM.
[10]
Richard Goodman,et al.
Some models of error in floating point multiplication
,
1981,
Computing.
[11]
Pat H. Sterbenz,et al.
Floating-point computation
,
1973
.
[12]
Murlan S. Corrington,et al.
Table erratum: Discrete variable methods in ordinary differential equations (Wiley, New York, 1962) by Peter Henrici
,
1970
.