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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » lammps计算结果xyz文件中每个数字代表的字符长短。
111/2返回列表12下一页
查看: 2143  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 梁成伟老师课题组欢迎你的加入 +9 一鸭鸭哟 2026-03-14 11/550 2026-03-19 17:22 by ！本暗一次！
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 6/300 2026-03-19 17:18 by fei626-918
[考研] 本人考085602 化学工程 专硕 +17 不知道叫什么！ 2026-03-15 19/950 2026-03-19 15:06 by 尽舜尧1
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 294求调剂材料与化工专硕 +11 陌の森林 2026-03-18 11/550 2026-03-19 13:22 by houyaoxu
[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 332求调剂 +3 ydfyh 2026-03-17 3/150 2026-03-19 10:14 by 功夫疯狂
[考研] 0817调剂 +3 没有答案_ 2026-03-14 3/150 2026-03-19 09:51 by Xu de nuo
[考研] 354求调剂 +4 Tyoumou 2026-03-18 7/350 2026-03-18 21:45 by Tyoumou
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 334求调剂 +3 志存高远意在机� 2026-03-16 3/150 2026-03-18 08:34 by lm4875102
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 机械专硕325，寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生，2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
信息提示
关闭
请填处理意见
关闭