<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">Stefano:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">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. </div></blockquote><div> </div><div>Indeed, we do not have this case. I'll create a branch from master and test your patch.</div><div><br></div><div>Hong</div></div><br></div></div>