24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 600  |  回复: 4

zhenyanwo

管理员

优秀!!有木有!!!优秀!!有木有!!!优秀!!有木有!!!优秀!!有木有!!!

[求助] 大家帮忙看看像这样的文件如何提取想要的内容。 已有2人参与

- That atom                    :               272
                 E-Eref  m_perp    ftot        fpar       fperp       eigen       delr  npart evalf   a1
                 ( eV )                     ( eV/Ang )              ( eV/Ang**2 )
   0  K=   0     0.0120  1 1      0.6542     -0.5858      0.2911      0.0000      0.054   0     2   0.00
   1  K=   1     0.2074  2 2      2.2216     -2.1115      0.6907      0.0000      0.290   1     5   0.00
   ...
SADDLE 1001 CONVERGED |ret  20134 |delta energy=    0.0255 |force_(tot,par,perp)=     0.0854    0.0217    0.0000 |eigenval=  -3.5023 |npart=    1 |delr=    0.343 |evalf=   327 |
  - Configuration stored in file :           sad1001
  - Total energy Saddle (eV)     : -4.1048037236E+03
RELAXATION
  ...
139  M=   5     0.0116  0 0      0.2256      0.1374      0.0000      0.0000      0.377   2   333   0.00
144  M=  10     0.0093  0 0      0.2595      0.1610      0.0000      0.0000      0.380   2   338   0.00
  ...
- Configuration stored in file :           min1001
  - Total energy Minimum (eV)    : -4.1052381485E+03
MINIMUM 1001 ACCEPTED |E(fin-ini)=   -0.4090 |E(fin-sad)=   -0.4344 |npart=   18 |delr=    1.186 |evalf=   395 |
- That atom                    :               744
                 E-Eref  m_perp    ftot        fpar       fperp       eigen       delr  npart evalf   a1
                 ( eV )                     ( eV/Ang )              ( eV/Ang**2 )
   0  K=   0     0.0110  1 1      0.5893     -0.5266      0.2644      0.0000      0.054   0     2   0.00
   1  K=   1     0.1973  2 2      2.0867     -1.9790      0.6615      0.0000      0.287   0     5   0.00
   2  K=   2     0.5864  3 4      3.1487     -2.5694      1.8201      0.0000      0.579  12    10   0.00
SADDLE 1002    FAILED |ret  60132 |delta energy=    0.1470 |force_(tot,par,perp)=     0.0861    0.0179    0.0000 |eigenval=   0.1918 |npart=   20 |delr=    1.050 |evalf=   288 |
我想实现,从文件中提取这些信息
1、如果saddle xxxx 显示的是CONVERGED 那么记录下当前原子标号 ,例如:that atom : 272
并记录delta energy 例如:0.0255
2、这个文件特点是文件中saddle 显示的是CONVERGED,下面就是 RELAXATION过程
要记录E(fin-ini) ,E(fin-sad)
最后想要得到的是信息:272   0.0255   -0.4090   -0.4344
3、如果 saddle 显示的是  FAILED 那么就跳过不记录任何信息

不知道如何实现,是不是在原子标号处做上标记,在看saddle的情况,如果是CONVERGED ,就记录信息。如果不是再往下走,如何写成代码,还是没想好,希望高人指点!谢谢,我只有五个金币,全部奉上。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

karon1988

实习版主

优秀!!有木有!!!优秀!!有木有!!!优秀!!有木有!!!优秀!!有木有!!!

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
zhenyanwo: 金币+2, ★★★很有帮助, 谢谢,办法很好! 2015-07-23 09:44:56
你这需求, 用正则表达式吧

[ 发自小木虫客户端 ]
2楼2015-07-22 09:56:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jerkwin

新虫

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
zhenyanwo: 金币+3, ★★★很有帮助, 谢谢!办法很受用! 2015-07-23 09:45:33
awk吧, 或者perl
这些脚本语言都可以
3楼2015-07-22 21:07:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhenyanwo

兑换贵宾

优秀!!有木有!!!优秀!!有木有!!!优秀!!有木有!!!优秀!!有木有!!!

引用回帖:
3楼: Originally posted by jerkwin at 2015-07-22 21:07:12
awk吧, 或者perl
这些脚本语言都可以

谢谢,我先试试看。
4楼2015-07-23 08:24:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhenyanwo

超级版主

优秀!!有木有!!!优秀!!有木有!!!优秀!!有木有!!!优秀!!有木有!!!

引用回帖:
2楼: Originally posted by karon1988 at 2015-07-22 09:56:44
你这需求, 用正则表达式吧

谢谢!
5楼2015-07-23 08:42:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhenyanwo 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿西安交大材料学硕(英一数二)347,求调剂到高分子/材料相关专业 +5 zju51 2026-03-31 7/350 2026-04-01 00:47 by fmesaito
[考研] 070300化学专业279调剂 +10 哈哈哈^_^ 2026-03-31 10/500 2026-03-31 23:13 by liu823948201
[考研] 267求调剂 +9 uiybh 2026-03-31 9/450 2026-03-31 22:15 by 天倌赐福
[考研] 085600,321分求调剂 +7 大馋小子 2026-03-31 7/350 2026-03-31 18:02 by 无际的草原
[考研] 085601 329分调剂 +6 yzsa12 2026-03-31 6/300 2026-03-31 15:23 by yanflower7133
[考研] 学硕274求调剂 +17 Li李鱼 2026-03-26 17/850 2026-03-31 15:19 by 客尔美德
[考研] 求调剂 +4 研s. 2026-03-25 4/200 2026-03-31 10:05 by fmesaito
[考研] 一志愿食品科学与工程083200求调剂 +4 XQTJZ 2026-03-30 4/200 2026-03-31 04:10 by fmesaito
[考研] 0703化学321分求调剂 +10 三dd. 2026-03-30 11/550 2026-03-30 19:24 by markhwc
[考研] 297求调剂 +17 田洪有 2026-03-26 18/900 2026-03-30 18:32 by nothing投稿中
[考研] 292求调剂 +13 是妍子也是研子 2026-03-30 13/650 2026-03-30 18:01 by 小徐0109
[考研] 317分 一志愿南理工材料工程 本科湖工大 求调剂 +12 芋泥小铃铛 2026-03-28 12/600 2026-03-30 17:06 by wangjy2002
[考研] 298求调剂 +3 什么是胖头鱼 2026-03-30 5/250 2026-03-30 14:41 by 青海小西牛
[考研] 本科双非材料,跨考一志愿华电085801电气,283求调剂,任何专业都可以 +6 芝士雪baoo 2026-03-28 8/400 2026-03-29 08:16 by 松花缸1201
[考研] 081200-11408-276学硕求调剂 +6 崔wj 2026-03-26 6/300 2026-03-29 01:11 by hanserlol
[考研] 11408软件工程求调剂 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 21:50 by zhq0425
[考研] 266分求材料化工冶金矿业等专业的调剂 +4 哇呼哼呼哼 2026-03-26 4/200 2026-03-27 17:02 by zhyzzh
[考研] 315调剂 +4 0860求调剂 2026-03-26 5/250 2026-03-27 11:23 by wangjy2002
[考研] 调剂 +4 柚柚yoyo 2026-03-26 4/200 2026-03-26 20:43 by fmesaito
[考研] 环境专硕324分求调剂推荐 +5 轩小宁—— 2026-03-26 5/250 2026-03-26 12:05 by i_cooler
信息提示
请填处理意见