 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 8 | 1/1 | 返回列表 |
| 查看: 4389 | 回復(fù): 7 | |||
zhuhongaaa金蟲 (正式寫手)
|
[交流]
【轉(zhuǎn)帖】【ZZ】使用xdrfile讀取gromacs的trr和xtc文件 已有7人參與
|
|
xdrfile是gmx開發(fā)的一個用來讀取trr文件和xtc文件的庫，可以在gmx的主頁上下載安裝。安裝不需要一共gmx的安裝路徑。 安裝后生成一個trr2xtc可執(zhí)行文件�？梢院唵蔚陌裻rr文件轉(zhuǎn)化問xtc文件來減小體積，節(jié)約空間。 更為主要的是我們可以利用xdrfile中的庫自行編寫程序讀取gmx的trr文件和xtc文件，用于軌跡分析，相當(dāng)方便。 xdrfile/include/xdrfile/目錄下面有三個用C寫的頭文件，可以include到我們自己的C程序和C++程序中。 下面是一個讀取trr文件的C程序示例： /* file name is trr.c */ #include #include #include int natoms,step,natom,read_return; float time,lambda; matrix box; rvec *x,*v,*f; XDRFILE *trr; main () { trr=xdrfile_open (“traj.trr”,”r”); read_return=read_trr_natoms (“traj.trr”,&natoms); x = calloc(natoms, sizeof(x[0])); while (1) { read_return=read_trr (trr,natoms,&step,&time,box,lambda,x,v,f); if (read_return!=0) { break; } for (natom=1;natom<=natoms;natom++) { printf (“%d %f %d %f %f %f\n”,step,time,natom,x[natom-1][0],x[natom-1][1],x[natom-1][2]); } } xdrfile_close (trr); } ~ ~ ~ 可以用gcc 直接編譯運行。編譯的時候注意添加頭文件的路徑。 下面是用C++寫的讀取xtc文件的一個示例： /* file name is xtc.cpp */ #include extern “C” { #include } using namespace std; int main (int argc,char * argv[]) { char * filename; int natoms,step; filename=”traj.xtc”; float time; float p; matrix box; rvec *x; XDRFILE *xtc; xtc=xdrfile_open(filename,”r”); int read_return=read_xtc_natoms(filename,&natoms); x=(rvec * )calloc(natoms,sizeof(x[0])); while(1) { read_return=read_xtc(xtc,natoms,&step,&time,box,x,&p); if(read_return!=0) { break; } for(int natom=0;natom cout< } xdrfile_close(xtc); return 0; } ~ 可以用g++ 直接編譯運行。編譯的時候注意添加頭文件的路徑。 http://zhuhong.freetzi.com/wordpress/?p=57 有問題歡迎留言或mail。 轉(zhuǎn)載請注明。 |
回復(fù)此樓» 猜你喜歡
324分 085600材料化工求調(diào)劑
已經(jīng)有3人回復(fù)
-
317求調(diào)劑
已經(jīng)有6人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有11人回復(fù)
-
307求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有3人回復(fù)
-
求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分
已經(jīng)有3人回復(fù)
-
化學(xué)求調(diào)劑
已經(jīng)有4人回復(fù)
-
招收調(diào)劑碩士
已經(jīng)有9人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有8人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 請教一下，如何確定具體分子式，最好有詳細步驟，還有生成CIF文件的問題
已經(jīng)有7人回復(fù)
- ppt動畫（求職自我介紹、新的征途、華科）+其他ppt部分教程+ppt演示、使用操作大全
已經(jīng)有100人回復(fù)
- 如何用MS的納米管坐標(biāo)轉(zhuǎn)到gromacs可識別的文件
已經(jīng)有7人回復(fù)
- Guidelines for Design of Wind Turbines(Low Resolution) 單個文件
已經(jīng)有11人回復(fù)
- latex模板里的文件是什么意思，怎么用？
已經(jīng)有5人回復(fù)
- 關(guān)于matlab中循環(huán)讀取多個excel文件的問題？
已經(jīng)有3人回復(fù)
- linux下面使用castep跑出來結(jié)果文件中沒有xsd文件
已經(jīng)有10人回復(fù)
- gromacs模擬
已經(jīng)有4人回復(fù)
- GROMACS 小分子top文件
已經(jīng)有9人回復(fù)
- gromacs模擬離子液體，建立gro文件
已經(jīng)有5人回復(fù)
- GROMACS ： 把兩個*.gro文件合并
已經(jīng)有4人回復(fù)
- 【求助】如何使用gromacs中的constraint
已經(jīng)有8人回復(fù)
- 【求助】重金求助：gromacs中g(shù)_rms命令中怎樣應(yīng)用index.ndx文件
已經(jīng)有8人回復(fù)
- 【求助】求助Gromacs中關(guān)于g_hbond的out文件的問題
已經(jīng)有10人回復(fù)
piao666
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 10
- 帖子: 1
- 在線: 21分鐘
- 蟲號: 10452419
- 注冊: 2018-09-26
- 專業(yè): 其他無機非金屬材料
2樓2018-09-26 15:44:33
小城小暖
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 137.1
- 散金: 5
- 帖子: 49
- 在線: 4.4小時
- 蟲號: 7918487
- 注冊: 2018-02-05
- 性別: GG
- 專業(yè): 分離過程
3樓2019-12-01 20:20:47
★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
|
本帖內(nèi)容被屏蔽 |
4樓2019-12-31 16:21:25
★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
|
本帖內(nèi)容被屏蔽 |
5樓2020-01-02 17:25:31
★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
|
本帖內(nèi)容被屏蔽 |
6樓2020-01-04 11:20:13
|
本帖內(nèi)容被屏蔽 |
7樓2020-01-06 16:59:44
★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
|
本帖內(nèi)容被屏蔽 |
8樓2022-10-25 01:26:52
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zhuhongaaa 的主題更新
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 307求調(diào)劑 +5 | 冷笙123 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +4 | 吃吃吃才有意義 2026-03-19 | 4/200 |
|
|
|
[教師之家] 焦慮 +9 | 水冰月月野兔 2026-03-13 | 13/650 |
|
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +6 | 薛云鵬 2026-03-13 | 6/300 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機器學(xué)習(xí) 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|