24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wangqingbo

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 817
红花: 1
帖子: 276
在线: 93.1小时
虫号: 337314
注册: 2007-04-02
性别: GG
专业: 数论

[交流] 【求助】蒙卡编程问题已有12人参与

前段时间，用蒙特卡罗编程计算狭缝孔的吸附，遇到一个问题，请高手指点：
用的是LJ12-6作用势，在随机插入分子后，计算每个分子与其他分子及与孔壁作用的能量时，有的能量值比较大，将该能量带入到是否接受插入的判定法则后，导致数据溢出，用的是C++语言，已经声明long double 了还是溢出，不知如何解决。

哎，现在对编程都没信心了

回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

1楼 2010-06-05 11:33:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

mljphy

铁虫 (正式写手)

应助: 5 (幼儿园)
金币: 280
散金: 204
红花: 7
帖子: 311
在线: 89.8小时
虫号: 977325
注册: 2010-03-21
性别: GG
专业: 凝聚态物性 II ：电子结构

确实有这种情况

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

引用回帖:

Originally posted by wangqingbo at 2010-06-11 13:21:14:
谢谢各位指点，我感觉double精度应该就可以，但是在插入粒子的时候，计算得到的插入前后的能量值，如果该值大于零，然后将该值带入到插入的判定准则里面，这个时候判定准则有指数项，所以只要能量变化大于1000就 ...

我也出现这种情况，最后模拟之后反而能量有所升高。
确实是前后能量差太大，而玻尔兹曼因子中指数的分母太小，也不知道如何解决！