On Aug 25, 2011, at 5:32 PM, Mark F. Adams wrote: > 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).
petscmat.h is the public include file. matimpl.h is the private include file for data structures common to all matrices, then things like mpibaij.h are private to a single matrix type (or small number). Barry > > 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 >> >> >
