压力的泊松方程

singularity

我正在做PCG的流体模拟,以前用Jacobi方法做的.现在需要压力泊松的矩阵形式.但是我推倒得到的形式是:

-6r        r        0        0            p0        div0(vel)

   r     -6r        r         0            p1        div1(vel)
                                       x             =
   0       r     -6r         r             p2        div2(vel)
  
   0      0        r      -6r             p3        div3(vel)

不知道是否正确.
还有个疑问:
对于不可压缩流体r为常数.既然这样这个方程求解就相当容易了,那为什么存在求解困难等问题呢

[ 本帖最后由 singularity 于 2009-2-12 07:33 编辑 ]