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

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

24小時(shí)熱門版塊排行榜

返回列表

咔咔咔咔

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 11.6
帖子: 35
在線: 26.3小時(shí)
蟲號(hào): 1348572
注冊(cè): 2011-07-17
專業(yè): 半導(dǎo)體材料

[求助] 二維ising模型蒙特卡洛程序出錯(cuò)，求指點(diǎn)

#include
#include
#include
#include
#define L 150
#define preN 100000
#define lateN 5100000
#define deltaN (lateN-preN)
double recordEnergy[deltaN];
int recordMagnet[deltaN];
int lattice[L][L];
int temp[L][L];
double energy;
int magnet;
double T;
struct CX
{
double C;
double X;
};

void InitialLattice();
void SaveInitial();
void GetInitial();
void CaculateFirstMagnet();
void CaculateFirstEnergy();
void GenerateNext();
struct CX Caculate_C_and_X();

int main()
{
int i;
struct CX cx[200];
FILE *ft,*fx,*fc;
ft=fopen("T.txt","w"

;
fx=fopen("X.txt","w"

;
fc=fopen("C.txt","w"

;
srand(time(NULL));

T = 1.0;
for(i=0;T<4.0;i++)
{
  InitialLattice();
  cx=Caculate_C_and_X();
  fprintf(ft,"%f\n",T);
  fprintf(fx,"%f\n",cx.X);
  fprintf(fc,"%f\n",cx.C);
  printf("T=%f ",T);
  printf("  X=%f ",cx.X);
  printf("  C=%f\n",cx.C);
      T = T+0.05;
}

fclose(ft);
fclose(fx);
fclose(fc);
return 0;
}

void InitialLattice()
{
int i,j;
for(i=0;i {
      for(j=0;j   {
if((rand()%10)>5) lattice[j]=1;
else             lattice[j]=1;
      }
}
}
void SaveInitial()
{
int k,m;
for(k=0;k {
  for(m=0;m temp[k][m]=lattice[k][m];
}
}
void GetInitial()
{
int k,m;
for(k=0;k {
  for(m=0;m lattice[k][m]=temp[k][m];

}
}
void CaculateFirstMagnet()
{
int i,j;
int iMagnet=0;
for(i=0;i {
      for(j=0;j iMagnet += lattice[j];
}
magnet = iMagnet;
}
void CaculateFirstEnergy()
{

double dEnergy = 0.0;
int up,down,right,left,px,py;
for(px=0;px {
  for(py=0;py   {
up = (L+px-1)%L;
down = (px+1)%L;
right = (L+py-1)%L;
left = (py+1)%L;
dEnergy += 0.5*lattice[px][py]*(lattice[up][py]+lattice[down][py]+lattice[px][right]+lattice[px][left]);
  }
}

energy = dEnergy;
}
void GenerateNext()
{
int px,py;
int dE,sum;
int up,down,right,left;
px = rand()%L;
py = rand()%L;

up = (L+px-1)%L;
down = (px+1)%L;
right = (L+py-1)%L;
left = (py+1)%L;

sum = (lattice[px][py])*(lattice[up][py]+lattice[down][py]+lattice[px][right]+lattice[px][left]);
switch(sum)
{
case 4: dE= 8; break;
case 2: dE= 4; break;
case 0: dE= 0; break;
case -2: dE= -4; break;
case -4: dE= -8; break;
default: printf("---Erro!---\n"

; break;
}

if( dE>0 )
{
      if((rand()*1.0/RAND_MAX) <= exp((-1.0)*dE/T))
  {
lattice[px][py] = -lattice[px][py];
magnet += 2*lattice[px][py];
energy += dE;
  }

}
else
{
  lattice[px][py] = -lattice[px][py];
  magnet += 2*lattice[px][py];
  energy += dE;
}

}
struct CX Caculate_C_and_X()
{
struct CX cx;
int i;
double AvgEnergy=0.0;
double AvgSquareEnergy=0.0;
double AvgMagnet=0.0;
double AvgSquareMagnet=0.0;

CaculateFirstMagnet();
CaculateFirstEnergy();

for(i=0; i {
      GenerateNext();
      if(i>=preN)
      {
         recordEnergy[i-preN] = energy;
         recordMagnet[i-preN] = magnet;
  }
}

for(i=0;i {
      AvgEnergy += (double)recordEnergy/deltaN;
      AvgSquareEnergy += (double)recordEnergy*recordEnergy/deltaN;
      AvgMagnet += (double)recordMagnet/deltaN;
      AvgSquareMagnet += (double)recordMagnet*recordMagnet/deltaN;
}

cx.C = (AvgSquareEnergy - AvgEnergy*AvgEnergy)/T/T;
cx.X = (AvgSquareMagnet - AvgMagnet*AvgMagnet)/T;

return cx;
}

:\microsoft visual studio\myprojects\ising\ising.cpp(44) : error C2440: '=' : cannot convert from 'struct CX' to 'struct CX [200]'
      No user-defined-conversion operator available that can perform this conversion, or the operator cannot be called
d:\microsoft visual studio\myprojects\ising\ising.cpp(46) : error C2228: left of '.X' must have class/struct/union type
d:\microsoft visual studio\myprojects\ising\ising.cpp(47) : error C2228: left of '.C' must have class/struct/union type
d:\microsoft visual studio\myprojects\ising\ising.cpp(49) : error C2228: left of '.X' must have class/struct/union type
d:\microsoft visual studio\myprojects\ising\ising.cpp(50) : error C2228: left of '.C' must have class/struct/union type
d:\microsoft visual studio\myprojects\ising\ising.cpp(67) : error C2440: '=' : cannot convert from 'const int' to 'int [150]'
      There are no conversions to array types, although there are conversions to references or pointers to arrays
d:\microsoft visual studio\myprojects\ising\ising.cpp(68) : error C2440: '=' : cannot convert from 'const int' to 'int [150]'
      There are no conversions to array types, although there are conversions to references or pointers to arrays
d:\microsoft visual studio\myprojects\ising\ising.cpp(98) : error C2297: '+=' : illegal, right operand has type 'int [150]'
d:\microsoft visual studio\myprojects\ising\ising.cpp(187) : error C2440: 'type cast' : cannot convert from 'double [5000000]' to 'double'
      There is no context in which this conversion is possible
d:\microsoft visual studio\myprojects\ising\ising.cpp(188) : error C2440: 'type cast' : cannot convert from 'double [5000000]' to 'double'
      There is no context in which this conversion is possible
d:\microsoft visual studio\myprojects\ising\ising.cpp(189) : error C2440: 'type cast' : cannot convert from 'int [5000000]' to 'double'
      There is no context in which this conversion is possible
d:\microsoft visual studio\myprojects\ising\ising.cpp(190) : error C2440: 'type cast' : cannot convert from 'int [5000000]' to 'double'
      There is no context in which this conversion is possible
執(zhí)行 cl.exe 時(shí)出錯(cuò).

ising.obj - 1 error(s), 0 warning(s)
怎么辦

回復(fù)此樓

» 猜你喜歡

0854電子信息求調(diào)劑已經(jīng)有7人回復(fù)
0854人工智能方向招收調(diào)劑已經(jīng)有3人回復(fù)
調(diào)劑已經(jīng)有4人回復(fù)
一志愿北化315 求調(diào)劑已經(jīng)有3人回復(fù)
有機(jī)合成求助已經(jīng)有5人回復(fù)
材料學(xué)碩333求調(diào)劑已經(jīng)有3人回復(fù)
289求調(diào)劑已經(jīng)有7人回復(fù)
收08調(diào)劑生已經(jīng)有6人回復(fù)
材料專碩找調(diào)劑已經(jīng)有5人回復(fù)
085602 289分求調(diào)劑已經(jīng)有5人回復(fù)

1樓 2012-04-24 20:54:59

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

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

人生轉(zhuǎn)折點(diǎn)，真心求幫助！已經(jīng)有13人回復(fù)
關(guān)于ANSYS在體模型的表面施加面載荷的問(wèn)題。已經(jīng)有9人回復(fù)
用c語(yǔ)言編程怎來(lái)表達(dá)時(shí)滯！謝謝！已經(jīng)有6人回復(fù)
請(qǐng)高手指點(diǎn)一下怎么確定一個(gè)符號(hào)矩陣中的等值元素已經(jīng)有6人回復(fù)
OPEN CV 幫助文檔參考手冊(cè) 已經(jīng)有31人回復(fù)
上海藥物所悲慘348分，求接收！�。。。。。。。。。�！已經(jīng)有36人回復(fù)
GPC（凝膠滲透色譜）譜圖求助已經(jīng)有9人回復(fù)
求高手指點(diǎn)：對(duì)于保守性較低的基因可以用簡(jiǎn)并引物來(lái)擴(kuò)嗎？已經(jīng)有9人回復(fù)
matlab程序運(yùn)行出錯(cuò)，求指點(diǎn)修改已經(jīng)有4人回復(fù)
山東大學(xué)有機(jī)333，求調(diào)劑已經(jīng)有8人回復(fù)
求助��！求一程序，用matlab程序做，用蒙特卡洛方法模擬已經(jīng)有10人回復(fù)
【求助】醇沉已經(jīng)有7人回復(fù)
【求助】藥學(xué)基礎(chǔ)知識(shí) 已經(jīng)有10人回復(fù)
求助C語(yǔ)言有限元一維、二維、矩形網(wǎng)格以及三角形網(wǎng)格程序已經(jīng)有4人回復(fù)

點(diǎn)擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

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