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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » lammps计算结果xyz文件中每个数字代表的字符长短。
51/1返回列表
查看: 1994  |  回复: 10
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

herarysara

木虫 (小有名气)

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

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

» 猜你喜欢

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

herarysara

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

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的回帖

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的回帖
查看全部 11 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