Hi All, I am solving a nonlinear system whose Jacobian matrix has a block structure. More precisely, there is a mesh, and for each vertex there are 11 variables associated with it. I am using BAIJ.
I thought block ILU(k) should be more efficient than the point-wise ILU(k). After some numerical experiments, I found that the block ILU(K) is much slower than the point-wise version. Any thoughts? Fande,