| 4 | 1/1 | 返回列表 |
| 查看: 657 | 回復: 3 | |||
_yixiao鐵蟲 (小有名氣)
|
[交流]
【求助】建立體心晶胞結構代碼,運行不出結果! 已有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的晶格常數(shù) 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方向的重復次數(shù),其中M任意整數(shù) 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); } 這是頭文件中定義的變量及函數(shù): #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 ] |

木蟲 (著名寫手)
鐵蟲 (小有名氣)

木蟲 (著名寫手)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856調劑,是學校就去 +6 | sllhht 2026-03-19 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 294求調劑材料與化工專碩 +13 | 陌の森林 2026-03-18 | 13/650 |
|
|
[考研] 材料專碩274一志愿陜西師范大學求調劑 +8 | 薛云鵬 2026-03-13 | 8/400 |
|
|
[考研] 085601材料工程專碩求調劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 材料與化工求調劑 +7 | 為學666 2026-03-16 | 7/350 |
|
|
[考研] 085600材料與化工調劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 346求調劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|
[考研] 344求調劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 295求調劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
[考研] 085600材料與化工 +5 | 安全上岸! 2026-03-16 | 5/250 |
|
|
[考研] 302求調劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 生物學071000 329分求調劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
[考研] 301求調劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 326求調劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 308求調劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 283求調劑 +3 | 聽風就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 一志愿211 0703方向310分求調劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 321求調劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 085601材料工程315分求調劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|