24小时热门版块排行榜

返回列表

当前主题已经存档。

lansesanji

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 1287
散金: 349
红花: 1
帖子: 124
在线: 21.6小时
虫号: 606371
注册: 2008-09-18
专业: 油气储运工程

[交流] c语言求助

各位精通C语言的大虾，我有个问题，我的程序定义了将近10个数组，但是数组的维数如果定义大了就不能正常运行，维数如果小的话就没有问题，到底是怎么回事啊？帮帮忙啊，不胜感激啊……

回复此楼

» 猜你喜欢

参与限项已经有3人回复
假如你的研究生提出不合理要求已经有7人回复
实验室接单子已经有4人回复
全日制（定向）博士已经有4人回复
对氯苯硼酸纯化已经有3人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有12人回复
不自信的我已经有12人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复
北核录用已经有3人回复

蓝色三井

1楼 2009-12-23 14:53:32

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

shenao_07

木虫 (著名写手)

应助: 0 (幼儿园)
金币: 1581.3
散金: 726
红花: 2
帖子: 1350
在线: 146.9小时
虫号: 926332
注册: 2009-12-12
性别: GG
专业: 能源化工

这个你仔细检查下你定义的组数对不，还有语言编写有错误不，应该还是定义的时候出问题了

赞一下

回复此楼

2楼2009-12-23 17:54:05

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

sun_sun_sun

金虫 (小有名气)

博学EPI: 1
应助: 0 (幼儿园)
金币: 777.8
帖子: 135
在线: 28.5小时
虫号: 874679
注册: 2009-10-16
性别: GG
专业: 电气

建议LZ把程序贴上来，大家才能帮你分析，有时问题可能出在别的语句上。

赞一下

回复此楼

3楼2009-12-23 18:23:53

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

lansesanji

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 1287
散金: 349
红花: 1
帖子: 124
在线: 21.6小时
虫号: 606371
注册: 2008-09-18
专业: 油气储运工程

现在贴上源程序，大家帮帮忙了！

