24小时热门版块排行榜

返回列表

本帖产生 2 个模拟EPI ，点击这里进行查看

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

qphll

金虫 (正式写手)

模拟EPI: 10
应助: 18 (小学生)
金币: 2708.7
散金: 2294
红花: 8
帖子: 554
在线: 227.1小时
虫号: 89654
注册: 2005-08-29
性别: GG
专业: 碳素材料与超硬材料

[交流] 【分享】尝试lammps中, 分享中... 已有9人参与

Lesson 1
Loop inside Lammps input script

想要尝试在lammps的input script里面做循环, 结果因为一个小问题, 捣鼓了我几个小时, 这才完全通过测试. 分享一下.

(1) input script里面的循环块

include tempfile
include tempfile2

variable i loop 10
label loopa
   fix 2 all nvt temp ${mytemp} ${mytemp2} 100.0
      run 200
      unfix 2

   next mytemp
   next mytemp2
next i
jump SELF loopa

这个的SELF是让程序执行到这里, 跳回自己, 然后从标签 loopa开始执行. 当然loopa是随便取的, 你可以用CHN来做label.

另外, 在某些情况下, 如果要让c++ rewind, 那么最好在执行脚本里面这样写:

lmp -in script

而不是

lmp < script

否则, 你人品不好的时候, 会出问题, 哈哈.

(2) tempfile 和 tempfile2是在同目录下的另外两个文件. tempfile的文件内容是:

variable mytemp index 500.0 700.0 900.0 1100.0 1300.0 1500.0 1700.0 1900.0 2100.0 2300.0

注意, 只有一行!  至于lammps能读多长的一行, 我还没有测试.

tempfile2的文件内容, 也是一行:

variable mytemp2 index 700.0 900.0 1100.0 1300.0 1500.0 1700.0 1900.0 2100.0 2300.0 2500.0

折腾我的问题是, 我原先在这两个tempfile文件中, 数值之间用逗号分隔, 但是事实上, 是需要用空格分隔的.

(3)  如果你需要做的循环不是很多, 那么不需要额外准备tempfile 和 tempfile2文件. 而只是需要在 input script中这样做:

variable mytemp index 500.0 700.0 900.0 1100.0 1300.0 1500.0 1700.0 1900.0 2100.0 2300.0
#variable mytemp2 index 700.0 900.0 1100.0 1300.0  1500.0 1700.0 1900.0 2100.0 2300.0 2500.0

variable i loop 10
label loopa
   fix 2 all nvt temp ${mytemp} ${mytemp2} 100.0
      run 200
      unfix 2

   next mytemp
   next mytemp2
next i
jump SELF loopa

当然, 这里不需要原先的这两句 include语句了.

总结一下你需要熟悉的命令:

variable, include, jume, next

End of Lesson 1.

ENJOY.

[ Last edited by qphll on 2010-12-3 at 11:54 ]

回复此楼