Mixed finite element discretization with petsc?

Jed Brown jed at 59A2.org
Mon Oct 5 02:47:51 CDT 2009


On Mon, Oct 5, 2009 at 09:39, Thomas Witkowski
<thomas.witkowski at tu-dresden.de> wrote:
> From my point the main problem is that I cannot dived
> the overall matrix row by row to the processes.

Why not?

The usual thing is to partition your mesh and then resolve ownership
of the interface nodes.  You will insert element matrices where some
of the nodes are not owned (i.e. set values for rows that are now
owned) but PETSc handles this so it is transparent to you.

Jed


More information about the petsc-users mailing list