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
