monte carlo 模拟光子传输

xuanyz · 发表于 2019-6-4 20:19:54

蒙特卡洛方法模拟算法测试可以运行大致思路时一次一个光子随机传输，重复很多次并统计结果。

当光子数为1千万个光子时，程序运行正常大约耗时10s左右，依次增加光子数 1.1---1.2 ----1.3----1.4----1.5---1.6千万均计算正常耗时性应增加也基本在十几秒左右，但增加到1.7千万或者更大时程序就永远无法运行完毕，不知道是什么问题，fortran对大量循环次数有限制吗，我采用fortran语言编写的程序，望高手指教，随便什么建议都可以，感谢感谢了

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

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

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

xuanyz · 发表于 2019-6-4 21:20:14

本帖最后由 xuanyz 于 2019-6-5 11:53 编辑

附件中文程序代码

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

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

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

xuanyz · 发表于 2019-6-4 23:07:45

本帖最后由 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已经更新

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

pm，pn 改成整数就好了

xuanyz · 发表于 2019-6-5 09:31:52

感谢vvt,
已经按照你的思路将pm,pn改成整数
修改如下 INTEGER ZSTART,ZFINSH,BOTTOM,pm,pn
测试光子数pm=16000000正常
但光子数增加到pm=17000000时程序依然永远无法运行完毕？？？？
难道是我的系统或编译器的限制
我用的是win7 64位 vs2010 结合Intel® Parallel Studio XE 2011编译fortran的, 还请指教再次感谢

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

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

xuanyz · 发表于 2019-6-5 09:53:33

再次感谢 vvt ,这次
pm=17000000测试通过了，
还想弱弱的问一下，为何将pm，pn 改成整数就可以呢，原因是什么呢，方便解释一下嘛谢谢

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

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

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

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

规矩勋章

元老勋章

新人勋章

水王勋章

热心勋章