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语言编写的程序，望高手指教，随便什么建议都可以，感谢感谢了

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

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

附件中文程序代码

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已经更新

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的, 还请指教再次感谢

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

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

xuanyz · 发表于 2019-6-5 10:15:19

再次感谢vvt，多谢了

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