Do you have examples of doing unstructured grid finite element computations (FEM) with PETSc?
There are at least two ways to write a finite element code using PETSc 1) use the Sieve construct in PETSc, this is a high level approach that uses a small number of abstractions to help you manage distributing the grid data structures and computing the elements into the matrices. 2) manage the grid data structure yourself and use PETSc IS and VecScatter to communicate the required ghost point communication.