[petsc-dev] Changing Mat block size?

[Barry Smith]

  Ok, then don’t add bs to MatCreateShell and AIJ just go through the several steps to create and set up the matrices


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

> Barry Smith <bsmith at mcs.anl.gov> writes:
> 
>>  blocks are always square. 
>> 
>>  If the user wants non-square blocks then they go though the long involved construction process
> 
> non-square blocks show up in shell matrices for interpolation and
> restriction, which is not exactly a weird application.  I still think
> block size is an extraneous concept to most users, thus I don't think
> that adding it to all the convenience wrappers is a welcome or
> clearly-beneficial change.
> 
> What problems would we be inviting by making MatGetBlockSize return
> PetscAbs(rmap->bs)?  Obviously there are a number of places where we
> access the field directly, but those should use the interface anyway.