[petsc-users] BILU for AIJ matrix

[Barry Smith]

> On Mar 17, 2015, at 2:32 PM, Chung-Kan Huang <ckhuangf at gmail.com> wrote:
> 
> Hi,
> 
> It is true that PCILU constructs only pointed ILU when AIJ is used.
> 
> Will it matter if I use MATSetBlockSize to a AIJ with block size > 1?
> 
> Is it true that the only benefit I will get from using MATSetBlockSize is to be able to use MatSetValuesBlocked?

   Yes, for AIJ that is pretty much the only benefit. But if you just start with a BAIJ matrix you get 1) the point block ILU, 2) less memory usage, 3) faster matrix operations so all around BAIJ wins when you have a natural block size > 1. 

  Barry

> 
> 
> Kan
> 
> 
> -- 
> Cheers
>