 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 8 | 1/1 | 返回列表 |
| 查看: 4391 | 回復(fù): 7 | |||
zhuhongaaa金蟲 (正式寫手)
|
[交流]
【轉(zhuǎn)帖】【ZZ】使用xdrfile讀取gromacs的trr和xtc文件 已有7人參與
|
|
xdrfile是gmx開發(fā)的一個用來讀取trr文件和xtc文件的庫，可以在gmx的主頁上下載安裝。安裝不需要一共gmx的安裝路徑。 安裝后生成一個trr2xtc可執(zhí)行文件。可以簡單的把trr文件轉(zhuǎn)化問xtc文件來減小體積，節(jié)約空間。 更為主要的是我們可以利用xdrfile中的庫自行編寫程序讀取gmx的trr文件和xtc文件，用于軌跡分析，相當方便。 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ù)此樓» 猜你喜歡
307求調(diào)劑
已經(jīng)有9人回復(fù)
-
288求調(diào)劑
已經(jīng)有15人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有5人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有14人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有18人回復(fù)
-
265求調(diào)劑
已經(jīng)有3人回復(fù)
-
296求調(diào)劑
已經(jīng)有6人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有13人回復(fù)
-
085600材料與化工 求調(diào)劑
已經(jīng)有15人回復(fù)
-
生物學(xué)調(diào)劑招人�。�！
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 請教一下，如何確定具體分子式，最好有詳細步驟，還有生成CIF文件的問題
已經(jīng)有7人回復(fù)
- ppt動畫（求職自我介紹、新的征途、華科）+其他ppt部分教程+ppt演示、使用操作大全
已經(jīng)有100人回復(fù)
- 如何用MS的納米管坐標轉(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ā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +5 | Doleres 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 0703化學(xué)調(diào)劑 ，六級已過，有科研經(jīng)歷 +12 | 曦熙兮 2026-03-15 | 12/600 |
|
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +3 | Mqqqqqq 2026-03-19 | 3/150 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 化學(xué)求調(diào)劑 +3 | 臨澤境llllll 2026-03-17 | 4/200 |
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +5 | z1z2z3879 2026-03-18 | 5/250 |
|
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +7 | 手機用戶 2026-03-17 | 7/350 |
|
|
|
[考研]
 0854，計算機類招收調(diào)劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風就是雨； 2026-03-16 | 3/150 |
|
|
|
[基金申請]
今年的國基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|