24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1873  |  回复: 23

Ballendtoo

银虫 (初入文坛)

[求助] 怎么样在Lammps里面固定输出间隔,然后自动输出? 已有1人参与

比如400000步的程序,我希望每跑间隔50000步自动输出一个当时的盒子大小x, y, z 该怎么实现呢?

有做过的同学吗
谢谢大家了
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

ysshj

金虫 (小有名气)


月只蓝: 金币+1, 鼓励交流! 2014-04-16 18:37:17
引用回帖:
3楼: Originally posted by Ballendtoo at 2014-04-11 14:48:05
你好,为什么跑程序出来的结果,lx ly lz会输出 -nan -nan -nan 是势能设置出了问题吗? 求救~...

我只能说这种原因很复杂,导致出现除数很小的情况。
势函数如果按参考文献上的取,感觉一般不会有大错。
要不把in文件,data文件发上来大家看看?
一个人要像一支队伍。
4楼2014-04-11 16:08:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ysshj

金虫 (小有名气)

引用回帖:
18楼: Originally posted by Ballendtoo at 2014-04-22 16:52:03
用了
thermo          10
thermo_style custom step lx ly lz xlo xhi ylo yhi zlo zhi
thermo_modify        flush yes
这一段 每10step 输出lxlylz  但是为什么不输出toteng能量了呢?...

thermo_style 就是根据自己的需要修改输出内容:step lx ly lz xlo xhi ylo yhi zlo zhi,后面加上ke pe etotal就也出能量了。默认是thermo_style custom step
temp epair emol etotal press。详情见手册。

» 本帖已获得的红花(最新10朵)

一个人要像一支队伍。
19楼2014-04-22 17:47:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

ysshj

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
Ballendtoo: 金币+2, ★★★很有帮助, 十分感谢~~~~~我先试试看,确实是需要输出lx ly lz 2014-04-03 19:05:55
月只蓝: 金币+2, 感谢指导! 2014-04-07 09:54:16
Ballendtoo: 金币+3, ★★★★★最佳答案 2014-04-11 19:32:53
thermo         50000
thermo_style custom step lx ly lz xlo xhi ylo yhi zlo zhi
一个人要像一支队伍。
2楼2014-04-03 18:39:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ballendtoo

银虫 (初入文坛)

引用回帖:
2楼: Originally posted by ysshj at 2014-04-03 18:39:21
thermo         50000
thermo_style custom step lx ly lz xlo xhi ylo yhi zlo zhi

你好,为什么跑程序出来的结果,lx ly lz会输出 -nan -nan -nan 是势能设置出了问题吗? 求救~
3楼2014-04-11 14:48:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ballendtoo

银虫 (初入文坛)

引用回帖:
4楼: Originally posted by ysshj at 2014-04-11 16:08:30
我只能说这种原因很复杂,导致出现除数很小的情况。
势函数如果按参考文献上的取,感觉一般不会有大错。
要不把in文件,data文件发上来大家看看?...

in 文件应该没有问题,因为跑一个结构可以正确输出结果,但是我把结构从单层改成双层,就输出nan了,具体的如下:

Ewald initialization ...
  G vector (1/distnace) = -nan
  estimated absolute RMS force accuracy = -nan
  estimated relative force accuracy = -nan
  KSpace vectors: actual max1d max3d = 0 1 13
Memory usage per processor = 5.82992 Mbytes
Step Lx Ly Lz Xlo Xhi Ylo Yhi Zlo Zhi
       0         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan
      10         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan
      20         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan
      30         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan
      40         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan
      50         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan
      60         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan
      70         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan
      80         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan
      90         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan
     100         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan         -nan
Loop time of 0.03101 on 16 procs for 100 steps with 11508 atoms

Ewald initialization ...
  G vector (1/distance) = 0.374164
  estimated absolute RMS force accuracy = 1.64594e-05
  estimated relative force accuracy = 1.14304e-06
  KSpace vectors: actual max1d max3d = 38056 28 92596
Setting up run ...
Memory usage per processor = 43.2059 Mbytes
Step Temp E_pair E_mol TotEng Press Volume
       0         1373         -nan            0         -nan         -nan    182551.45
      50         -nan         -nan            0         -nan         -nan         -nan
Loop time of 0.211481 on 16 procs for 50 steps with 11508 atoms

就是这两段出现了问题,但是跑9000个原子 相同的条件下 就可以输出
5楼2014-04-11 19:29:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ballendtoo

银虫 (初入文坛)

引用回帖:
4楼: Originally posted by ysshj at 2014-04-11 16:08:30
我只能说这种原因很复杂,导致出现除数很小的情况。
势函数如果按参考文献上的取,感觉一般不会有大错。
要不把in文件,data文件发上来大家看看?...

