[petsc-dev] PC PBJACOBI

Barry Smith bsmith at mcs.anl.gov
Thu Aug 25 15:38:30 CDT 2011


   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