[petsc-users] Block ILU with AIJ matrix?

[Gong Ding]

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