I tried getting rid of this code: #include <../src/mat/impls/baij/mpi/mpibaij.h>
but got errors when I tried to look into a Mat object. Is there a generic include to get a mat object? (I tried petscmat.h). Mark 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 > >