I just got it in totalview and "diag" is 0, so I'll try to trace back and see whats up. M On Aug 25, 2011, at 4:38 PM, Barry Smith wrote:
> > Mark, > > I added support for AIJ for you, I will remove the outdated comments. > > valgrind it and let me know where the memory corruption is if you can't > figure out the bug. > > Barry > > On Aug 25, 2011, at 3:32 PM, Mark F. Adams wrote: > >> I see this in pbjacobi.c: >> /* >> Currently only implemented for baij matrices and directly access baij >> data structures. >> */ >> #include <../src/mat/impls/baij/mpi/mpibaij.h> >> #include <../src/mat/blockinvert.h> >> >> But the code looks general and I see what looks like MatInvertBlockDiagonal >> code in AIJ: >> >> But I'm getting a segv in when using and AIJ matirix. I will look into this >> but is this comment and the BAIJ include statement wrong? >> >> Mark > >