[petsc-users] Block ILU with AIJ matrix?

[Matthew Knepley]

2011/10/31 Gong Ding <gdiso at ustc.edu>

> 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?


I do not believe any sparse factorization package handles varying block
size.

   Matt


>
> Gong Ding

-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111031/063865e0/attachment.htm>