#2
林月儿2015-05-24 16:19
|
我编写的是一个子程序代码,程序结构大体上是这样的:
fo 我编写的是一个子程序代码,程序结构大体上是这样的:
for(int i=M;i<M+1;i++) {
for(int j=N;j<N+1;j++) {
for(int k=P-1;k<P;k++) {
udutox=u[i][j][k]/(1*dx)*(u[i][j][k]-u[i-1][j][k]);
vdutoy=v[i][j][k]/(1*dy)*(u[i][j][k]-u[i][j-1][k]);
wdutoz=w[i][j][k]/(2*dz)*(u[i][j][k+1]-u[i][j][k-1]);
dutotox=(u[i][j][k]-2*u[i-1][j][k]+u[i-2][j][k])/(dx*dx);
dutotoy=(u[i][j][k]-2*u[i][j-1][k]+u[i][j-2][k])/(dy*dy);
dutotoz=(u[i][j][k+1]-2*u[i][j][k]+u[i][j][k-1])/(dz*dz);
ustar[i][j][k]=u[i][j][k]-dt*(udutox+vdutoy+wdutoz-1/Re*(dutotox+dutotoy+dutotoz));
}
}
}
功能就是为了求解三维数组ustar[i][j][k]的值,一共是有125段这样的fou循环,这个是第124个循环,编译错误就指示在这个循环上,请问可能是什么问题?