Barry Smith <[email protected]> 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.
pgpluIznchC_7.pgp
Description: PGP signature
