[petsc-dev] PC PBJACOBI

Mark F. Adams mark.adams at columbia.edu
Thu Aug 25 15:48:38 CDT 2011


I just got it in totalview and "diag" is 0, so I'll try to trace back and see whats up.
M
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
> 
> 




More information about the petsc-dev mailing list