[petsc-dev] Changing Mat block size?

[Barry Smith]

  blocks are always square. 

  If the user wants non-square blocks then they go though the long involved construction process

   Barry

  Don’t make life so damn hard


On Feb 14, 2014, at 12:08 PM, Jed Brown <jed at jedbrown.org> wrote:

> Barry Smith <bsmith at mcs.anl.gov> writes:
> 
>>  Change MatCreateShell() to take the block size 
>> 
>>  Note that we should probably also change MatCreateSeqAIJ() and MatCreateAIJ() to take the block size now before a release.
> 
> What if the blocks are not square?  This just seems like an extraneous
> concept for most users.  Would we really be shooting ourselves in the
> foot if we just checked PetscAbs(map->bs) in any compatibility
> comparisons?  Then old code should work and users wouldn't have to deal
> with the concept unless it was relevant to them.