Stefano: The bug is fixed on branch hzhang/fix-ptap-seqaij https://bitbucket.org/petsc/petsc/commits/4bcca364924c81402b2f7ad375877c5546b7966b
Satish, Can you help merge this branch to next for testing? Hong I have found a bug in MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy. > Attached is a MWE (using a modified version of mat/examples/tests/ex94.c) > to reproduce it. > > As far as I can tell, the symbolic step computes an rmax for C equal to 4. > Instead, the numerical step tries to write 1+rmax at apj in that function. > Don't know what would be the proper fix. > > > [szampini@KW14141 tests]$ ./ex94 -f0 Amat -f1 Amat -Pmat Pmat > [0]PETSC ERROR: PetscTrFreeDefault() called from > MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy() line 270 in > /home/szampini/Sources/petsc/src/mat/impls/aij/seq/matptap.c > [0]PETSC ERROR: Block [id=0(16)] at address 0x22399e0 is corrupted > (probably write past end of array) > [0]PETSC ERROR: Block allocated in MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy() > line 208 in /home/szampini/Sources/petsc/src/mat/impls/aij/seq/matptap.c > [0]PETSC ERROR: --------------------- Error Message > -------------------------------------------------------------- > [0]PETSC ERROR: Memory corruption: http://www.mcs.anl.gov/petsc/ > documentation/installation.html#valgrind > [0]PETSC ERROR: Corrupted memory > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html > for trouble shooting. > [0]PETSC ERROR: Petsc Development GIT revision: v3.7.2-1249-ga13144f GIT > Date: 2016-08-25 12:12:30 +0300 > [0]PETSC ERROR: ./ex94 on a debug named KW14141 by szampini Mon Sep 26 > 10:25:06 2016 > [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 > [0]PETSC ERROR: #1 PetscTrFreeDefault() line 298 in > /home/szampini/Sources/petsc/src/sys/memory/mtr.c > [0]PETSC ERROR: #2 MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy() line 270 in > /home/szampini/Sources/petsc/src/mat/impls/aij/seq/matptap.c > [0]PETSC ERROR: #3 MatPtAP_SeqAIJ_SeqAIJ() line 42 in > /home/szampini/Sources/petsc/src/mat/impls/aij/seq/matptap.c > [0]PETSC ERROR: #4 MatPtAP() line 9257 in /home/szampini/Sources/petsc/ > src/mat/interface/matrix.c > [0]PETSC ERROR: #5 main() line 259 in /home/szampini/Sources/petsc/ > src/mat/examples/tests/ex94.c > [0]PETSC ERROR: PETSc Option Table entries: > [0]PETSC ERROR: -f0 Amat > [0]PETSC ERROR: -f1 Amat > [0]PETSC ERROR: -matload_block_size 1 > [0]PETSC ERROR: -Pmat Pmat > [0]PETSC ERROR: ----------------End of Error Message -------send entire > error message to [email protected] > application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0 > > > > > -- > Stefano >
