Miorelli: Fixed https://bitbucket.org/petsc/petsc/commits/a560ef987c2a1e86047c71b5951614168aab22f9
After our regression tests, I'll merge it to petsc-release. Thanks for report it! Hong Dear All, > > > > I think there may be a bug in PETSc 3.7.0, function > MatDuplicate_MPIAIJ_MatPtAP. > > Whenever I call MatDuplicate on this type of matrix I'm getting SIGSEGV. > > > > The bug is easily reproduced with the attached file. This is just ex2f > from KSP, where I added two additional lines to create an A^T A matrix > which is then duplicated. > > > > call MatTransposeMatMult(A, A, MAT_INITIAL_MATRIX, 1.d0, AtA, ierr) > > call MatDuplicate(AtA,MAT_COPY_VALUES,AtA2,ierr) > > > > > > [1]PETSC ERROR: > ------------------------------------------------------------------------ > > [1]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, > probably memory access out of range > > [1]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger > > [1]PETSC ERROR: or see > http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind > > [1]PETSC ERROR: or try http://valgrind.org on GNU/linux and Apple Mac OS > X to find memory corruption errors > > [1]PETSC ERROR: likely location of problem given in stack below > > [1]PETSC ERROR: --------------------- Stack Frames > ------------------------------------ > > [0]PETSC ERROR: Note: The EXACT line numbers in the stack are not > available, > > [0]PETSC ERROR: INSTEAD the line number of the start of the function > > [0]PETSC ERROR: is given. > > [1]PETSC ERROR: [1] MatDuplicate_MPIAIJ_MatPtAP line 78 > /state/std2/FEMI/PETSc/petsc-3.7.0/src/mat/impls/aij/mpi/mpiptap.c > > [1]PETSC ERROR: [1] MatDuplicate line 4324 > /state/std2/FEMI/PETSc/petsc-3.7.0/src/mat/interface/matrix.c > > > > > > Thanks, > > > > Federico > > > > > > *______* *______* *______* > > Federico Miorelli > > > > Senior R&D Geophysicist > > *Subsurface Imaging - General Geophysics **Italy* > > > > CGG Electromagnetics (Italy) Srl > > > *This email and any accompanying attachments are confidential. If you > received this email by mistake, please delete it from your system. Any > review, disclosure, copying, distribution, or use of the email by others is > strictly prohibited.* >
