版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

xcid

銀蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 454.9
散金: 11
帖子: 188
在線: 52.6小時
蟲號: 1093658
注冊: 2010-09-09
性別: GG
專業(yè): 基礎(chǔ)物理學(xué)

[求助] 大家能不能幫我看看這個程序，我用的是DEV-C++，但是一運(yùn)行就報錯

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <iostream.h>
#include <fstream.h>

#define NKPTS_MAX 10000
#define NBANDS_MAX 200

int main ( int argc, char *argv[])
{
if( argc < 4)
{
printf("should be like: \n"

;
printf("./a.out nkx nky nkz \n"

;
exit(0);
  }
  int nkx, nky,nkz;
  nkx=atoi(argv[1]);
  nky=atoi(argv[2]);
  nkz=atoi(argv[3]);

  fstream fid;
  fid.open("OUTCAR",ios::in);
  if ( !fid  )
  {
printf("cant open OUTCAR\n"

;
exit(0);
}

char str[255], stmp[255];

//keyword1[]="NKPTS";
int NKPTS,NBANDS;
do {
fid.getline(str,255);
if(strstr(str,"NKPTS"

!= NULL )
{
printf("str=%s\n", str);
   sscanf(str,"%s%s%s%d",stmp,stmp,stmp,&NKPTS);
   sscanf(str,"%s%s%s%s%s%s%s%s%s%s%s%s%s%s%d",stmp,stmp,stmp,stmp,stmp,stmp,stmp,stmp,stmp,stmp,stmp,stmp,stmp,stmp,&NBANDS);
   printf("NKPTS: %d, NBANDS: %d\n",NKPTS,NBANDS);
   break;
}
}while( !fid.eof()) ;

float A[3][3], B[3][3];
do {
fid.getline(str,255);
if(strstr(str,"reciprocal lattice vectors"

!= NULL )
{
   for(int i=0;i<3;i++)
   {
fid.getline(str,255);
      sscanf(str,"%f%f%f%f%f%f",&A[0],&A[1],&A[2],&B[0],&B[1],&B[2]);
   }
   printf("A1: %f %f %f \n",A[0][0],A[0][1],A[0][2]);
   printf("A2: %f %f %f \n",A[1][0],A[1][1],A[1][2]);
   printf("A3: %f %f %f \n",A[2][0],A[2][1],A[2][2]);
   printf("B1: %f %f %f \n",B[0][0],B[0][1],B[0][2]);
   printf("B2: %f %f %f \n",B[1][0],B[1][1],B[1][2]);
   printf("B3: %f %f %f \n",B[2][0],B[2][1],B[2][2]);

   break;
}
}while( !fid.eof());

float Ef;
do {
fid.getline(str,255);
if(strstr(str,"E-fermi" != NULL )
{
      printf("str=%s\n", str);
   sscanf(str,"%s%s%f",stmp,stmp,&Ef);
   printf("Ef: %f\n",Ef);
   break;
}
}while( !fid.eof()) ;

fid.getline(str,255);
fid.getline(str,255);

float kpt[NKPTS_MAX][3];
float Egv[NKPTS_MAX][NBANDS_MAX];
for(int ik =0; ik < NKPTS; ik++)
{
      // k-point 1: kx ky kz
      fid.getline(str,255);
      sscanf(str,"%s%s%s%f%f%f",stmp,stmp,stmp,&kpt[ik][0],&kpt[ik][1],&kpt[ik][2]);
      fid.getline(str,255);

for(int ib=0;ib<NBANDS;ib++)
{
int bandno,occ;
      fid.getline(str,255);
            sscanf(str,"%d%f%f",&bandno,&Egv[ik][ib],&occ);
   //       printf("%d %d : %f\n", ik,ib,Egv[ik][ib]);
}

      fid.getline(str,255);
}

  fid.close();

  // find those cross Ef

int Nf=0;
static int nstore[NBANDS_MAX];
for( int ib =0 ; ib < NBANDS; ib ++)
  {
float vmin=10000, vmax=-10000;
      for( int ik =0 ; ik < NKPTS ; ik++ )
      {
if ( Egv[ik][ib] < vmin)
               { vmin = Egv[ik][ib];}
if( Egv[ik][ib] > vmax )
  { vmax= Egv[ik][ib];}
}
if( vmax > Ef && vmin < Ef)
{ nstore[Nf++]=ib; }

  }
  int ib;
  float Egvf[NKPTS][Nf];
  for(int i=0; i<Nf; i++)
  {
ib=nstore;
      for( int ik=0 ; ik < NKPTS ; ik ++ )
{
      Egvf[ik]=Egv[ik][ib];
}
  }

//save to bxsf

FILE *fid2;
fid2=fopen("FS.bsxf","w";
if( fid2 == NULL)
{
   printf( "can't create fout: FS.bsxf\n";
   exit(0);
}

fprintf(fid2,"BEGIN_INFO\n";
fprintf(fid2,"  Fermi Energy: %7.5f\n",Ef);
fprintf(fid2,"END_INFO\n";

fprintf(fid2,"\n";

fprintf(fid2,"BEGIN_BLOCK_BANDGRID_3D\n";
fprintf(fid2,"  one_word_comment \n";
fprintf(fid2,"  BEGIN_BANDGRID_3D\n";
fprintf(fid2,"  %d \n",Nf);
fprintf(fid2," %d %d %d \n",nkx,nky,nkz);
fprintf(fid2," 0 0 0\n"; // origin Gamma
fprintf(fid2," %f %f %f \n", B[0][0],B[0][1],B[0][2]);//Vec(1,);
fprintf(fid2," %f %f %f \n", B[1][0],B[1][1],B[1][2]);
fprintf(fid2," %f %f %f \n", B[2][0],B[2][1],B[2][2]);

int ik = 0;
for(int i=0; i < Nf; i++)
{
fprintf(fid2, " BAND: %d \n",nstore+1);
for(int ix=0 ; ix < nkx; ix ++)
{
for(int  iy=0 ; iy < nky; iy ++)
      {
         for(int iz=0 ; iz <nkz; iz++)
{
            fprintf(fid2,"  %f  ",Egvf[ik++]);
         }
         fprintf(fid2,"\n";
      }
      fprintf(fid2,"\n";
}

}

fprintf(fid2,"  END_BANDGRID_3D\n";
fprintf(fid2,"END_BLOCK_BANDGRID_3D\n";

fclose(fid2);

}

回復(fù)此樓

» 猜你喜歡

材料與化工328分調(diào)劑已經(jīng)有3人回復(fù)
279分求調(diào)劑一志愿211 已經(jīng)有20人回復(fù)
293求調(diào)劑已經(jīng)有3人回復(fù)
一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271 已經(jīng)有6人回復(fù)
085404電子信息284分求調(diào)劑已經(jīng)有3人回復(fù)
306求0703調(diào)劑一志愿華中師范已經(jīng)有9人回復(fù)
一志愿吉大化學(xué)322求調(diào)劑已經(jīng)有6人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)
344求調(diào)劑已經(jīng)有3人回復(fù)
一志愿哈工大，085400，320，求調(diào)劑已經(jīng)有3人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

哪位大仙幫我看下程序，我用的matlab7.1，運(yùn)行提示du未定義已經(jīng)有11人回復(fù)
1stopt運(yùn)行程序，總是出現(xiàn)字符函數(shù)錯誤的提示，請高手幫助已經(jīng)有5人回復(fù)
配置一臺1萬元以內(nèi)的計(jì)算機(jī)跑linux系統(tǒng)，運(yùn)行C語言分子動力學(xué)模擬已經(jīng)有9人回復(fù)
求助matlab一程序運(yùn)行的問題，幫忙看看哪里不對已經(jīng)有4人回復(fù)
MC一個簡單的程序，運(yùn)行不了，哪位大俠幫看看怎么回事。多謝！已經(jīng)有7人回復(fù)
程序沒報錯中斷，但停止運(yùn)行已經(jīng)有3人回復(fù)
linux下運(yùn)行python報錯，不知原因已經(jīng)有12人回復(fù)
OpenGl 程序運(yùn)行報錯已經(jīng)有4人回復(fù)
為什么這個程序運(yùn)行時有延遲，在VC++6.0中，就是點(diǎn)運(yùn)行要過一會兒才提示輸入已經(jīng)有8人回復(fù)
【求助】dev c++編譯結(jié)束，運(yùn)行后怎么樣輸出結(jié)果？已經(jīng)有6人回復(fù)
【求助】為什么我這個程序編譯都沒有問題，就是運(yùn)行不出來？【已完結(jié)】已經(jīng)有9人回復(fù)

小菜一碟

1樓 2013-05-24 10:15:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ioadong

木蟲 (著名寫手)

應(yīng)助: 55 (初中生)
金幣: 4486
紅花: 20
帖子: 1656
在線: 674.1小時
蟲號: 928807
注冊: 2009-12-15
性別: GG
專業(yè): 有機(jī)合成

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

main函數(shù)這么長？嘗試先拆成功能小函數(shù)，這樣好定位錯誤。

贊一下

回復(fù)此樓

以科學(xué)的理性代替天然的非理性。

2樓2013-05-24 21:53:44

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 xcid 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 資源與環(huán)境調(diào)劑申請(333分) +6	holy J 2026-03-21	6/300	2026-03-24 09:51 by hengsmith
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3	開開心心沒煩惱 2026-03-23	4/200	2026-03-24 00:06 by 開開心心沒煩惱
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +8	皮皮皮123456 2026-03-22	8/400	2026-03-23 23:47 by 熱情沙漠
[考研] 一志愿國科過程所081700，274求調(diào)劑 +3	三水研0水立方 2026-03-23	3/150	2026-03-23 23:11 by MajorWen
[考研] 0854 考研調(diào)劑招生了！AI 方向 +4	pk3725069 2026-03-19	16/800	2026-03-23 23:09 by 汪��？！
[考研] 316求調(diào)劑 +7	梁茜雯 2026-03-19	7/350	2026-03-23 16:21 by lingjue
[考研] 323求調(diào)劑 +6	洼小桶 2026-03-18	6/300	2026-03-23 00:29 by king123！
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +4	沿岸有貝殼6 2026-03-21	4/200	2026-03-22 07:21 by ilovexiaobin
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +4	lemonzzn 2026-03-17	8/400	2026-03-21 22:49 by lemonzzn
[考研] 材料學(xué)碩301分求調(diào)劑 +7	Liyouyumairs 2026-03-21	7/350	2026-03-21 22:31 by peike
[考研] 材料工程專碩 348分求調(diào)劑 +3	冬辭. 2026-03-17	5/250	2026-03-21 18:47 by 學(xué)員8dgXkO
[考研] 297求調(diào)劑 +11	戲精丹丹丹 2026-03-17	12/600	2026-03-21 17:47 by ColorlessPI
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 274求調(diào)劑 +10	S.H1 2026-03-18	10/500	2026-03-20 23:51 by JourneyLucky
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9	Ncdx123456 2026-03-19	9/450	2026-03-20 23:41 by lovewei0727
[考研] 330求調(diào)劑 +4	小材化本科 2026-03-18	4/200	2026-03-20 23:13 by JourneyLucky
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4	@taotao 2026-03-19	4/200	2026-03-20 22:14 by JourneyLucky
[考研] 295復(fù)試調(diào)劑 +8	簡木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4	xbxudjdn 2026-03-18	4/200	2026-03-20 09:07 by 不168
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频