[petsc-dev] Bug in MatGetSubMatrix

[Hong]

Stefano:

> There's a bug in MatGetSubMatrix; the case MAT_REUSE_MATRIX is not covered
> by the interface when  isrow == iscol and isrow is a stride 1 IS
> representing the whole set of rows. A tentative patch is attached.
>

Indeed, we do not have this case. I'll create a branch from master and test
your patch.

Hong
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20150309/d2615d82/attachment.html>