<div dir="ltr"><span style="font-size:12.8px">MatCreateSubMatrix_MPIAIJ_</span><span style="font-size:12.8px">Same<wbr>RowColDist() is a new addition to petsc for a scalable implementation of </span><span style="font-size:12.8px">MatCreateSubMatrix_MPIAIJ()<wbr>, thus it is still buggy.  </span><span style="font-size:12.8px">I'm fixing a bug appeared with IS blocksize. </span><div style="font-size:12.8px"><span style="font-size:12.8px"><br></span></div><div style="font-size:12.8px"><span style="font-size:12.8px">I cannot reproduce this valgrind error, neither from our nightly tests.</span></div><div style="font-size:12.8px"><span style="font-size:12.8px">Can I have your test that shows this error?</span></div><div style="font-size:12.8px"><span style="font-size:12.8px"><br></span></div><div style="font-size:12.8px">Hong</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jun 19, 2017 at 5:38 AM, Lisandro Dalcin <span dir="ltr"><<a href="mailto:dalcinl@gmail.com" target="_blank">dalcinl@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Hong, any chance this is related to your work around commit<br>
8b3fa1f7362a032c26473e2c80cb80<wbr>f09a5a21c7<br>
<br>
This test is working just fine with PETSc 3.5, 3.6, and 3.7, but fails<br>
with petsc/master<br>
<br>
$ git  clone <a href="https://github.com/dalcinl/petiga" rel="noreferrer" target="_blank">https://github.com/dalcinl/<wbr>petiga</a><br>
$ cd petiga<br>
$ make<br>
$ cd test<br>
$ make IGAInputOutput<br>
$ mpiexec -n 4 valgrind -q ./IGAInputOutput -iga_dim 2<br>
==26589== Invalid write of size 8<br>
==26589==    at 0x5A7B8E7: ISGetSeqIS_SameColDist_Private (mpiaij.c:3026)<br>
==26589==    by 0x5A7CFFF: MatCreateSubMatrix_MPIAIJ_<wbr>SameRowColDist<br>
(mpiaij.c:3130)<br>
==26589==    by 0x5A7E1DC: MatCreateSubMatrix_MPIAIJ (mpiaij.c:3243)<br>
==26589==    by 0x5B64553: MatCreateSubMatrix (matrix.c:7894)<br>
==26589==    by 0x4E9A704: MatView_MPI_IGA (petigamat.c:75)<br>
==26589==    by 0x5B19527: MatView (matrix.c:986)<br>
==26589==    by 0x402F5D: main (IGAInputOutput.c:119)<br>
<br>
<br>
<br>
<br>
--<br>
Lisandro Dalcin<br>
============<br>
Research Scientist<br>
Computer, Electrical and Mathematical Sciences & Engineering (CEMSE)<br>
Extreme Computing Research Center (ECRC)<br>
King Abdullah University of Science and Technology (KAUST)<br>
<a href="http://ecrc.kaust.edu.sa/" rel="noreferrer" target="_blank">http://ecrc.kaust.edu.sa/</a><br>
<br>
4700 King Abdullah University of Science and Technology<br>
al-Khawarizmi Bldg (Bldg 1), Office # 0109<br>
Thuwal 23955-6900, Kingdom of Saudi Arabia<br>
<a href="http://www.kaust.edu.sa" rel="noreferrer" target="_blank">http://www.kaust.edu.sa</a><br>
<br>
Office Phone: <a href="tel:%2B966%2012%20808-0459" value="+966128080459">+966 12 808-0459</a><br>
</blockquote></div><br></div>