24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3060  |  回复: 6

舞昂琦VoVo

木虫 (小有名气)

[求助] 如何在lammps里加入新的势函数形式?

想在lammps里加入一种势函数形式,但看了一下lammps里已有的势函数源代码,感觉比较复杂,有点无从下手,不知道哪位有类似的经验,能够指教一下,万分感激~
回复此楼
我不想停止
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zh1987hs: 金币+5, 谢谢 2012-04-10 08:06:43
舞昂琦VoVo: 金币+2, ★★★很有帮助 2012-04-10 13:51:55
lammps的势函数已经非常全面了,你仔细遇到以下lammps手册,说不定你想添加的,人家已经有了呢。http://lammps.sandia.gov/doc/Manual.html

如果没有的话,建议你首先把lammps里面的源程序里面的势函数看明白,打个比方说,你想添加一个新的bond style ,那么你把 源程序里面的bond_fene.cpp 和 bond_harmonic.cpp 两个程序看明白了,再去写你想要添加的应该就比较容易了。你还需要弄明白lammps的整个程序框架。
看看这个吧:http://lammps.sandia.gov/doc/Section_modify.html#mod_14
2楼2012-04-09 22:25:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

quanta

新虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zh1987hs: 金币+4, 谢谢 2012-04-10 08:06:58
舞昂琦VoVo: 金币+2, ★★★很有帮助 2012-04-10 13:52:03
lammps现在已经很庞大了,研究其框架可以找比较早期的版本。等大致了解了之后然后在最新版本上进行代码修改。lammps是用c++写的,虽然不算标准的模块化设计,但相对gromacs的代码来说比较模块了,添加一个新的势基本不影响其他的代码。
读代码从main.cpp开始看,了解主要框架,然后重点参考你和你要加势的最像近倒那一类,比如bond_harmonic,pair_lj等。在这个基础上完成自己的修改。一般是2个文件,一个.h,一个.cpp。你需要完成的方法在bond.h,angle.h,pair.h..里面有。
然后还要修改1个文件: style.h,让lammps知道你的势。
然后就make,debug,。。。
3楼2012-04-10 00:20:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

舞昂琦VoVo

木虫 (小有名气)

引用回帖:
2楼: Originally posted by jiaoyixiong at 2012-04-09 22:25:13:
lammps的势函数已经非常全面了,你仔细遇到以下lammps手册,说不定你想添加的,人家已经有了呢。http://lammps.sandia.gov/doc/Manual.html

如果没有的话,建议你首先把lammps里面的源程序里面的势函数看明白, ...

恩恩,多谢前辈指教,现在还是看程序的能力不够啊~您是怎么一步步修炼看程序的能力的呀,再追问一下~
我不想停止
4楼2012-04-10 13:54:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

舞昂琦VoVo

木虫 (小有名气)

引用回帖:
3楼: Originally posted by quanta at 2012-04-10 00:20:22:
lammps现在已经很庞大了,研究其框架可以找比较早期的版本。等大致了解了之后然后在最新版本上进行代码修改。lammps是用c++写的,虽然不算标准的模块化设计,但相对gromacs的代码来说比较模块了,添加一个新的势基 ...

恩,早期的版本是指fortran版的么,现在准备开始看源代码了,才发现有点没头绪,读代码的能力有待提高啊,恩,多谢指教。。
我不想停止
5楼2012-04-10 13:58:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

quanta

新虫 (小有名气)

★ ★ ★ ★
zh1987hs: 金币+4, 谢谢 2012-04-10 15:16:30
引用回帖:
5楼: Originally posted by 舞昂琦VoVo at 2012-04-10 13:58:00:
恩,早期的版本是指fortran版的么,现在准备开始看源代码了,才发现有点没头绪,读代码的能力有待提高啊,恩,多谢指教。。

当然必须是C++的。如果找不到的话也没太大关系,看到无关紧要的就直接无视。我有07年的老版本,有些命令和现在不兼容,所以直接看现在的版本也无大碍,懂的取舍就好,没必要也没精力看完所有的代码。
另外一个很好的参考就是src下那些独立的包,比如dpd,use_xxx什么的。
另外读代码的能力真的是如你所说是“修炼”出来的,真经就是“多读多写,别无它法”。希望能对你有所帮助。
6楼2012-04-10 14:45:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cgzhang_gg

铜虫 (著名写手)

引用回帖:
2楼: Originally posted by jiaoyixiong at 2012-04-09 22:25:13
lammps的势函数已经非常全面了,你仔细遇到以下lammps手册,说不定你想添加的,人家已经有了呢。http://lammps.sandia.gov/doc/Manual.html

如果没有的话,建议你首先把lammps里面的源程序里面的势函数看明白,打 ...

你好,帮忙看一下这个求助帖子 http://muchong.com/bbs/viewthread.php?tid=5946644  现在,我想就利用lammps里面的命令,先不改程序,是否就可以达到目的?
7楼2013-05-30 08:44:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 舞昂琦VoVo 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +6 yexuqing 2026-04-19 7/350 2026-04-23 12:32 by yexuqing
[基金申请] 国自然面上和省基金B类撒花 +18 花田半亩~白 2026-04-21 18/900 2026-04-23 11:31 by 12021227
[考研] 有没有学校收留 +3 蒋昌鹏qtj 2026-04-20 3/150 2026-04-22 20:25 by 学员JpLReM
[考研] 312求调剂 +3 山河似你温柔 2026-04-22 3/150 2026-04-22 20:17 by 学员JpLReM
[考博] 华师大读博 +3 xq83 2026-04-22 5/250 2026-04-22 10:42 by xq83
[论文投稿] 急需审稿人!!! +3 陆小果画大饼 2026-04-21 3/150 2026-04-21 23:54 by jzy_123456
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 8/400 2026-04-21 16:26 by 啃面包的小书虫
[考研] 295分求调剂 +6 ?要上岸? 2026-04-17 6/300 2026-04-21 08:18 by Equinoxhua
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
信息提示
请填处理意见