Chris, just a note,
Perhaps I missed something here but do something similar to you, eg
overlapping partitions, and PETSc is setup very well to be intrusive
in your code (I sometimes write a little 'addvalue' wrapper) and
assemble your FE matrices directly into a global matrix. You use the
On Oct 20, 2009, at 10:13 PM, Chris Kees wrote:
> Thanks. Just to make sure I'm following, when I create the matrix I
> do:
>
> MatCreate(PETSC_COMM_WORLD,&self->m);
> MatSetSizes(self->m,n,n,N,N);
> MatSetType(self->m,MATMPIAIJ);
> MatSetFromOptions(self->m);
> MatSeqA
gov/pipermail/petsc-dev/attachments/20091021/770a2325/attachment.html>