24小时热门版块排行榜

返回列表

咔咔咔咔

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 11.6
帖子: 35
在线: 26.3小时
虫号: 1348572
注册: 2011-07-17
专业: 半导体材料

[求助] 二维ising模型蒙特卡洛程序出错，求指点

#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
执行 cl.exe 时出错.

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

回复此楼

» 猜你喜欢

三无产品还有机会吗已经有4人回复
投稿返修后收到这样的回复，还有希望吗已经有7人回复
压汞仪和BET测气凝胶孔隙率已经有4人回复
博士申请都是内定的吗？已经有14人回复
谈谈两天一夜的“延安行” 已经有13人回复
氨基封端PDMS和HDI反应快速固化已经有11人回复
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？已经有11人回复
论文投稿求助已经有4人回复
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛已经有3人回复
投稿精细化工已经有6人回复

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主咔咔咔咔的主题更新

返回列表