| 5 | 1/1 | 返回列表 |
| 查看: 590 | 回復(fù): 4 | |||
zhenyanwo鐵蟲(chóng) (初入文壇)
|
[求助]
大家?guī)兔纯聪襁@樣的文件如何提取想要的內(nèi)容。 已有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 | 我想實(shí)現(xiàn),從文件中提取這些信息 1、如果saddle xxxx 顯示的是CONVERGED 那么記錄下當(dāng)前原子標(biāo)號(hào) ,例如:that atom : 272 并記錄delta energy 例如:0.0255 2、這個(gè)文件特點(diǎn)是文件中saddle 顯示的是CONVERGED,下面就是 RELAXATION過(guò)程 要記錄E(fin-ini) ,E(fin-sad) 最后想要得到的是信息:272 0.0255 -0.4090 -0.4344 3、如果 saddle 顯示的是 FAILED 那么就跳過(guò)不記錄任何信息 不知道如何實(shí)現(xiàn),是不是在原子標(biāo)號(hào)處做上標(biāo)記,在看saddle的情況,如果是CONVERGED ,就記錄信息。如果不是再往下走,如何寫(xiě)成代碼,還是沒(méi)想好,希望高人指點(diǎn)!謝謝,我只有五個(gè)金幣,全部奉上。 |
專(zhuān)家顧問(wèn) (正式寫(xiě)手)
![]() |
專(zhuān)家經(jīng)驗(yàn): +14 |
鐵蟲(chóng) (初入文壇)
鐵蟲(chóng) (初入文壇)
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料080500調(diào)劑求收留 +5 | 一顆meteor 2026-03-13 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +3 | 申子申申 2026-03-19 | 7/350 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 288求調(diào)劑 +15 | 于海海海海 2026-03-19 | 15/750 |
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分,有科研經(jīng)歷,求調(diào)劑 +3 | 一瓶苯甲酸 2026-03-14 | 3/150 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
[考研] 317求調(diào)劑 +3 | 申子申申 2026-03-19 | 6/300 |
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
[考研] 材料專(zhuān)碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[基金申請(qǐng)] 被我言中:新模板不強(qiáng)調(diào)格式了,假專(zhuān)家開(kāi)始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 材料專(zhuān)碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專(zhuān)碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專(zhuān)利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語(yǔ) 2026-03-15 | 5/250 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|