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>
