<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">Ce Qin:</div><div class="gmail_quote">This is fixed in branch <span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:14px;white-space:pre-wrap">hzhang/matmatops_mpiaij4np1.</span></div><div class="gmail_quote"><font color="#333333"><span style="font-size:14px;white-space:pre-wrap">I consider it as a new feature - support mpiaij for np=1, not a bugfix.</span></font></div><div class="gmail_quote"><font color="#333333"><span style="font-size:14px;white-space:pre-wrap"><br></span></font></div><div class="gmail_quote"><font color="#333333"><span style="font-size:14px;white-space:pre-wrap">Once it passes our regression tests, I'll merge it to petsc master branch.</span></font></div><div class="gmail_quote"><font color="#333333"><span style="font-size:14px;white-space:pre-wrap">You may test it. Let us know you see any problem.</span></font></div><div class="gmail_quote"><font color="#333333"><span style="font-size:14px;white-space:pre-wrap">Thanks for reporting it.</span></font></div><div class="gmail_quote"><font color="#333333"><span style="font-size:14px;white-space:pre-wrap"><br></span></font></div><div class="gmail_quote"><font color="#333333"><span style="font-size:14px;white-space:pre-wrap">Hong</span></font></div><div class="gmail_quote"><font color="#333333"><span style="font-size:14px;white-space:pre-wrap"><br></span></font><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">Dear all,<div><br></div><div>I encountered a problem that MatMatMult crashes using only one processor. The matrix is MATMPIAIJ. I am using the master branch. The same program can run in parallel without problem. </div><div><br></div><div>Valgrind gives the following stack trace:</div><div><br></div><div><div><font face="monospace, monospace">==5543== Invalid read of size 8</font></div><div><font face="monospace, monospace">==5543==    at 0xDC094CE: MatGetBrowsOfAoCols_MPIAIJ (mpiaij.c:4509)</font></div><div><font face="monospace, monospace">==5543==    by 0xDC3FB95: MatMatMultSymbolic_MPIAIJ_MPIA<wbr>IJ_nonscalable (mpimatmatmult.c:218)</font></div><div><font face="monospace, monospace">==5543==    by 0xDC45193: MatMatMult_MPIAIJ_MPIAIJ (mpimatmatmult.c:61)</font></div><div><font face="monospace, monospace">==5543==    by 0xD9C657D: MatMatMult (matrix.c:9450)</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">==5543== </font></div><div><font face="monospace, monospace">==5543== Invalid read of size 4</font></div><div><font face="monospace, monospace">==5543==    at 0xDC094D7: MatGetBrowsOfAoCols_MPIAIJ (mpiaij.c:4510)</font></div><div><font face="monospace, monospace">==5543==    by 0xDC3FB95: MatMatMultSymbolic_MPIAIJ_MPIA<wbr>IJ_nonscalable (mpimatmatmult.c:218)</font></div><div><font face="monospace, monospace">==5543==    by 0xDC45193: MatMatMult_MPIAIJ_MPIAIJ (mpimatmatmult.c:61)</font></div><div><font face="monospace, monospace">==5543==    by 0xD9C657D: MatMatMult (matrix.c:9450)</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">==5543== </font></div><div><font face="monospace, monospace">==5543== Invalid read of size 4</font></div><div><font face="monospace, monospace">==5543==    at 0xDC0793B: MatGetBrowsOfAoCols_MPIAIJ (mpiaij.c:4513)</font></div><div><font face="monospace, monospace">==5543==    by 0xDC3FB95: MatMatMultSymbolic_MPIAIJ_MPIA<wbr>IJ_nonscalable (mpimatmatmult.c:218)</font></div><div><font face="monospace, monospace">==5543==    by 0xDC45193: MatMatMult_MPIAIJ_MPIAIJ (mpimatmatmult.c:61)</font></div><div><font face="monospace, monospace">==5543==    by 0xD9C657D: MatMatMult (matrix.c:9450)</font></div><div><font face="monospace, monospace"><br></font></div></div><div><div><font face="monospace, monospace">==5543== Invalid read of size 4</font></div><div><font face="monospace, monospace">==5543==    at 0xDC0794E: MatGetBrowsOfAoCols_MPIAIJ (mpiaij.c:4520)</font></div><div><font face="monospace, monospace">==5543==    by 0xDC3FB95: MatMatMultSymbolic_MPIAIJ_MPIA<wbr>IJ_nonscalable (mpimatmatmult.c:218)</font></div><div><font face="monospace, monospace">==5543==    by 0xDC45193: MatMatMult_MPIAIJ_MPIAIJ (mpimatmatmult.c:61)</font></div><div><font face="monospace, monospace">==5543==    by 0xD9C657D: MatMatMult (matrix.c:9450)</font></div><div><br></div></div><div>This problem may be related to this thread <a href="http://lists.mcs.anl.gov/pipermail/petsc-users/2014-September/022864.html" target="_blank">http://lists.mcs.anl.gov/<wbr>pipermail/petsc-users/2014-<wbr>September/022864.html</a>. If you need an example that reproduces this problem, please let me know.<br></div><div><br></div><div>Best,</div><div>Ce Qin</div></div>
</blockquote></div><br></div></div>