| 查看: 578 | 回复: 3 | |||
[交流]
【求助】建立体心晶胞结构代码，运行不出结果！已有2人参与
|
|
用VC++6.0编的一个建立体心立方晶胞结构的代码，编译成功了，但是运行不出结果，找不到错误所在，请高手帮忙一下啊！十分感谢！ void CBCC_CrystalDlg::OnCalculate() { // TODO: Add your control notification handler code here FILE * fileFe_xyz = fopen("Fe_xyz.txt", "w+"  ;fprintf(fileFe_xyz, "i xbas ybas zbas\n" ;a=2.8664; //Fe的晶格常数 xunit[0]=0.0; //以下为单胞原子坐标 yunit[0]=0.0; zunit[0]=0.0; xunit[1]=0.50; yunit[1]=0.50; zunit[1]=0.50; ncellx=M; //ncellx,ncelly,ncellz是在xyz方向的重复次数，其中M任意整数 ncelly=M; ncellz=M; iatoms=0; natoms=2*ncellx*ncelly*ncellz; for(int i=0;i xbas=0; ybas=0; zbas=0; } for(ix=1;i<=ncellx;ix++) { for(iy=1;iy<=ncelly;iy++) { for(iz=1;iz<=ncellz;iz++) { for(ibas=0;ibas<2;ibas++) { xbas[iatoms]=(float(ix-1)+xunit[ibas])*a; ybas[iatoms]=(float(iy-1)+yunit[ibas])*a; zbas[iatoms]=(float(iz-1)+zunit[ibas])*a; iatoms=iatoms+1; } } } } for(i=0;i fprintf(fileFe_xyz, "%6d %10.5f %10.5f %10.5f\n", i,xbas,ybas,zbas); } fclose(fileFe_xyz); } 这是头文件中定义的变量及函数： #define M 5 public: int i; int ix; int iy; int iz; int ncellx; int ncelly; int ncellz; int iatoms; int natoms; int ibas; double xunit[2]; double yunit[2]; double zunit[2]; double a; double xbas[4000]; double ybas[4000]; double zbas[4000]; void OnCalculate(); [ Last edited by _yixiao on 2010-6-3 at 20:01 ] |
回复此楼» 猜你喜欢
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有5人回复
-
博士读完未来一定会好吗
已经有29人回复
-
博士申请都是内定的吗？
已经有5人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有12人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
ICSD晶体结构数据的晶胞和文献上的不一样
已经有9人回复- 为什么优化晶胞时两次结果不同？
已经有11人回复
- 不同的结晶溶剂化物（晶胞参数相近）算不同晶体吗？
已经有5人回复
- CASTEP优化结构时超晶胞2x1x2和2x2x1几何优化结果一样吗？为什么
已经有4人回复
- 请问晶体结构一致其晶胞参数是否也一致？
已经有7人回复
- 求晶胞结构（六方相二氧化锗）
已经有10人回复
- 为什么晶胞体积大了就不能收结构？
已经有9人回复
- 如何看一个晶胞里面都几个分子组成？需要什么软件或者在线的查询网址？谢谢
已经有4人回复
- 【求助】怎样从建立三维周期晶胞，然后再切割成纳米结构？？
已经有4人回复
- 【求助】如何快速进行晶胞和原胞之间的结构参数的转换
已经有9人回复
- 【求助】用VMD怎样显示3D结构的晶胞呢？
已经有6人回复
- 【求助】体心立方结构间隙问题
已经有6人回复
- 【求助】大家帮帮忙——晶胞骨架结构和晶胞参数等
已经有4人回复
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
2楼2010-06-03 19:17:52
3楼2010-06-03 20:03:50
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
4楼2010-06-04 10:45:44