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 <[email protected]> wrote: > 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.
