| 查看: 2437 | 回复: 10 | ||
herarysara木虫 (小有名气)
|
[求助]
lammps计算结果xyz文件中每个数字代表的字符长短。
|
| lammps计算结果xyz文件，想提出文件中的数据，用c++的文件指针，可是不知道每个数据的字符长短，应该怎么确定，求赐教！ |
回复此楼» 猜你喜欢
基金啊基金
已经有3人回复
-
科研孤儿太难了
已经有15人回复
-
让我中一个面上吧！
已经有10人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
-
重要消息，中午系统在维护
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Lammps计算所用的势函数与晶体结构有关吗？
已经有4人回复
- lammps 计算径向分布函数并作图-求助
已经有12人回复
- Lammps计算SiO2熔点，原子丢失错误
已经有10人回复
- lammps 计算粘度 lj单位 转化为 metal单位 出错
已经有7人回复
- 【转】lammps使用（Michael博客）
已经有316人回复
- lammps计算碳纳米管热导率 建模问题和热导率计算交流
已经有5人回复
- lammps中速度自相关函数的计算
已经有11人回复
- VMD建模用于Lammps计算后，结果在MS和VMD中显示不一样，怎么回事呢？
已经有0人回复
- lammps计算晶格热导
已经有2人回复
- 求助LAMMPS两次运算结果不一致问题！
已经有11人回复
- 请教大神如何在64bit的win7系统下正确安装Lammps？计算机是4核的处理器
已经有14人回复
- lammps计算msd的结果如何输出？
已经有2人回复
- lammps计算粘度 GK法和rNEMD法
已经有4人回复
- LAMMPS计算无反应，nan
已经有4人回复
- LAMMPS运行结果出现nan
已经有20人回复
- LAMMPS并行计算的问题（cpu——time关系）
已经有17人回复
2楼2012-12-04 22:12:40
herarysara
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 3303.7
- 散金: 23
- 红花: 2
- 帖子: 210
- 在线: 89.7小时
- 虫号: 972766
- 注册: 2010-03-16
- 性别: GG
- 专业: 化学反应工程
3楼2012-12-05 09:49:33
herarysara
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 3303.7
- 散金: 23
- 红花: 2
- 帖子: 210
- 在线: 89.7小时
- 虫号: 972766
- 注册: 2010-03-16
- 性别: GG
- 专业: 化学反应工程
4楼2012-12-05 11:47:01
jiangtao9571
新虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 238.2
- 散金: 1046
- 红花: 1
- 帖子: 454
- 在线: 236.2小时
- 虫号: 706996
- 注册: 2009-02-23
- 专业: 高分子材料结构与性能
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
herarysara: 金币+5 2012-12-06 17:07:58
zh1987hs: 金币+3, 谢谢 2012-12-09 14:59:30
感谢参与，应助指数 +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
herarysara
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 3303.7
- 散金: 23
- 红花: 2
- 帖子: 210
- 在线: 89.7小时
- 虫号: 972766
- 注册: 2010-03-16
- 性别: GG
- 专业: 化学反应工程
6楼2012-12-06 17:07:34
jiangtao9571
新虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 238.2
- 散金: 1046
- 红花: 1
- 帖子: 454
- 在线: 236.2小时
- 虫号: 706996
- 注册: 2009-02-23
- 专业: 高分子材料结构与性能
涉及几十万帧 加一个循环即可：7楼2012-12-07 09:56:01
herarysara
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 3303.7
- 散金: 23
- 红花: 2
- 帖子: 210
- 在线: 89.7小时
- 虫号: 972766
- 注册: 2010-03-16
- 性别: GG
- 专业: 化学反应工程
8楼2012-12-07 16:54:34
9楼2012-12-08 10:25:32
herarysara
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 3303.7
- 散金: 23
- 红花: 2
- 帖子: 210
- 在线: 89.7小时
- 虫号: 972766
- 注册: 2010-03-16
- 性别: GG
- 专业: 化学反应工程
10楼2012-12-09 10:56:50
40