<div dir="ltr">Dear Hong,<div><br></div><div>I have tested the new feature, it works now.</div><div>Thanks.</div><div><br></div><div>Best,</div><div>Ce Qin</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2017-03-17 23:58 GMT+08:00 Hong <span dir="ltr"><<a href="mailto:hzhang@mcs.anl.gov" target="_blank">hzhang@mcs.anl.gov</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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_<wbr>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><span class="HOEnZb"><font color="#888888"><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></font></span><div><div class="h5"><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/piper<wbr>mail/petsc-users/2014-Septembe<wbr>r/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></div></div>
</blockquote></div><br></div>