Hi, I have several materials, each has different dofs. It is a block matrix with nonuniform block size. However, I know it is impossible to set different block size within a BAIJ matrix -- sigh. I can only use AIJ matrix here. My question is let a ILU preconditioner of AIJ matrix do block factorization?
I guess the block ILU can be 1) more efficient since it is cache friendly. 2) more accurate since it has a "matched" structure. Or any suggestion? Gong Ding
