Fortran中实数的精确问题

Question

我一直在尝试使用Fortran作为我的研究项目,使用GNU Fortran编译器(gfortran),最新版本,但我在处理实数时遇到了一些问题.如果您有代码:

program test
implicit none

real :: y = 23.234, z

z = y * 100000
write(*,*) y, z

end program

你会得到输出:

23.23999    2323400.0 

我发现这很奇怪.有人能告诉我这里到底发生了什么吗？看看z我可以看到y确实保留了它的精度,所以对于我认为应该不成问题的计算.但是为什么y的输出与我指定的值不完全相同,我能做些什么才能使它完全相同？

Answer 1

这不是问题 - 你看到的只是计算机中数字的浮点表示.计算机无法准确处理实数,只能接近它们的近似值.关于这一点的一个很好的解读可以在这里找到:每个计算机科学家应该知道浮点运算.