Fortran无法在865398.78和-865398.78之间添加？为什么答案是-0.03？

Question

在下面的代码中我加在一起865398.78和-865398.78。我期望得到0，但是我得到-0.03。

源代码：

​​

结果：

x= -865398.780000000  
x+865398.78= -3.000000002793968E-002 

我错了以 “读” 代码或别的什么用途？

Answer 1

数字865398.78在您的代码中以单精度表示。单精度可以通过写

x=x+865398.78_8 

Answer 2

处理大约7显著位数，而你的电话号码有8可以让双精度我会做出这样的回答一个big assumption：那real(8)对应double precision。

您可能会认为您的865398.78在任何地方都表示相同的东西。在真实的源代码中：它是一个默认的真实文字常量，接近865398.78。

当你

x=x+865398.78 

为x双精度，则默认实常数被转换为双精度值。

然而，在读声明

read(10,*)x 

给定输入“-865398.78”然后x需要一个双精度逼近该值。

你的非零答案来自一个事实，即将一个默认的实数/单精度逼近转换为双精度值并不普遍，并且在这种情况下，与初始双精度逼近不同。

最后这个事实在其他问题中有更详细的解释。正如使用x=x+865398.78_8（或更好的，不要使用8作为种类值）的解决方案。