monte carlo 模拟光子传输

vvt · 发表于 2019-6-4 20:40:30

fortran语法没有任何限制。
但可能来自于你使用的操作系统、平台、编译器等。

具体需给出代码以便测试分析。

vvt · 发表于 2019-6-4 21:54:10

没有输入文件，没法测试。

你试试把 pm ， pn 定义为整数。

vvt · 发表于 2019-6-5 09:03:10

pm，pn 改成整数就好了

vvt · 发表于 2019-6-5 09:42:07

110 PN=PN+1.0
改成
110 PN=PN+1

vvt · 发表于 2019-6-5 09:57:02

不定义成整数的话，就是浮点数。
浮点数是有误差的。
在数学上，17000000.0 + 1.0 = 17000001.0
但是浮点数只有7位有效数字，不能表达 17000001.0，就导致结果还是 17000000.0

所以你不管怎么加，都不能满足 IF(PN.GT.PM)
该用整数的时候，一定要用整数类型。因为整数类型是没有误差的（但是有表达的最大范围）
像计个数这种量，应该用整数类型。

[通用算法] monte carlo 模拟光子传输

规矩勋章

元老勋章

新人勋章

水王勋章

热心勋章