[petsc-dev] Changing Mat block size?

Barry Smith bsmith at mcs.anl.gov
Fri Feb 14 12:46:34 CST 2014

  blocks are always square. 

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


  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.

More information about the petsc-dev mailing list