[petsc-dev] MetSetBlockSize issue

Jed Brown jedbrown at mcs.anl.gov
Wed Apr 25 18:15:38 CDT 2012


On Wed, Apr 25, 2012 at 17:58, Mark F. Adams <mark.adams at columbia.edu>wrote:

> The blocks are not always dense, so by what appears to be your definition
> of block size it is not (always) a 'column block size'.   But I think that
> it is  semantically a blocked matrix and hence it has a column block size.


I think we shouldn't try to encode the near-null space being sparse. For
elasticity, we would only reduce the 18 matrix entries to 12, which is only
about break-even on storage since the dense null space can use fewer column
indices (this is ignoring the dense case also being more efficient).
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120425/0d3a2187/attachment.html>


More information about the petsc-dev mailing list