Fortran语言中的浮点算法
我继承了一些我试图理解的Fortran代码。它在很多地方使用实际变量，我认为，它不应该是这样的，但是也许我误解了FORTRAN中的工作原理（与我更熟悉的C++相比），因此这个问题。 因此，所讨论的变量本质上是“分类值”、“因子”或“枚举”，这取决于您如何看待/想要调用它。它们是数据类型REAL，只能接受有限数量的预定义整数值。所以说变量a只能是值1、2或3。这些值从外部文件中读取；在这些外部文件中，它们表示为整数，因此不存在“外部数据源中的舍入问题”之类的情况 然而，在代码中，它从不进行直接比较，总是进行大于/小于检查。因此，与其这样做Fortran语言中的浮点算法,fortran,Fortran,我继承了一些我试图理解的Fortran代码。它在很多地方使用实际变量，我认为，它不应该是这样的，但是也许我误解了FORTRAN中的工作原理（与我更熟悉的C++相比），因此这个问题。 因此，所讨论的变量本质上是“分类值”、“因子”或“枚举”，这取决于您如何看待/想要调用它。它们是数据类型REAL，只能接受有限数量的预定义整数值。所以说变量a只能是值1、2或3。这些值从外部文件中读取；在这些外部文件中，它们表示为整数，因此不存在“外部数据源中的舍入问题”之类的情况 然而，在代码中，它从不进行直接比较
if (a == 1) then
if (a > 0.9 .and. a < 1.1) then
如果（a>0.9.和.a<1.1），则
b = 1.5
a = REAL(INT(b))
if (a > 0.9 .and. a < 1.1) then
b=1.5
a=实（INT（b））
如果（a>0.9.和.a<1.1），则
有意义吗？只有在不使用实际值执行任何操作的情况下（只需指定一个值并将相等性与指定的相同文字、相同种类的参数进行比较），才不需要公差 问题是，对于实变量a，例如：
Real a
a = 2
a == 2
.true.a / b * b == 2
（或任何其他操作）不保证为.true。只有在不使用实值执行任何操作的情况下（只需指定一个值并将相等性与指定的相同文字、相同种类的参数进行比较），才不需要公差 问题是，对于实变量a，例如：
Real a
a = 2
a == 2
.true.a / b * b == 2
（或任何其他操作）不保证为.true。如果出于某种原因变量必须保持为实，则可以使用 比较中的内在Fortran函数NINT（最接近的整数）：
if( nint(a) == 1 ) then
....
如果出于某种原因，变量必须保持为实，则可以使用 比较中的内在Fortran函数NINT（最接近的整数）：
if( nint(a) == 1 ) then
....
你是对的：“永远不要在没有公差的情况下比较实数。”IEEE浮点数不能是精确的。你不能用二进制表示0.1，也不能用十进制表示1/3。FORTRAN中的实数是IEEE浮点数；与整数非常不同。FORTRAN的实数遵循模型形式
s*b^e*\sum{k=1}{p}f\k*b^{-k}s*b^e*\sum{k=1}{p}f_k*b^{-k}