我结构那个文件太长了
6楼2014-04-11 19:33:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ysshj

金虫 (小有名气)


月只蓝: 金币+1, 鼓励交流! 2014-04-19 14:30:37
引用回帖:
6楼: Originally posted by Ballendtoo at 2014-04-11 19:33:36
我结构那个文件太长了...

信息不详啊,data文件变了,in文件一般也得改。
层间作用的势函数是按文献设置的吗?
可以把in文件和data文件当附件上传的,这个就随意了。
一个人要像一支队伍。
7楼2014-04-11 20:45:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ballendtoo

银虫 (初入文坛)

引用回帖:
7楼: Originally posted by ysshj at 2014-04-11 20:45:43
信息不详啊,data文件变了,in文件一般也得改。
层间作用的势函数是按文献设置的吗?
可以把in文件和data文件当附件上传的,这个就随意了。...

好的,我把结构检查一下,等看好了发上来,多谢了
8楼2014-04-12 15:38:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ballendtoo

银虫 (初入文坛)

引用回帖:
7楼: Originally posted by ysshj at 2014-04-11 20:45:43
信息不详啊,data文件变了,in文件一般也得改。
层间作用的势函数是按文献设置的吗?
可以把in文件和data文件当附件上传的,这个就随意了。...

那个结构我放弃了。。太折磨了。。
9楼2014-04-16 14:09:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ysshj

金虫 (小有名气)

引用回帖:
9楼: Originally posted by Ballendtoo at 2014-04-16 14:09:37
那个结构我放弃了。。太折磨了。。...

感觉调试时是如履薄冰啊,指不定那就出错了。祝顺利。
一个人要像一支队伍。
10楼2014-04-16 14:22:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Ballendtoo 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 22专硕求调剂 +6 haoyun上岸 2026-04-11 8/400 2026-04-11 23:21 by labixiaoqiao
[考研] 一志愿211,0703化学305分求调剂 +26 严西西戏 2026-04-06 33/1650 2026-04-11 23:01 by 314126402
[考研] 293求调剂 +8 勇远库爱314 2026-04-06 8/400 2026-04-11 20:25 by 蓝云思雨
[考研] 271求调剂 +21 2261744733 2026-04-11 21/1050 2026-04-11 20:19 by Angchemist
[考研] 085410 273分调剂 +4 X1999 2026-04-09 4/200 2026-04-11 13:05 by pies112
[考研] 考研调剂 +26 硕星赴 2026-04-09 27/1350 2026-04-10 22:24 by 猪会飞
[考研] 263能源动力专硕求调剂 +3 加大号饭盒袋 2026-04-10 3/150 2026-04-10 22:23 by 286640313
[考研] 083200 305分 求二轮调剂 不接受跨专业 +9 Claireyyyy 2026-04-09 10/500 2026-04-10 21:21 by Claireyyyy
[考研] 289 分105500药学专硕求调剂(找B区学校) +6 白云123456789 2026-04-09 8/400 2026-04-10 21:13 by zhouxiaoyu
[考研] 计算机类求调剂,22408-274分 +7 上岸de小虫 2026-04-09 8/400 2026-04-10 19:56 by fxue1114
[考研] 284求调剂 +9 让我上岸吧阿西 2026-04-09 11/550 2026-04-10 19:18 by 靖jing
[考研] 0856专硕求调剂 希望是a区院校 +21 好好休息好不好 2026-04-09 24/1200 2026-04-10 16:58 by luoyongfeng
[考研] 机械专368 有去处吗 +4 种大树 2026-04-10 4/200 2026-04-10 15:31 by jiajinhpu
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +27 努力奋斗112 2026-04-07 30/1500 2026-04-10 15:06 by Kilig0317
[考研] 材料调剂 +5 hzhahg 2026-04-06 5/250 2026-04-10 10:10 by may_新宇
[考研] 求调剂材料科学与工程一志愿985初试365分 +5 材化李可 2026-04-08 5/250 2026-04-09 17:00 by Lilly_Li
[考研] 调剂 +12 月@163.com 2026-04-08 12/600 2026-04-09 14:27 by rl1980
[考研] 259求调剂 +5 就爱吃土豆呀呀 2026-04-07 5/250 2026-04-07 22:40 by JourneyLucky
[考研] 调剂 +4 mcbbc 2026-04-06 5/250 2026-04-07 12:33 by upczlm1989
[考研] 考研调剂 +3 Wwwwwww哇 2026-04-06 3/150 2026-04-06 20:55 by lbsjt
信息提示
请填处理意见