/* Note:Your choice is C IDE */
#include
#include
#define k 6
#define numb 67
void main()
{
  double T[numb][numb],detar,detax,aw,ae,an,as,ap,xw,xe,rn,rs,e1=0,e2=0,b[numb][numb],a[numb][numb]={0};
  double r[numb],x[numb],t[numb][numb];
  double T0[7][numb][numb]={0},a0[7][numb][numb]={0};
  int i,j,m,n,Nk;
  Nk=pow(2,k)+1;
  for(i=0;i=100;
  for(j=0;j   for(i=1;i<=Nk;i++) T[Nk]=300;
  for(j=1;j<=Nk;j++) T[j][0]=400;
  for(i=1;i    {for(j=1;j         T[j]=200;
   }
for(i=0;i<=Nk;i++)
{for(j=0;j<=Nk;j++)
      T0[1][j]=T[j];
   }
detar=1.0/Nk;
detax=1.0/Nk;
for(i=0;i<=Nk;i++) r=x=detar*i;
do{
e1=0;
e2=0;
for(i=0;i<=Nk;i++)
   {for(j=0;j<=Nk;j++)
         t[j]=T[j];
         }
for(n=1;n<7;n++)
      {
      for(m=0;m<3;m++)
         {for(i=1;i             {for(j=1;j                    {xw=x-x[i-1];
                  xe=x[i+1]-x;
                  rn=r[j+1]-r[j];
                  rs=r[j]-r[j-1];
                  ap=2/xw/xe+2/rs/rn+1/r[j]/rn;
                      aw=2/(xw*xw+xe*xw);
                      ae=2/(xe*xe+xe*xw);
                      as=2/(rs*rs+rs*rn);
                      an=2/(rn*rn+rs*rn)+1/r[j]/rn;
                      T[j]=aw/ap*T[i-1][j]+as/ap*T[j-1]+ae/ap*T[i+1][j]+an/ap*T[j+1]+a0[n-1][j]/ap;
                      T0[n][j]=T[j];
                      }
               }
         }

      for(i=1;i          {for(j=1;j                {xw=x-x[i-1];
               xe=x[i+1]-x;
               rn=r[j+1]-r[j];
               rs=r[j]-r[j-1];
               ap=2/xw/xe+2/rs/rn+1/r[j]/rn;
                   aw=2/(xw*xw+xe*xw);
                   ae=2/(xe*xe+xe*xw);
                   as=2/(rs*rs+rs*rn);
                   an=2/(rn*rn+rs*rn)+1/r[j]/rn;
                   b[j]=a0[n-1][j]-ap*T[j]+aw*T[i-1][j]+as*T[j-1]+ae*T[i+1][j]+an*T[j+1];
               }
            }

         Nk=(Nk+1)/2;
         r[0]=x[0]=0;
         r[1]=x[1]=pow(2,(n-2))*detar+0.5*detar;
         for(i=2;i=x=r[1]+(i-1)*(1-2*r[1])/(Nk-2);
         r[Nk]=x[Nk]=1;
         for(i=1;i             {for(j=1;j                {a[j]=(b[2*i-1][2*j]+b[2*i-1][2*j-1]+b[2*i][2*j]+b[2*i][2*j-1])/4;
               a0[n][j]=a[j];
               }

            }
         for(i=0;i<=Nk;i++)
            {for(j=0;j<=Nk;j++)
            T[j]=0;
         }
}

Nk=Nk*2-1;
for(n=5;n>0;n--)
   {
      Nk=Nk*2-1;
      if (Nk==pow(2,k)+1) for(i=0;i<=Nk;i++) r=x=detar*i;
      else {r[0]=x[0]=0;
            r[1]=x[1]=pow(2,(n-2))*detar+0.5*detar;
            for(i=2;i=x=r[1]+(i-1)*(1-2*r[1])/(Nk-2);
            r[Nk]=x[Nk]=1;
      }
      for(i=1;i<(Nk/2);i++)
         {for(j=1;j<(Nk/2);j++)
         {T0[n][2*i-1][2*j]=T0[n+1][j]+T0[n][2*i-1][2*j];
         T0[n][2*i-1][2*j-1]=T0[n+1][j]+T0[n][2*i-1][2*j-1];
         T0[n][2*i][2*j]=T0[n+1][j]+T0[n][2*i][2*j];
         T0[n][2*i][2*j-1]=T0[n+1][j]+T0[n][2*i][2*j-1];
   }
         }

      for(m=0;m<3;m++)
         {for(i=1;i             {for(j=1;j                    {xw=x-x[i-1];
                  xe=x[i+1]-x;
                  rn=r[j+1]-r[j];
                  rs=r[j]-r[j-1];
                  ap=2/xw/xe+2/rs/rn+1/r[j]/rn;
                      aw=2/(xw*xw+xe*xw);
                      ae=2/(xe*xe+xe*xw);
                      as=2/(rs*rs+rs*rn);
                      an=2/(rn*rn+rs*rn)+1/r[j]/rn;
                      T0[n][j]=aw/ap*T0[n][i-1][j]+as/ap*T0[n][j-1]+ae/ap*T0[n][i+1][j]+an/ap*T0[n][j+1]+a0[n][j]/ap;
                   }
            }
         }
      for(i=0;i<=Nk;i++)
         {for(j=0;j<=Nk;j++)
               T[j]=T0[1][j];
         }

      }
for(i=0;i<=Nk;i++)
   {for(j=0;j<=Nk;j++)
         {e2=fabs(t[j]-T[j]);
            if (e2>e1) e1=e2;
            }
      }
}while (e1>0.001);
for(i=0;i<=Nk;i++)
{for(j=0;j<=Nk;j++)
(printf("T[%d][%d]=%lf ",i,j,T[j]));
}
}