[petsc-users] Block ILU with AIJ matrix?

[behzad baghapour]

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
paralleled.

http://www-users.cs.umn.edu/~saad/software/ITSOL/index.html

B.B.

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
https://sites.google.com/site/behzadbaghapour
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>