版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

下一頁(yè)

herarysara

木蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 3337.7
散金: 6
紅花: 2
帖子: 182
在線: 83.8小時(shí)
蟲號(hào): 972766
注冊(cè): 2010-03-16
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

[求助] lammps計(jì)算結(jié)果xyz文件中每個(gè)數(shù)字代表的字符長(zhǎng)短。

lammps計(jì)算結(jié)果xyz文件，想提出文件中的數(shù)據(jù)，用c++的文件指針，可是不知道每個(gè)數(shù)據(jù)的字符長(zhǎng)短，應(yīng)該怎么確定，求賜教！

回復(fù)此樓

» 猜你喜歡

26考研求調(diào)劑已經(jīng)有6人回復(fù)
求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以已經(jīng)有5人回復(fù)
085600材料與化工求調(diào)劑已經(jīng)有3人回復(fù)
311求調(diào)劑已經(jīng)有7人回復(fù)
303求調(diào)劑已經(jīng)有3人回復(fù)
有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～已經(jīng)有7人回復(fù)
梁成偉老師課題組歡迎你的加入已經(jīng)有10人回復(fù)
一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) 已經(jīng)有8人回復(fù)
070300化學(xué)319求調(diào)劑已經(jīng)有3人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有4人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

Lammps計(jì)算所用的勢(shì)函數(shù)與晶體結(jié)構(gòu)有關(guān)嗎？已經(jīng)有4人回復(fù)
lammps 計(jì)算徑向分布函數(shù)并作圖-求助已經(jīng)有12人回復(fù)
Lammps計(jì)算SiO2熔點(diǎn)，原子丟失錯(cuò)誤已經(jīng)有10人回復(fù)
lammps 計(jì)算粘度 lj單位轉(zhuǎn)化為 metal單位出錯(cuò) 已經(jīng)有7人回復(fù)
【轉(zhuǎn)】lammps使用（Michael博客）已經(jīng)有315人回復(fù)
lammps計(jì)算碳納米管熱導(dǎo)率建模問(wèn)題和熱導(dǎo)率計(jì)算交流已經(jīng)有5人回復(fù)
lammps中速度自相關(guān)函數(shù)的計(jì)算已經(jīng)有11人回復(fù)
VMD建模用于Lammps計(jì)算后，結(jié)果在MS和VMD中顯示不一樣，怎么回事呢？已經(jīng)有0人回復(fù)
lammps計(jì)算晶格熱導(dǎo) 已經(jīng)有2人回復(fù)
求助LAMMPS兩次運(yùn)算結(jié)果不一致問(wèn)題！已經(jīng)有11人回復(fù)
請(qǐng)教大神如何在64bit的win7系統(tǒng)下正確安裝Lammps？計(jì)算機(jī)是4核的處理器已經(jīng)有14人回復(fù)
lammps計(jì)算msd的結(jié)果如何輸出？已經(jīng)有2人回復(fù)
lammps計(jì)算粘度 GK法和rNEMD法已經(jīng)有4人回復(fù)
LAMMPS計(jì)算無(wú)反應(yīng)，nan 已經(jīng)有4人回復(fù)
LAMMPS運(yùn)行結(jié)果出現(xiàn)nan 已經(jīng)有20人回復(fù)
LAMMPS并行計(jì)算的問(wèn)題（cpu——time關(guān)系）已經(jīng)有17人回復(fù)

1樓 2012-12-04 19:55:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jiaoyixiong

榮譽(yù)版主 (職業(yè)作家)

模擬EPI: 6
應(yīng)助: 487 (碩士)
貴賓: 0.501
金幣: 14945
散金: 6107
紅花: 74
帖子: 3750
在線: 1129.1小時(shí)
蟲號(hào): 1108016
注冊(cè): 2010-09-26
性別: GG
專業(yè): 流體力學(xué)
管轄: 分子模擬

【答案】應(yīng)助回帖

★ ★
感謝參與，應(yīng)助指數(shù) +1
zh1987hs: 金幣+2, 謝謝 2012-12-05 08:13:34

1、lammps是開源軟件，你為什么不去看看源程序呢？

2、你為什么不試試用lammps輸出那么一小段xyz格式的結(jié)果，查看一下呢？

