主题
Search

浮点运算


简单来说,浮点运算是指通过任何数量的自动化设备对浮点表示执行的算术运算。

传统上,这个定义的措辞仅适用于对实数的浮点表示(即,浮点数集合的有限元素)执行的算术运算,尽管几种额外的浮点数据类型,包括有符号无穷大和 NaN,也通常被允许作为此类函数的输入。

尽管定义简洁,但值得注意的是,最广泛采用的计算标准将几乎整个浮点理论都归为“浮点运算”之下。这种广度的原因之一源于这样一个事实:任何浮点表示都只能表示实数连续统 R 的有限子集;这种有限性带来了一系列意想不到的障碍,其中最主要的是实数算术的某些性质(例如,加法的结合律)有时对浮点数不成立(IEEE 计算机协会 2008)。因此,任何对浮点运算和/或代数的全面处理都必须在讨论实际操作本身之前,解决包括浮点数表示、舍入等在内的众多注意事项。

截至 2014 年,最常用的浮点运算标准是 IEEE 754-2008 浮点算术标准(简写为 IEEE 754-2008,此后简称为 IEEE 754)。该框架是对其前身 IEEE 754-1985 的大规模修订,包含一套内置的指南,几乎详细说明了浮点理论的方方面面。特别是,IEEE 754 详细讨论了浮点理论的以下方面:

1. 浮点表示和格式。

2. 浮点表示的属性,包括浮点数的舍入。

3. 浮点表示的算术和代数运算。

4. 无穷大、非数字 (NaNs)、符号和异常。

上述许多主题在标准文档的多个章节中进行了讨论(IEEE 计算机协会 2008)。

IEEE 754 对浮点表示定义的“必需”算术运算是加法、减法、乘法、除法、平方根和融合乘加(由 (x,y,z)|->(x×y)+z 定义的三元运算);这些是必需的,因为遵守该框架要求这些运算在整个过程中都支持正确的舍入。该框架内还提供了许多其他“推荐”运算,其中一些本质上是算术运算;这些是推荐的,因为框架并未严格要求支持它们。最后,请注意,该框架既包含一组实用函数(也可以被认为是算术函数),即 copynegateabs,也包含许多为向量值输入定义的密切相关函数(IEEE 计算机协会 2008,第 46-47 页)

运算函数可能的异常
expe^x[-infty,+infty]溢出;下溢
expm1e^x-1[-infty,+infty]溢出;下溢
exp22^x[-infty,+infty]溢出;下溢
exp2m12^x-1[-infty,+infty]溢出;下溢
exp1010^x[-infty,+infty]溢出;下溢
exp10m110^x-1[-infty,+infty]溢出;下溢
loglog_e(x)=ln(x)[0,+infty]除以零(如果 x=0);无效操作(如果 x<0
log2log_2(x)[0,+infty]除以零(如果 x=0);无效操作(如果 x<0
log10log_10(x)=log(x)[0,+infty]除以零(如果 x=0);无效操作(如果 x<0
logp1log_e(x+1)[-1,+infty]除以零(如果 x=-1);无效操作(如果 x<-1);下溢
log2p1log_2(x+1)[-1,+infty]除以零(如果 x=-1);无效操作(如果 x<-1);下溢
log10p1log_10(x+1)[-1,+infty]除以零(如果 x=-1);无效操作(如果 x<-1);下溢
hypot(x,y)sqrt(x^2+y^2)[-infty,+infty]^2溢出;下溢
rSqrt1/sqrt(x)[0,+infty]无效操作(如果 x<0);除以零(如果 x=+/-0
compound(x,n)(1+x)^n[-1,+infty]×Z无效操作(如果 x<-1
rootn(x,n)x^(1/n)[-infty,+infty]×Z无效操作(如果 n=0x<0n 为偶数);溢出/下溢(如果 n=-1
pown(x,n)x^n[-infty,+infty]×Z多种情况
pow(x,y)x^y[-infty,+infty]^2多种情况
powr(x,y)x^y[0,+infty]×[-infty,+infty]多种情况

上表总结了 IEEE 754 中推荐的算术运算。请注意,标记为“多种情况”的异常的具体细节在 IEEE 754 文档中详细说明(IEEE 计算机协会 2008,第 43-45 页)。

如上所述,即使是一些基本的要求的算术运算符,在浮点表示和舍入的情况下,其行为也是不可预测的。这源于这样一个事实:IEEE 754 中假定“正常”算术运算具有无限精度,而浮点加法、减法、乘法和除法的值(分别符号地写为  direct sum ,  circleminus ,  tensor , 和 ⌀)是通过对以公共指数表示的浮点数执行“正常”运算 +, -, ×, 和 /,然后将结果舍入到固定位数的有效数字(通过所谓的首选指数)。因此,在计算的算术和/或舍入步骤中都可能发生精度损失、溢出和下溢。例如,x=1.234567×10^5y=9.876543×10^(-3) 相加的结果正好是

 x+y=(1.234567+0.0000009876543)×10^5=1.23456709876543×10^5.
(1)

另一方面,在基数b=10 且精度为 7 位的框架中,浮点加法返回的值将是

 x direct sum y=1.234567×10^5=x.
(2)

类似地,给定 x=1.234571467×10^5y=1.23457129×10^5,可以得到

 x-y=(1.234571467-1.234571419)×10^5=0.000000048×10^5=0.0048
(3)

使用上面假设的 7 位精度。但是,可以得到

 x circleminus y=0.0,
(4)

因此导致完全没有精度。请注意,在像这样的极端情况下,实现 IEEE 754 的系统不会实际产生 +/-0 作为结果:特别是,这种情况会触发下溢警告。关于上述其他算术函数的详细信息和注意事项可以在文档中找到(IEEE 计算机协会 2008,第 5 节和第 9 节)。

浮点运算在本质上是算术运算和代数/三角运算之间存在一些区别:后一种运算通常归为浮点代数的范畴。


另请参阅

算术, 偏移指数, 浮点代数, 浮点指数, 浮点规格化数, 浮点数, 浮点首选指数, 浮点量子, 浮点表示, IEEE 754-2008, 区间算术, NaN, 安静 NaN, 信令 NaN, 有效数, 次正规数

此条目由 Christopher Stover 贡献

使用 Wolfram|Alpha 探索

参考文献

Goldberg, D. "What Every Computer Scientist Should Know About Floating-Point Arithmetic." ACM Comput. Surv. 23, 5-48, March 1991. http://docs.sun.com/source/806-3568/ncg_goldberg.html.Hauser, J. R. "Handling Floating-Point Exceptions in Numeric Programs." ACM Trans. Program. Lang. Sys. 18, 139-174, 1996. http://www.jhauser.us/publications/HandlingFloatingPointExceptions.html.IEEE Computer Society. "IEEE Standard for Floating-Point Arithmetic: IEEE Std 754-2008 (Revision of IEEE Std 754-1985)." 2008. http://ieeexplore.ieee.org/xpl/articleDetails.jsp?arnumber=4610935.Severance, C. (Ed.). "IEEE 754: An Interview with William Kahan." Computer, 114-115, Mar. 1998.Stevenson, D. "A Proposed Standard for Binary Floating-Point Arithmetic: Draft 8.0 of IEEE Task P754." IEEE Comput. 14, 51-62, 1981.

在 Wolfram|Alpha 上引用

浮点运算

请引用本文为

Stover, Christopher. "浮点运算。" 来自 MathWorld--Wolfram Web 资源,由 Eric W. Weisstein 创建。 https://mathworld.net.cn/Floating-PointArithmetic.html

主题分类