| 查看: 1586 | 回复: 7 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
LAMMPS源代码的求助 已有1人参与
|
||||
|
给位老师，大家好，最近我在修改lammps中bond_harmonic的源代码，对其中部分不是很清楚，在这里请教一下： bond_harmonic的势函数为E=K（r-r0)^2, 在源代码中 rsq = delx*delx + dely*dely + delz*delz; r = sqrt(rsq); dr = r - r0[type]; rk = k[type] * dr; 其势函数是这样写的，在这里，我不是很懂程序中[type]的意思，各位有知道的能讲解一下吗？万分感谢！ |
回复此楼» 猜你喜欢
全日制（定向）博士
已经有5人回复
-
假如你的研究生提出不合理要求
已经有10人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- lammps势函数问题
已经有2人回复
- lammps 计算粘度 lj单位 转化为 metal单位 出错
已经有7人回复
- 流固耦合方向德国博士，回国的话，国内哪些高校(或中科院)比较对口？
已经有92人回复
- lammps中能否在源代码中自己编辑公式
已经有1人回复
- 关于LAmmps的windows下源程序编译调试问题
已经有7人回复
- 关于towhee 回答zhang_jaj的问题
已经有44人回复
- 如何在lammps里加入新的势函数形式？
已经有6人回复
- 【求助】更改源代码 并行出问题
已经有9人回复
4楼2016-02-23 15:35:20
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
2楼2016-02-23 15:27:18
3楼2016-02-23 15:33:43
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
【答案】应助回帖
|
没法找到，lammps代码不可能每一行都给你注释，每个变量的具体含义只能根据变量的名字和在程序中的调用，结合注释才能推测。 像ebond， fbond都还是比较简单的能根据名字能推测出含义的变量，也就是键的能量和力。如果碰到其他你看了半天代码也推测不出来含义的变量，可以到lammps mailing list去请教别人。lammps developer guide（http://lammps.sandia.gov/doc/Developer.pdf）可以帮助你了解代码的结构，但不会告诉你每个变量每个函数的具体含义。总之要了解代码是要花一些时间的。 |
5楼2016-02-24 00:26:58