贊一下

回復(fù)此樓

2樓2012-12-04 22:12:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

herarysara

木蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 3337.7
散金: 6
紅花: 2
帖子: 182
在線: 83.8小時(shí)
蟲號(hào): 972766
注冊(cè): 2010-03-16
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

2樓: Originally posted by jiaoyixiong at 2012-12-04 22:12:40
1、lammps是開源軟件，你為什么不去看看源程序呢？

2、你為什么不試試用lammps輸出那么一小段xyz格式的結(jié)果，查看一下呢？

謝謝你的回復(fù)，其實(shí)我就是想知道他的輸出文件每一行所占的字符數(shù)，單純從文件上看到的字節(jié)數(shù)是不一樣的。

贊一下

回復(fù)此樓

3樓2012-12-05 09:49:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

herarysara

木蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 3337.7
散金: 6
紅花: 2
帖子: 182
在線: 83.8小時(shí)
蟲號(hào): 972766
注冊(cè): 2010-03-16
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

既是從xyz文件中，用fgets命令取走每個(gè)坐標(biāo)如fgets(str,20,from)中間的20指每個(gè)坐標(biāo)的字符長(zhǎng)度，我想知道的就是這個(gè)數(shù)字應(yīng)該怎么確定。xyz文件每行對(duì)應(yīng)多少個(gè)字符。

贊一下

回復(fù)此樓

4樓2012-12-05 11:47:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jiangtao9571

新蟲 (正式寫手)

應(yīng)助: 4 (幼兒園)
金幣: 232.4
散金: 1046
紅花: 1
帖子: 452
在線: 236.2小時(shí)
蟲號(hào): 706996
注冊(cè): 2009-02-23
專業(yè): 高分子材料結(jié)構(gòu)與性能

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
herarysara: 金幣+5 2012-12-06 17:07:58
zh1987hs: 金幣+3, 謝謝 2012-12-09 14:59:30

我試了一下得到的xyz中數(shù)據(jù)的格式如下：
ITEM: TIMESTEP
0
ITEM: NUMBER OF ATOMS
10000
ITEM: BOX BOUNDS
0 22.66
0 22.66
0 22.66
ITEM: ATOMS id type xs ys zs
5943 1 0.0583092 0.053962 0.061988
412 1 0.0225291 0.0612129 0.058574
2582 1 0.0497735 0.0260711 0.0323743

這些都是數(shù)據(jù)，不要把它當(dāng)成字符啊，這樣不存在長(zhǎng)短了啊，只要輸出用浮點(diǎn)型，保證幾位小數(shù)就行，我寫的小程序供你參考：
struct bankdata
{
int num1;
int num2;
double x;
double y;
double z;
}data[200000],datastation;

struct boxsize_infomation
{

double a;
double b;
}boxsize[3];

