<div dir="ltr"><div>I have found a bug in  MatPtAPNumeric_SeqAIJ_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_SparseAxpy() line 270 in /home/szampini/Sources/petsc/src/mat/impls/aij/seq/matptap.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_SparseAxpy() line 208 in /home/szampini/Sources/petsc/src/mat/impls/aij/seq/matptap.c</div><div>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div><div>[0]PETSC ERROR: Memory corruption: <a href="http://www.mcs.anl.gov/petsc/documentation/installation.html#valgrind">http://www.mcs.anl.gov/petsc/documentation/installation.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">http://www.mcs.anl.gov/petsc/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/src/sys/memory/mtr.c</div><div>[0]PETSC ERROR: #2 MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy() line 270 in /home/szampini/Sources/petsc/src/mat/impls/aij/seq/matptap.c</div><div>[0]PETSC ERROR: #3 MatPtAP_SeqAIJ_SeqAIJ() line 42 in /home/szampini/Sources/petsc/src/mat/impls/aij/seq/matptap.c</div><div>[0]PETSC ERROR: #4 MatPtAP() line 9257 in /home/szampini/Sources/petsc/src/mat/interface/matrix.c</div><div>[0]PETSC ERROR: #5 main() line 259 in /home/szampini/Sources/petsc/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----------</div><div>application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0</div></div><div><br></div><div><br></div><br clear="all"><div><br></div>-- <br><div class="gmail_signature">Stefano</div>
</div>