monte carlo 模拟光子传输
蒙特卡洛方法 模拟算法测试可以运行大致思路时一次一个光子随机传输，重复很多次并统计结果。当光子数为1千万个光子时，程序运行正常大约耗时10s左右，依次增加光子数 1.1---1.2 ----1.3----1.4----1.5---1.6千万均计算正常耗时性应增加也基本在十几秒左右，但增加到1.7千万或者更大时程序就永远无法运行完毕，不知道是什么问题，fortran对大量循环次数有限制吗，我采用fortran语言编写的程序，望高手指教，随便什么建议都可以，感谢感谢了
fortran语法没有任何限制。
但可能来自于你使用的操作系统、平台、编译器等。
具体需给出代码以便测试分析。 本帖最后由 xuanyz 于 2019-6-5 11:53 编辑
附件中文程序代码 没有输入文件，没法测试。
你试试把 pm ， pn 定义为整数。 本帖最后由 xuanyz 于 2019-6-5 08:41 编辑
抱歉，是我的疏忽，不需要输入文件，注释或者删掉掉代码中的这几行即可
DO 2 i=1,72
read(10,111) SCATT(i)
111 format(F7.5)
2 CONTINUE
DO 4 I=1,36
SINTH(I)=SCATT(I)
CSINTH(I)=SCATT(I+36)
4 CONTINUE
附件Upscat.for已经更新
pm，pn 改成整数就好了 感谢vvt,
已经按照你的思路将pm,pn改成整数
修改如下INTEGER ZSTART,ZFINSH,BOTTOM,pm,pn
测试光子数pm=16000000正常
但光子数增加到pm=17000000时 程序依然永远无法运行完毕？？？？
难道是我的系统或编译器的限制
我用的是win7 64位 vs2010 结合Intel® Parallel Studio XE 2011编译fortran的, 还请指教 再次感谢 110 PN=PN+1.0
改成
110 PN=PN+1 再次感谢 vvt ,这次
pm=17000000测试通过了，
还想弱弱的问一下，为何将pm，pn 改成整数就可以呢，原因是什么呢，方便解释一下嘛 谢谢 不定义成整数的话，就是浮点数。
浮点数是有误差的。
在数学上，17000000.0 + 1.0 = 17000001.0
但是浮点数只有7位有效数字，不能表达 17000001.0，就导致结果还是 17000000.0
所以你不管怎么加，都不能满足 IF(PN.GT.PM)
该用整数的时候，一定要用整数类型。因为整数类型是没有误差的（但是有表达的最大范围）
像计个数这种量，应该用整数类型。
页:
[1]
2