void readdata()
{
int i=0;
char line[80];
printf("Reading atoms.xyz.........\n"

;
fgets(line,80,fp1);
fgets(line,80,fp1); sscanf(line,"%d\n",×tep);
fgets(line,80,fp1);
fgets(line,80,fp1); sscanf(line,"%d\n",&atom_tol);
fgets(line,80,fp1);
fgets(line,80,fp1); sscanf(line,"%lf %lf\n",&boxsize[0].a,&boxsize[0].b);
fgets(line,80,fp1); sscanf(line,"%lf %lf\n",&boxsize[1].a,&boxsize[1].b);
fgets(line,80,fp1); sscanf(line,"%lf %lf\n",&boxsize[2].a,&boxsize[2].b);
fgets(line,80,fp1);

printf("step=%d,timestep=%d\n",step,timestep);

while(i {
fgets(line,80,fp1);
sscanf(line,"%d %d %lf %lf %lf\n",&datastation.num1,&datastation.num2,&datastation.x,&datastation.y,&datastation.z);

data.num1=datastation.num1;
data.num2=datastation.num2;
data.x=datastation.x;
data.y=datastation.y;
data.z=datastation.z;
i++;
}
}

void writefile()
{
fprintf(fp,"ITEM: TIMESTEP\n";
fprintf(fp,"%d\n",timestep);
fprintf(fp,"ITEM: NUMBER OF ATOMS\n";
fprintf(fp,"%d\n",atom_tol);
fprintf(fp,"ITEM: BOX BOUNDS\n";
fprintf(fp,"%g %g\n",boxsize[0].a,boxsize[0].b);
fprintf(fp,"%g %g\n",boxsize[1].a,boxsize[1].b);
fprintf(fp,"%g %g\n",boxsize[2].a,boxsize[2].b);
fprintf(fp,"ITEM: ATOMS id type xs ys zs\n";

for (i=0; i {
fprintf(fp,"%d %d %lf %lf %lf \n",data.num1,data.num2,data.x,data.y,data.z);
}
}

這個(gè)很詳細(xì)了哦，仔細(xì)看看哦，估計(jì)你能看懂，如果不行的話要去補(bǔ)習(xí)C++了。

贊一下

回復(fù)此樓

5樓2012-12-06 16:38:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

herarysara

木蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 3337.7
散金: 6
紅花: 2
帖子: 182
在線: 83.8小時(shí)
蟲號(hào): 972766
注冊(cè): 2010-03-16
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

確實(shí)沒(méi)咋看懂，我是用文件指針形式逐個(gè)讀取xyz文件中的數(shù)據(jù)，因?yàn)樯婕暗杰壽E文件幾十萬(wàn)幀，感覺這樣實(shí)現(xiàn)不了。不過(guò)還是非常感謝你的回復(fù)。

贊一下

回復(fù)此樓

6樓2012-12-06 17:07:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jiangtao9571

新蟲 (正式寫手)

應(yīng)助: 4 (幼兒園)
金幣: 232.4
散金: 1046
紅花: 1
帖子: 452
在線: 236.2小時(shí)
蟲號(hào): 706996
注冊(cè): 2009-02-23
專業(yè): 高分子材料結(jié)構(gòu)與性能

★ ★
zh1987hs: 金幣+2, 謝謝 2012-12-09 14:59:46

引用回帖:

6樓: Originally posted by herarysara at 2012-12-06 17:07:34
確實(shí)沒(méi)咋看懂，我是用文件指針形式逐個(gè)讀取xyz文件中的數(shù)據(jù)，因?yàn)樯婕暗杰壽E文件幾十萬(wàn)幀，感覺這樣實(shí)現(xiàn)不了。不過(guò)還是非常感謝你的回復(fù)。

你看不出來(lái)我也是用文件指針形式逐個(gè)讀取xyz文件中的數(shù)據(jù)嗎？？？

涉及幾十萬(wàn)幀加一個(gè)循環(huán)即可：
如下：step = timestep/10000；
step=0;
while(step<20)
{
readdata();
mol=atom_tol/atom_per;
writefile();
}

贊一下

回復(fù)此樓

7樓2012-12-07 09:56:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

herarysara

木蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 3337.7
散金: 6
紅花: 2
帖子: 182
在線: 83.8小時(shí)
蟲號(hào): 972766
注冊(cè): 2010-03-16
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

7樓: Originally posted by jiangtao9571 at 2012-12-07 09:56:01
你看不出來(lái)我也是用文件指針形式逐個(gè)讀取xyz文件中的數(shù)據(jù)嗎？？？

涉及幾十萬(wàn)幀加一個(gè)循環(huán)即可：
如下：step = timestep/10000；
step=0;
while(step<20)
{
readdata();
mol=atom_tol/ato ...

你上面的程序，是有指定選取80個(gè)字符嗎？如果是，我就是想知道讀取這個(gè)字符長(zhǎng)度該怎樣確定。

贊一下

回復(fù)此樓

8樓2012-12-07 16:54:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jiaoyixiong

榮譽(yù)版主 (職業(yè)作家)

模擬EPI: 6
應(yīng)助: 487 (碩士)
貴賓: 0.501
金幣: 14945
散金: 6107
紅花: 74
帖子: 3750
在線: 1129.1小時(shí)
蟲號(hào): 1108016
注冊(cè): 2010-09-26
性別: GG
專業(yè): 流體力學(xué)
管轄: 分子模擬

【答案】應(yīng)助回帖

★ ★
zh1987hs: 金幣+2, 謝謝 2012-12-09 15:00:01

不知道你的問(wèn)題解決了沒(méi)？
xyz文件內(nèi)，主要就是存儲(chǔ)的是原子的坐標(biāo)信息吧？
如果把xyz軌跡文件導(dǎo)入到VMD，使用命令：
set sel [atomselect top "the atoms you need"]
$sel get {x y z}
就可以讀取出你所需要的原子的坐標(biāo)了
如果只需x坐標(biāo)，就寫成 $sel get {x }
即可
這比自己寫那么長(zhǎng)的程序簡(jiǎn)單多了，樓主不妨試試

贊一下

回復(fù)此樓

9樓2012-12-08 10:25:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

herarysara

木蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 3337.7
散金: 6
紅花: 2
帖子: 182
在線: 83.8小時(shí)
蟲號(hào): 972766
注冊(cè): 2010-03-16
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

9樓: Originally posted by jiaoyixiong at 2012-12-08 10:25:32
不知道你的問(wèn)題解決了沒(méi)？
xyz文件內(nèi)，主要就是存儲(chǔ)的是原子的坐標(biāo)信息吧？
如果把xyz軌跡文件導(dǎo)入到VMD，使用命令：
set sel
$sel get {x y z}
就可以讀取出你所需要的原子的坐標(biāo)了
如果只需x坐標(biāo)，就寫成 ...

謝謝你的回復(fù)，我不是要單純得到某個(gè)原子的坐標(biāo)數(shù)字，我是想以字符串的的形式讀取。

贊一下

回復(fù)此樓

10樓2012-12-09 10:56:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 herarysara 的主題更新

返回列表

下一頁(yè)

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 085600材料與化工求調(diào)劑 +3	緒幸與子 2026-03-17	3/150	2026-03-17 15:29 by 霧散后相遇lc
[考研] 301求調(diào)劑 +3	yy要上岸呀 2026-03-17	3/150	2026-03-17 14:04 by laoshidan
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +3	share_joy 2026-03-16	3/150	2026-03-17 13:41 by houyaoxu
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 286求調(diào)劑 +3	lemonzzn 2026-03-16	5/250	2026-03-16 20:43 by lemonzzn
[考研] 化學(xué)調(diào)劑0703 +8	啊我我的 2026-03-11	8/400	2026-03-16 17:23 by 我的船我的海
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 265求調(diào)劑 +4	威化餅07 2026-03-12	4/200	2026-03-14 17:23 by userper
[考研] 255求調(diào)劑 +3	李嘉慧， 2026-03-12	4/200	2026-03-14 16:58 by 有只貍奴
[考研] 308 085701 四六級(jí)已過(guò)求調(diào)劑 +7	溫喬喬喬喬 2026-03-12	14/700	2026-03-14 10:49 by JourneyLucky
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6	邱gl 2026-03-12	7/350	2026-03-13 23:24 by 邱gl
[考研] 材料371求調(diào)劑 +9	鱷魚? 2026-03-11	11/550	2026-03-13 22:53 by JourneyLucky
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8	dream…… 2026-03-12	8/400	2026-03-13 22:17 by 星空星月
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月
[考研] 290求調(diào)劑 +9	ADT 2026-03-11	9/450	2026-03-13 21:55 by JourneyLucky
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8	努力奮斗112 2026-03-12	9/450	2026-03-13 15:41 by JourneyLucky
[考研] 290求調(diào)劑 +7	ADT 2026-03-12	7/350	2026-03-13 15:17 by JourneyLucky
[考研] 0856化學(xué)工程280分求調(diào)劑 +4	shenzxsn 2026-03-11	4/200	2026-03-13 11:55 by ymwdoctor
[考研] 296求調(diào)劑 +3	大口吃飯身體健 2026-03-13	3/150	2026-03-13 10:31 by 學(xué)員8dgXkO

24小時(shí)熱門版塊排行榜

[求助] lammps計(jì)算結(jié)果xyz文件中每個(gè)數(shù)字代表的字符長(zhǎng)短。

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] lammps計(jì)算結(jié)果xyz文件中每個(gè)數(shù)字代表的字符長(zhǎng)短。

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: