24小时热门版块排行榜    

查看: 3315  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

舞昂琦VoVo

木虫 (小有名气)

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

想在lammps里加入一种势函数形式,但看了一下lammps里已有的势函数源代码,感觉比较复杂,有点无从下手,不知道哪位有类似的经验,能够指教一下,万分感激~
回复此楼
我不想停止
已阅   回复此楼   关注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的回帖
查看全部 7 个回答

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的回帖

舞昂琦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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 7/350 2026-08-25 03:24 by BZKMTicpDhFj
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-25 01:47 by BZKMTicpDhFj
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见