[petsc-users] how to switch ILU to its block version?

[Barry Smith]

  Fande,

    The MATSEQBAIJ matrix automatically uses a point block ILU; with the block size of the matrix. So if your matrix naturally has a block size bs if you use the MATSEQBAIJ matrix type it will be point block ILU of block size bs.

    We don't have a variable point block ILU in PETSc for SeqAIJ. It is possible to write but would be rather involved.

   Barry

> On Mar 1, 2015, at 1:47 PM, Fande Kong <fd.kong at siat.ac.cn> wrote:
> 
> Hi all,
> 
> I have a code with ILU as subdomain solver, and it works well. Now I want to try a block version of ILU, there are some easy ways to implement this?
> 
> Fande,