[petsc-users] Block ILU with AIJ matrix?

behzad baghapour behzad.baghapour at gmail.com
Mon Oct 31 00:37:57 CDT 2011

I think you may choose the maximum block size as a BAIJ dimension. However
ITSOL has variable BLOCK ILU routine which may help you BUT it is not



On Mon, Oct 31, 2011 at 7:37 AM, Matthew Knepley <knepley at gmail.com> wrote:

> 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

Behzad Baghapour
Ph.D. Candidate, Mechecanical Engineering
University of Tehran, Tehran, Iran
Fax: 0098-21-88020741
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111031/b19020c2/attachment.htm>

More information about the petsc-users mailing list