Jed, I assume you mean src/ksp/ksp/examples/tutorials/ex2 I have not yet been able to duplicate the error with this example. Do you know of any other test cases that might be appropriate?
My PyLith simulation with -malloc and -malloc_debug flags chokes on line 933 of src/mat/impls/baij/mpi/baijjob.c with a EXEC_BAD_ACCESS. 933 for (l=0; l<nzA; l++) cols[imark+l] = cstart + cworkA[l]; When the memory error happens at this line, the debugger shows l = 1400100 imark = 0 cstart = 0 cworkA[l] = 0 cols[imark+l] = -17958193 The address of the bad memory access matches the address of cols[imark+l]. Brad On 11/4/11 8:39 PM, Jed Brown wrote: > On Fri, Nov 4, 2011 at 21:30, Brad Aagaard <baagaard at usgs.gov > <mailto:baagaard at usgs.gov>> wrote: > > Barry, > > Sorry for the hg goof. Here is the changeset that I think is related > to the problem. > > http://petsc.cs.iit.edu/petsc/__petsc-dev/rev/7782710e3f96 > <http://petsc.cs.iit.edu/petsc/petsc-dev/rev/7782710e3f96> > > > If bs=1, then this doesn't change anything. As far as I know, it is not > possible to create an SBAIJ(1) with map->bs != 1. Can you find a way for > me to reproduce this? I see no problem with this, for example: > > mpiexec -n 4 ./ex2 -mat_type sbaij -pc_type asm
