[petsc-dev] PC PBJACOBI

[Mark F. Adams]

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
> 
>