<div dir="ltr">MatCreateSubMatrices_xxx() calloc an array of size n+1, n=ismax and submat[n]=null as default.<div><br><div>For MatCreateSubMatrices_MPIAIJ(), I calloc ismax+nstages and use </div><div>submat[n] ... to store some datastructures for reuse.</div><div><br></div><div>When submat[n] != null, these data structures are destroyed.</div><div><br></div><div>Hong<br><div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jul 12, 2017 at 11:29 AM, Stefano Zampini <span dir="ltr"><<a href="mailto:stefano.zampini@gmail.com" target="_blank">stefano.zampini@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Is line number 6890 of matrix.c correct???<div><br></div><div><a href="https://bitbucket.org/petsc/petsc/src/b02a13eda23b9fabea7dbfcadcfb98e3258ee49d/src/mat/interface/matrix.c?at=master&fileviewer=file-view-default#matrix.c-6890" target="_blank">https://bitbucket.org/petsc/<wbr>petsc/src/<wbr>b02a13eda23b9fabea7dbfcadcfb98<wbr>e3258ee49d/src/mat/interface/<wbr>matrix.c?at=master&fileviewer=<wbr>file-view-default#matrix.c-<wbr>6890</a><span class="gmail-HOEnZb"><font color="#888888"><br clear="all"><div><br></div>-- <br><div class="gmail-m_4119919259534487923gmail_signature">Stefano</div>
</font></span></div></div>
</blockquote></div><br></div></div></div></div></div>