<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">Stefano:</div><div class="gmail_quote">The bug is fixed on branch hzhang/fix-ptap-seqaij</div><div class="gmail_quote"><a href="https://bitbucket.org/petsc/petsc/commits/4bcca364924c81402b2f7ad375877c5546b7966b">https://bitbucket.org/petsc/petsc/commits/4bcca364924c81402b2f7ad375877c5546b7966b</a><br></div><div class="gmail_quote"><br></div><div class="gmail_quote">Satish,</div><div class="gmail_quote">Can you help merge this branch to next for testing?</div><div class="gmail_quote"><br></div><div class="gmail_quote">Hong</div><div class="gmail_quote"><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"><div>I have found a bug in  MatPtAPNumeric_SeqAIJ_<wbr>SeqAIJ_SparseAxpy.</div><div><div>Attached is a MWE (using a modified version of mat/examples/tests/ex94.c) to reproduce it.</div><div><br></div><div>As far as I can tell, the symbolic step computes an rmax for C equal to 4.</div></div><div>Instead, the numerical step tries to write 1+rmax at apj in that function.</div><div>Don't know what would be the proper fix.</div><div><br></div><div><br></div><div><div>[szampini@KW14141 tests]$ ./ex94 -f0 Amat -f1 Amat -Pmat Pmat</div><div>[0]PETSC ERROR: PetscTrFreeDefault() called from MatPtAPNumeric_SeqAIJ_SeqAIJ_<wbr>SparseAxpy() line 270 in /home/szampini/Sources/petsc/<wbr>src/mat/impls/aij/seq/matptap.<wbr>c</div><div>[0]PETSC ERROR: Block [id=0(16)] at address 0x22399e0 is corrupted (probably write past end of array)</div><div>[0]PETSC ERROR: Block allocated in MatPtAPNumeric_SeqAIJ_SeqAIJ_<wbr>SparseAxpy() line 208 in /home/szampini/Sources/petsc/<wbr>src/mat/impls/aij/seq/matptap.<wbr>c</div><div>[0]PETSC ERROR: --------------------- Error Message ------------------------------<wbr>------------------------------<wbr>--</div><div>[0]PETSC ERROR: Memory corruption: <a href="http://www.mcs.anl.gov/petsc/documentation/installation.html#valgrind" target="_blank">http://www.mcs.anl.gov/petsc/<wbr>documentation/installation.<wbr>html#valgrind</a></div><div>[0]PETSC ERROR: Corrupted memory</div><div>[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank">http://www.mcs.anl.gov/petsc/<wbr>documentation/faq.html</a> for trouble shooting.</div><div>[0]PETSC ERROR: Petsc Development GIT revision: v3.7.2-1249-ga13144f  GIT Date: 2016-08-25 12:12:30 +0300</div><div>[0]PETSC ERROR: ./ex94 on a debug named KW14141 by szampini Mon Sep 26 10:25:06 2016</div><div>[0]PETSC ERROR: Configure options --download-c2html --download-ctetgen --download-exodusii --download-hdf5 --download-hypre --download-metis --download-ml --download-mumps --download-netcdf --download-p4est --download-parmetis --download-pastix --download-ptscotch --download-scalapack --download-suitesparse --download-triangle --with-debugging=1 PETSC_ARCH=debug</div><div>[0]PETSC ERROR: #1 PetscTrFreeDefault() line 298 in /home/szampini/Sources/petsc/<wbr>src/sys/memory/mtr.c</div><div>[0]PETSC ERROR: #2 MatPtAPNumeric_SeqAIJ_SeqAIJ_<wbr>SparseAxpy() line 270 in /home/szampini/Sources/petsc/<wbr>src/mat/impls/aij/seq/matptap.<wbr>c</div><div>[0]PETSC ERROR: #3 MatPtAP_SeqAIJ_SeqAIJ() line 42 in /home/szampini/Sources/petsc/<wbr>src/mat/impls/aij/seq/matptap.<wbr>c</div><div>[0]PETSC ERROR: #4 MatPtAP() line 9257 in /home/szampini/Sources/petsc/<wbr>src/mat/interface/matrix.c</div><div>[0]PETSC ERROR: #5 main() line 259 in /home/szampini/Sources/petsc/<wbr>src/mat/examples/tests/ex94.c</div><div>[0]PETSC ERROR: PETSc Option Table entries:</div><div>[0]PETSC ERROR: -f0 Amat</div><div>[0]PETSC ERROR: -f1 Amat</div><div>[0]PETSC ERROR: -matload_block_size 1</div><div>[0]PETSC ERROR: -Pmat Pmat</div><div>[0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint@mcs.anl.gov-------<wbr>---</div><div>application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0</div></div><span class="gmail-HOEnZb"><font color="#888888"><div><br></div><div><br></div><br clear="all"><div><br></div>-- <br><div>Stefano</div>
</font></span></div>
</blockquote></div><br></div></div>