|
|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
for(i=0;i<900;i++)
t=t_h;
epsilon=1e-4;
do
{for(i=0;i<900;i++)
{t1=c;
for(j=0;j<i;j++)
t1=t1-a[i+j*900]*t1[j];
for(j=i+1;j<900;j++)
t1=t1-a[i+j*900]*t[j];
t1=t1/a[i+j*900];
}
norm=0.0;
for(i=0;i<900;i++)
norm=norm+(t1-t)*(t1-t);
norm=sqrt(norm);
for(i=0;i<900;i++)
t=t1;
}
while(norm>epsilon);
用高斯塞的尔迭代.求温度场
t1=t1-a[i+j*900]*t1[j];这句出毛病了显示Floating point overflow
是怎么回事呢.请高手帮忙.
|
|