应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » lammps计算结果xyz文件中每个数字代表的字符长短。
北京石油化工学院2026年研究生招生接收调剂公告
111/2返回列表12下一页
查看: 2167  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

herarysara

木虫 (小有名气)

[求助] lammps计算结果xyz文件中每个数字代表的字符长短。

lammps计算结果xyz文件，想提出文件中的数据，用c++的文件指针，可是不知道每个数据的字符长短，应该怎么确定，求赐教！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-12-04 19:55:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
zh1987hs: 金币+2, 谢谢 2012-12-05 08:13:34
1、lammps是开源软件，你为什么不去看看源程序呢？

2、你为什么不试试用lammps输出那么一小段xyz格式的结果，查看一下呢？
赞一下 回复此楼
2楼2012-12-04 22:12:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

herarysara

木虫 (小有名气)
引用回帖:
2楼: Originally posted by jiaoyixiong at 2012-12-04 22:12:40
1、lammps是开源软件，你为什么不去看看源程序呢？

2、你为什么不试试用lammps输出那么一小段xyz格式的结果，查看一下呢？

谢谢你的回复，其实我就是想知道他的输出文件每一行所占的字符数，单纯从文件上看到的字节数是不一样的。
赞一下 回复此楼
3楼2012-12-05 09:49:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

herarysara

木虫 (小有名气)
既是从xyz文件中，用fgets命令取走每个坐标如fgets(str,20,from)中间的20指每个坐标的字符长度，我想知道的就是这个数字应该怎么确定。xyz文件每行对应多少个字符。
赞一下 回复此楼
4楼2012-12-05 11:47:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiangtao9571

新虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
herarysara: 金币+5 2012-12-06 17:07:58
zh1987hs: 金币+3, 谢谢 2012-12-09 14:59:30
我试了一下得到的xyz中数据的格式如下：
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

这些都是数据，不要把它当成字符啊，这样不存在长短了啊，只要输出用浮点型，保证几位小数就行，我写的小程序供你参考：
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);
                }
}

这个很详细了哦，仔细看看哦，估计你能看懂，如果不行的话要去补习C++了。
赞一下 回复此楼
5楼2012-12-06 16:38:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

herarysara

木虫 (小有名气)
确实没咋看懂，我是用文件指针形式逐个读取xyz文件中的数据，因为涉及到轨迹文件几十万帧，感觉这样实现不了。不过还是非常感谢你的回复。
赞一下 回复此楼
6楼2012-12-06 17:07:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiangtao9571

新虫 (正式写手)
★ ★
zh1987hs: 金币+2, 谢谢 2012-12-09 14:59:46
引用回帖:
6楼: Originally posted by herarysara at 2012-12-06 17:07:34
确实没咋看懂，我是用文件指针形式逐个读取xyz文件中的数据，因为涉及到轨迹文件几十万帧，感觉这样实现不了。不过还是非常感谢你的回复。

你看不出来我也是用文件指针形式逐个读取xyz文件中的数据吗？？？  涉及几十万帧  加一个循环即可：
如下：step = timestep/10000；
        step=0;
        while(step<20)
        {
        readdata();
        mol=atom_tol/atom_per;
        writefile();
        }
赞一下 回复此楼
7楼2012-12-07 09:56:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

herarysara

木虫 (小有名气)
引用回帖:
7楼: Originally posted by jiangtao9571 at 2012-12-07 09:56:01
你看不出来我也是用文件指针形式逐个读取xyz文件中的数据吗？？？  涉及几十万帧  加一个循环即可：
如下：step = timestep/10000；
        step=0;
        while(step<20)
        {
        readdata();
        mol=atom_tol/ato ...

你上面的程序，是有指定选取80个字符吗？如果是，我就是想知道读取这个字符长度该怎样确定。
赞一下 回复此楼
8楼2012-12-07 16:54:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★
zh1987hs: 金币+2, 谢谢 2012-12-09 15:00:01
不知道你的问题解决了没？
xyz文件内，主要就是存储的是原子的坐标信息吧？
如果把xyz轨迹文件导入到VMD，使用命令：
set sel [atomselect top "the atoms you need"]
$sel get {x y z}
就可以读取出你所需要的原子的坐标了
如果只需x坐标，就写成 $sel get {x }
即可
这比自己写那么长的程序简单多了，楼主不妨试试
赞一下 回复此楼
9楼2012-12-08 10:25:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

herarysara

木虫 (小有名气)
引用回帖:
9楼: Originally posted by jiaoyixiong at 2012-12-08 10:25:32
不知道你的问题解决了没？
xyz文件内，主要就是存储的是原子的坐标信息吧？
如果把xyz轨迹文件导入到VMD，使用命令：
set sel
$sel get {x y z}
就可以读取出你所需要的原子的坐标了
如果只需x坐标，就写成  ...

谢谢你的回复，我不是要单纯得到某个原子的坐标数字，我是想以字符串的的形式读取。
赞一下 回复此楼
10楼2012-12-09 10:56:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 herarysara 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程调剂 +7 小刘同学吖吖 2026-04-06 8/400 2026-04-08 17:18 by 努力向上，不断�
[考研] 一志愿211，化学学硕，310分，本科重点双非，求调剂 +13 努力奋斗112 2026-04-06 14/700 2026-04-08 15:07 by screening
[考研] 315求调剂 +17 欣喜777 2026-04-04 18/900 2026-04-08 13:54 by hangsimei
[考研] 280求调剂 +9 李rien 2026-04-04 9/450 2026-04-08 13:20 by lijunpoly
[考研] 304求调剂 +16 c297914 2026-04-05 17/850 2026-04-08 13:00 by grayjzr
[考研] 化工求调剂！ +26 RichLi_ 2026-04-06 26/1300 2026-04-08 10:28 by tj2m
[考研] 307求调剂 +14 超级伊昂大王 2026-04-06 14/700 2026-04-08 07:03 by 无际的草原
[考研] 一志愿211电子信息347求调剂 +3 554916 2026-04-03 3/150 2026-04-07 23:22 by 如若时光倒流
[考博] 博士申请 +3 IQwQl 2026-04-05 3/150 2026-04-07 20:31 by greychen00
[考研] 286求调剂 +20 Faune 2026-04-06 20/1000 2026-04-07 11:33 by 诗与自由
[考研] 085600，321分求调剂 +14 大馋小子 2026-04-04 15/750 2026-04-06 22:58 by qlm5820
[考研] 284求调剂 +7 徐同学_001 2026-04-04 13/650 2026-04-05 17:19 by yulian1987
[考研] 考研生物学考A区211，初试322，科目生化和生物综合，求调剂 +6 。。。54 2026-04-03 6/300 2026-04-05 14:54 by JOKER0401
[考研] 材料调剂 +12 一样YWY 2026-04-04 12/600 2026-04-05 08:24 by 544594351
[考研] 286求调剂 +3 草木不言 2026-04-04 3/150 2026-04-04 22:40 by lbsjt
[考研] 求调剂 +4 晟功? 2026-04-03 4/200 2026-04-04 21:58 by hemengdong
[考研] 321求调剂 +13 认真求上学 2026-04-02 13/650 2026-04-04 18:23 by macy2011
[考研] 297求调剂 +11 ljy20040718！ 2026-04-03 13/650 2026-04-04 09:23 by 来看流星雨10
[考研] 085601一志愿北理325分求调剂 +6 找调剂，， 2026-04-02 6/300 2026-04-03 22:20 by 柟�?
[考研] 293求调剂 +4 珂珂乐 2026-04-02 4/200 2026-04-02 20:10 by 6781022
信息提示
关闭
请填处理意见
关闭