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