【求助】代码问题

7楼

double rho[NX+1][NY+1],u[NX+1][NY+1][2],u0[NX+1][NY+1][2],f[NX+1][NY+1][Q]
设置的边界是 NX和 NY，再定义上面 [NX+1][NY+1][是什么意思？

sjtugxl

0-NX是NX+1个数

hellolegend

楼上正解

同问
何雅玲老师的顶盖驱动流动的代码中初始化函数里面有下面一句
for (i =0; i<NX; i++)
{
   for (j =0; j<NY; j++)
   {
       u[ i ][ NY ][0]=U;
    }
}
这句怎么表示顶盖的速率赋值为U? 好像表示的是右边的节点其x方向的速率为U

[ 本帖最后由 hellolegend 于 2011-5-9 08:17 编辑 ]

nanosmile

按C语言的格式，U[NX+1][NY+1]，表示NX+1行，NY+1列，
在XY平面，水平为X，垂直为Y
i=0--NX表示行数，即Y坐标，
j=0--NY表示列数，即X坐标，
跟3楼的观点相同。
但是在编程中，有一个思维转换的过程，你可以将整个U[NX+1][NY+1]的矩阵逆时针转动90°，就变成我们习惯的坐标系了。何老师的程序也是这么个意思。

erhu10

没看明白思维转换的那部分，想请教一下，何老师的那句代码是不是只用下边的表述就可以了，
for(i=0;i<=NX;i++)
{
    u[ i ][NY][0]=U;
}

[ 本帖最后由 erhu10 于 2012-4-8 15:47 编辑 ]