24小时热门版块排行榜    

查看: 3323  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +12 晴天加油 2026-08-26 13/650 2026-08-28 11:32 by skyjudy
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +10 Kittylucky 2026-08-27 11/550 2026-08-28 10:45 by Kittylucky
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见